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ABSTRACT 

We present the results of an investigation into the X-ray properties of radio-intermediate and radio- 
loud quasars (RIQs and RLQs, respectively). We combine large, modern optical (e.g., SDSS) and 
radio (e.g., FIRST) surveys with archival X-ray data from Chandra, XMM-Newton, and ROSAT 
to generate an optically selected sample that includes 188 RIQs and 603 RLQs. This sample is 
constructed independently of X-ray properties but has a high X-ray detection rate (85%); it provides 
broad and dense coverage of the l — z plane, including at high redshifts (22% of objects have z = 2 — 5), 
and it extends to high radio-loudness values (33% of objects have R* — 3 — 5, using logarithmic units). 
We measure the "excess" X-ray luminosity of RIQs and RLQs relative to radio-quiet quasars (RQQs) 
as a function of radio loudness and luminosity, and parameterize the X-ray luminosity of RIQs and 
RLQs both as a function of optical/UV luminosity and also as a joint function of optical/UV and 
radio luminosity. RIQs are only modestly X-ray bright relative to RQQs; it is only at high values of 
radio-loudness (R* >3.5) and radio luminosity that RLQs become strongly X-ray bright. We find no 
evidence for evolution in the X-ray properties of RIQs and RLQs with redshift (implying jet-linked 
IC/CMB emission does not contribute substantially to the nuclear X-ray continuum). Finally, we 
consider a model in which the nuclear X-ray emission contains both disk/corona-linked and jet-linked 
components and demonstrate that the X-ray jet-linked emission is likely beamed but to a lesser degree 
than applies to the radio jet. This model is used to investigate the increasing dominance of jet-linked 
X-ray emission at low inclinations. 
Subject headings: galaxies — active: quasars — general 



1. INTRODUCTION 

1.1. Radio-loud and radio-intermediate quasars 

Quasar emission is believed to result largely from ac- 
cretion onto a supermassive black hole (e.g., Lynden-Bcll 
1969). The bulk of the optical/UV continuum in radio- 
quiet quasars (RQQs) is associated with quasi-thermal 
emission originating in the accretion disk (e.g., Shields 
1978), while the X-ray emission is postulated to arise 
from Compton upscattering of disk photons occuring in 
a hot "corona" (e.g., Haardt & Maraschi 1991). This 
scenario leads naturally to a correlation between opti- 
cal/UV and X-ray luminosity. Extensive studies of RQQs 
(e.g., Avni & Tananbaum 1986; Strateva et al. 2005; 
Steffen et al. 2006; Just et al. 2007; Kelly et al. 2007; 
Green et al. 2009) have found that the optical/UV-to- 
X-ray spectral slope steepens (in the sense that objects 
become relatively less X-ray luminous) with increasing 
optical/UV luminosity. Intriguingly, most studies (see 
above references) find that there does not appear to be 
significant evolution with redshift in the spectral energy 
distributions of RQQs; despite the strong evolution in 
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the space density of quasars, these studies generally find 
that RQQs in the early universe appear to have sim- 
ilar optical/UV-to-X-ray spectral slopes to their local 
analogs. 

Radio-loud quasars (RLQs) are often defined to be the 
subset of quasars with a radio-loudness parameter sat- 
isfying R* > 1, where R* is the logarithmic ratio of 
monochromatic luminosities (with units of erg s _1 Hz -1 ) 
measured at (rest-frame) 5 GHz and 2500 A (e.g., Stocke 
et al. 1992)0; RQQs must minimally satisfy R* < 1 and 
often are found to have R* < 0. RLQs comprise ^10% of 
quasars, with this fraction apparently varying with both 
luminosity and redshift (e.g., Jiang et al. 2007 and refer- 
ences therein). The definitive physical trigger for radio- 
loudness remains elusive, but RLQs generally have more 
massive central black holes than RQQs (e.g., Laor 2000; 
Metcalf & Magliocchetti 2006; also, Shankar et al. 2010 
find this to be redshift-dependent), and it has also been 
suggested that RLQs host more rapidly spinning black 
holes than do RQQs (e.g., Wilson & Colbert 1995; Meier 
2001; but see also Garofalo et al. 2010). RLQs and 
RQQs are typically treated as distinct populations, in 
part due to the apparent relative scarcity of objects with 
R* w 1. The appropriateness of this canonical separa- 
tion has been questioned due to the discovery of numer- 
ous quasars of intermediate radio-loudness (e.g., White 
et al. 2000; Cirasuolo et al. 2003), which may outnumber 
RLQs (e.g., White et al. 2007), but there does appear 
to be a genuine bimodality of R* allowing fairly objec- 

7 Some authors measure at 4400 A instead, following Kellerman 
ct al. (1989); this method results in only a minor change (~ 0.1) 
in calculated R* values. Note that many authors prefer to define 
R* in linear units. 
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tive distinction between RQQs and RLQs (e.g., Ivezic 
et al. 2004; Zamfir et al. 2008). It should be noted 
that RQQs are not necessarily radio-silent; for example, 
Miller et al. (1993) found the radio emission from radio- 
detected RQQs to be dominated by a starburst-linked 
component p and interpreted radio- intermediate quasars 
(RIQs) as being RQQs in which a low-power and mildly 
relativistic jet is viewed at low inclinations (see also, e.g., 
Falcke et al. 1996; Zamfir et al. 2008). 

The observed properties of RLQs and their likely par- 
ent population of radio galaxies are dependent upon ori- 
entation to the observer's line of sight (e.g., Barthel 
1989; Urry & Padovani 1995). As with RQQs (e.g., 
Antonucci 1993), there is believed to be an obscuring 
"torus" present in RLQs that blocks the broad-line re- 
gion from view at large inclinations, but RLQs are fur- 
ther complicated by significant non-isotropic jet emis- 
sion. Radio jets have been measured to have relativistic 
bulk velocities on parsec scales from multi-epoch high- 
resolution radio imaging of moving knots in the inner 
jet, and the frequent lack of a detectable counterjet is 
consistent with Doppler beaming (e.g., Worrall & Birkin- 
shaw 2006 and references therein). RLQ jets terminate 
in hotspots within lobes, for which the velocities are typ- 
ically non-relativistic (e.g., Scheuer 1995) and so this 
emission is relatively isotropic. Both the ratio of radio 
core-to-lobe flux and the ratio of core radio-to-optical 
luminosities are observed to depend upon orientation, 
and these results suggest that both the lobe emission^ 
and the optical continuum in RLQs are correlated with 
intrinsic unbeamcd jet power (e.g., Wills & Brotherton 
1995). The luminosities of narrow emission lines appear 
to correlate directly with jet power, with the link plausi- 
bly coming from a mutual underlying dependence upon 
accretion rate and/or black-hole spin (e.g., Rawlings & 
Saunders 1991; Willott et al. 1999). Various unification 
models (e.g., Urry & Padovani 1995; Jackson & Wall 
1999) link narrow-line radio galaxies, broad-line radio 
galaxies, RLQs, and blazars by decreasing inclination. 
Our focus in the present study is restricted to broad-line 
quasars, but our results are of potential relevance to ra- 
dio galaxies and blazars in the context of such unification 
schemes. 

X-ray studies of RLQs strongly suggest that the nu- 
clear X-ray emission contains a significant jet-linked com- 
ponent. Zamorani et al. (1981) discovered that RLQs 
are more X-ray luminous than are RQQs with compara- 
ble optical/UV luminosities, by a typical factor of about 
three. Worrall et al. (1987) used Einstein data to show 
that the relative X-ray brightness is greater for RLQs 
that are more radio-luminous or have flatter radio spec- 
tra, and found no evidence for redshift evolution out 
to z ~ 3.5 in the properties of RLQs. Wilkes & Elvis 
(1987) and Shastri et al. (1993) uncovered X-ray spec- 
tral flattening with increasing radio loudness and ra- 
dio core dominance in samples of quasars observed with 

8 Weak radio emission from RQQs has also been suggested to be 
generated within magnetically heated coronae (Laor & Behar 2008) 
or slow and dense disk winds (Blundell & Kuncic 2007). 

9 The scatter within the correlation of core-to-lobe flux ratio to 
inclination is a factor of ~ 5—10 for a given inclination (Figure la of 
Wills &; Brotherton 1995); this scatter may reflect environmental 
effects, which can be sufficient to induce lobe asymmetries (e.g., 
Mackay et al. 1971; Gopal-Krishna & Wiita 2000). 



Einstein. Brinkmann et al. (1997) investigated a large 
sample of ROSAT-detected RLQs and found that both 
lobe-dominated and core-dominated RLQs show X-ray 
luminosity correlated with core radio luminosity, with 
the X-ray luminosity of core-dominated RLQs increas- 
ing more rapidly with increasing core radio luminosity 
(e.g., their Figure 15). It is also noteworthy in the con- 
text of unification schemes that FR II (see Fanaroff & 
Riley 1974 for description of the FR I and FR II classes) 
radio-galaxy X-ray spectra typically show both a domi- 
nant absorbed and a weaker unabsorbed component, ap- 
parently linked with the disk/corona and jet, respectively 
(e.g., Evans et al. 2006; Hardcastle et al. 2009). 

1.2. Aims of this work 

Recent wide-angle, overlapping surveys in the radio 
(e.g., Faint Images of the Radio Sky at Twenty-cm, or 
FIRST; Becker et al. 1995) and optical (e.g., the Sloan 
Digital Sky Survey, or SDSS; York et al. 2000) may be 
matched (e.g., Ivezic et al. 2002) to enable the selection 
of large, well-defined samples of RLQs, for which X-ray 
properties may be investigated. For example, Suchkov 
et al. (2006) present a catalog of SDSS Data Release 
Four (DR4) quasars matched to FIRST sources as well 
as X-ray sources from pointed ROSAT PSPC observa- 
tions. Jester et al. (2006a) matched a subset of SDSS 
DR5 quasars to FIRST sources and X-ray sources from 
the ROSAT All Sky Survey, finding radio loudness to 
be dependent upon both optical and X-ray luminosity. 
The improved capabilities of modern X-ray observato- 
ries such as Chandra and XMM-Newton have substan- 
tially advanced understanding of RLQs. For example, 
the high angular resolution and low background of Chan- 
dra enable the routine detection of X-ray emission from 
the knots of large-scale RLQ jets (e.g., Worrall 2009 
and references therein), while the broad bandpass and 
high throughput of XMM-Newton generate high signal- 
to-noise X-ray spectra useful for quantifying differences 
between RQQs and RLQs (e.g., Page et al. 2005; Young 
et al. 2009) or radio galaxies and RLQs (e.g., Belsole et 
al. 2006). Samples of SDSS quasars with X-ray coverage 
by Chandra or XMM-Newton that include subsamples 
matched to FIRST sources are presented and discussed 
by Green et al. (2009) and Young et al. (2009). 

Guided by the results described in §1.1 (and taking- 
advantage of our large sample size, which permits finer 
categorization), we consider three categories of quasars 
in this work: RQQs, RIQs, and RLQs (rather than just 
RQQs and RLQs), where we define RIQs to consist of 
objects with l<i?*< 2; consequently, the objects we clas- 
sify as RLQs satisfy R* > 2. The goal of this study is 
to quantify the optical-to-X-ray properties of RIQs and 
RLQs and to investigate the physical origin of their X-ray 
emission. Such an effort requires (1) consistent selection 
criteria that are unbiased with respect to the X-ray prop- 
erties we wish to investigate; (2) a large sample of quasars 
spanning a broad range of radio properties and possess- 
ing sensitive X-ray coverage; (3) radio imaging capable of 
resolving extended sources (multifrequency radio cover- 
age to calculate or constrain spectral indices is also use- 
ful); (4) a high fraction of X-ray detections along with 
proper statistical consideration of X-ray limits; and (5) 
effective, broad coverage of the luminosity-redshift plane 
to include the full population of RIQs and RLQs, and 
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to avoid degeneracies in regression analysis and other 
biases. We generate a large sample of RIQs and RLQs 
with archival X-ray coverage by matching the SDSS DR5 
quasar catalog of Schneider et al. (2007) and the photo- 
metrically selected quasars from Richards et al. (2009) 
to FIRST and to high-quality observations from Chan- 
dra, XMM-Newton, or ROSAT. We supplement this pri- 
mary sample with additional RLQs observed by Ein- 
stein, high-redshift RLQs, and low-luminosity RLQs de- 
tected in deep multiwavelength surveys. The full sample 
enables accurate parameterization of X-ray luminosity 
correlations across a wide range of radio properties, no- 
tably including the previously sparsely probed but well- 
populated RIQ regime. We are also able to take advan- 
tage of recent advances in statistical methods (e.g., Kelly 
2007) in our analysis and of newly established constraints 
on jet properties (e.g., Mullin & Hardcastle 2009) in our 
modeling. In addition, our use of modern accurate cos- 
mological parameters eliminates a source of systematic 
error present in some earlier analyses of luminosity cor- 
relations. 

The outline of this paper is as follows: in §2 we de- 
scribe the selection methods used to generate our sample, 
in §3 we discuss characteristics of the RIQs and RLQs 
studied here, in §4 we compare the X-ray properties of 
RIQs and RLQs to those of RQQs, in §5 we parameterize 
X-ray luminosity in RIQs and RLQs as a joint function 
of optical/UV and radio properties, in §6 we determine 
a plausible physical model for the spectral energy distri- 
butions of RIQs and RLQs, and in §7 we summarize our 
results. We adopt a standard cosmology with Hq = 70 
km s _1 Mpc -1 , VLm = 0.3, and 17a = 0.7 (e.g., Spergel 
et al. 2007) throughout. Radio, optical/UV, and X-ray 
monochromatic luminosities £ T , £ uv , and £ x are expressed 
as logarithmic quantities with units of erg s _1 Hz -1 
(suppressed hereafter), at rest- frame 5 GHz, 2500 A, 
and 2 keV, respectively. In these units the radio loud- 
ness is R* = £ T — £ uv and the useful quantity a ox , the 
optical/UV-to- X-ray spectral slope (e.g., Avni & Tanan- 
baum 1986), is a ox = 0.384x (£ K — £ uv ). Object names are 
typically taken from the SDSS DR5 spectroscopic quasar 
catalog of Schneider et al. (2007) or from the SDSS DR6 
photometric quasar catalog of Richards et al. (2009) and 
are J2000 throughout. 

2. SAMPLE SELECTION 

Our primary sample consists of 654 optically selected 
RIQs and RLQs with SDSS/FIRST observations and 
high-quality X-ray coverage from Chandra (171), XMM- 
Newton (202), or ROSAT (281). The primary sam- 
ple is split nearly evenly between spectroscopic (312) 
and high-confidence photometric (342) quasars. Most 
(562) of the primary sample objects possess serendipi- 
tous off-axis X-ray coverage, while the remainder (92) 
were targeted in the observations used in our sample. 
The X-ray detection fraction for the primary sample is 
84%; the detection fraction for those objects with Chan- 
dra/ XMM-Newton/ ROSAT coverage is 95%/92%/70% 
(typical ROSAT observations are comparatively less sen- 
sitive and have higher background) . Two additional sam- 
ples of spectroscopic SDSS quasars are constructed to 
examine the influence of SDSS targeting selection flags: 
a "QSO/HIZ" sample that contains 200 SDSS/FIRST 
RIQs and RLQs targeted as quasars from their SDSS 



colors (this is entirely a subsample of the spectroscopic 
primary sample), and a "FIRST" sample that contains 
180 SDSS/FIRST RIQs and RLQs targeted as quasars 
based on their radio emission; there is considerable over- 
lap between these two groups, but neither contain the 
quasars targeted by SDSS as "serendipitous" that we also 
include in the primary spectroscopic sample. 

We extend the primary sample of spcctro- 
scopic/photometric RIQs and RLQs with several 
supplemental samples to increase coverage of the £ — z 
plane: 93 luminous RLQs with Einstein coverage from 
Worrall et al. (1987), 13 high-redshift RLQs studied by 
Bassett et al. (2004) and Lopez et al. (2006), and 31 low- 
luminosity RLQs selected from deep multiwavelength 
surveys (see §2.2.3 for details and references) including 
the Extended Chandra Deep Field-South (E-CDF-S) 
and the Chandra Deep Field-South (CDF-S), the 
Chandra Deep Field-North (CDF-N), and the Cosmic 
Evolution Survey (COSMOS). These supplemental sam- 
ples are combined with the primary sample to produce 
the full sample. Almost all analyses are performed 
separately on the primary sample as well as on the full 
sample. 

The full sample consists of 791 quasars with R* > 1 
(with an X-ray detection fraction of 85%), of which 188 
are RIQs with R* < 2 and 603 are RLQs with R* > 2. 
The sky coverage of the full sample is shown in Figure 1. 
Properties for objects in the primary sample are provided 
in Table 1, properties for objects in the deep- fields sample 
are provided in Table 2, and characteristics of the various 
samples are provided in Table 3. In the remainder of this 
section, we provide details about the selection of all the 
samples used throughout this paper. 

2.1. Primary sample 
2.1.1. Spectroscopic sample 

The spectroscopic sample of 312 RIQs and RLQs is 
drawn from the SDSS DR5 Quasar Catalog of Schneider 
et al. (2007). The sky area covered by DR5 spectro- 
scopic observations is 5740 deg 2 near the north Galactic 
cap (Adelman-McCarthy et al. 2007). The Schneider et 
al. (2007) quasar catalog includes quasars targeted for 
matching a variety of (often overlapping) criteria (see 
Schneider et al. 2007 and Richards et al. 2002 for details). 
Of the 77429 objects in the quasar catalog, 51577 were 
targeted based on quasar-like optical colors and have 
BEST target flags of "QSO" or "HIZ" set (see Schneider 
et al. 2007 for description of these parameters). Most of 
the remaining quasars were targeted as "serendipitous" 
objects based on possessing non-stellar optical colors. A 
smaller number of quasars were targeted based on their 
radio (FIRST) or X-ray (RASS) emission and/or were 
photometrically categorized as stars or galaxies. Quasars 
targeted as "QSO" or "HIZ" or "serendipitous" were 
considered for inclusion in the optically-selected spectro- 
scopic sample; matching to the FIRST radio survey and 
then to archival high-quality X-ray coverage provides an 
initial list of 345 such RIQs and RLQs. As described 
in Appendix A, we remove from this initial list 22 BAL 
quasars, 8 highly reddened quasars, and 3 GHz-peaked 
spectrum sources. The remaining 312 objects constitute 
our spectroscopic sample of optically-selected RIQs and 
RLQs. 
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Figure 1. The sky coordinates of objects in the full sample of 
791 RIQs and RLQs. The filled symbols are X-ray detections, 
and the open symbols are X-ray upper limits. The symbol size is 
linearly proportional to _R*; larger symbols correspond to quasars 
with greater radio loudness values. Primary sample SDSS spectro- 
scopic/photometric quasars (n = 654) are plotted in blue/cyan, 
and the supplemental samples of Einstein (n = 93) and high- 
redshift (n = 13) RLQs are plotted in green and red, respectively. 
The deep-field sample (n = 31) is indicated by the locations of 
the COSMOS, E-CDF-S, and CDF-N surveys (marked with purple 
squares and labels; the square size does not indicate the solid-angle 
coverage of these surveys). 

The radio properties of these quasars are determined 
from the 1.4 GHz FIRST survey, which has a resolution 
of ~5", a 5cr limiting flux density of ~1 mJy, and 9030 
deg 2 of sky coverage, much of which overlaps the SDSS 
area (Becker et al. 1995). Objects were retained as RIQs 
or RLQs if the summed luminosity of their constituent 
components satisfied £ T > 31.0 (motivated by, e.g., Zam- 
fir et al. 2008, who define 31 < £1.4 < 31.6 as RIQs and 
£ 1A > 31.6 as RLQs) along with R* > 1.0 (with £ T and R* 
defined as in §1.2). Although the requirements of optical 
selection and a joint SDSS/FIRST detection necessarily 
limit the completeness of our sample, the depths of the 
SDSS and FIRST surveys are well matched for detect- 
ing RIQs and RLQs. An mi=19.1 quasar (the limit for 
z < 3 candidate "QSO" SDSS spectroscopic targeting) 
with R* = 1 and a typical spectral slopes (a r = —0.5) 
would have a 1.4 GHz flux density of ~0.9 mJy, near 
the FIRST point-source detection limit. Details of the 
optical screening, radio matching, and selection of X-ray 
observations are provided in Appendix A. 

In addition to the optically-selected ("QSO/HIZ" or 
"serendipitous") spectroscopic sample of 312 RIQs and 
RLQs, we construct and separately analyze a slightly 
smaller sample of 200 "QSO/HIZ" targeted spectroscopic 
quasars, and we also generate a sample of 180 quasi- 
radio-selected "FIRST" targeted spectroscopic quasars. 
Such objects were targeted by the SDSS on the (not nec- 
essarily exclusive) basis of being likely optical counter- 
parts to FIRST radio sources. This does not constitute 
a true radio-selected sample because the SDSS retains 



Figure 2. Color-magnitude plot of RIQs and RLQs, where the 
relative color A(g—i) is the g—i color (corrected for Galactic extinc- 
tion) for a given object less the median quasar color at that redshift. 
Bluer/redder objects have negative/positive values of A(cj — i). 
The filled symbols are X-ray detections and the open symbols are 
X-ray upper limits. Larger symbols correspond to more radio-loud 
quasars. Primary sample SDSS spectroscopic/photometric quasars 
are plotted in blue/cyan, and the supplemental samples of Ein- 
stein, high-redshift, and deep-field RLQs are plotted in green, red, 
and purple, respectively The A(g — i) and mj values for the deep 
field RLQs have been calculated by tranforming UBVRI magni- 
tudes to ugriz equivalents. 

optical magnitude limits for FIRST sources and because 
lobe-dominated radio sources without a FIRST core com- 
ponent will not be targeted by SDSS as FIRST sources, 
but this approach provides a useful basis for broad com- 
parison to RIQs and RLQs targeted on the basis of their 
optical colors. There is considerable overlap between the 
categories of "QSO/HIZ" and "FIRST" targeted objects; 
163 of the 180 objects (91%) with the "FIRST" flag set 
also have either the "QSO" or "HIZ" flags set. The me- 
dian properties of these samples are similar, with median 
A(.g - i) = 0.06/0.08, median t T = 33.14/33.16, median 
R* = 2.50/2.41, and median a ox = -1.41/ - 1.43 for the 
objects selected (as "QSO/HIZ") /(as "FIRST"). These 
results do not mandate that a complete radio-selected 
sample of RIQs and RLQs would have properties consis- 
tent with those of color-selected SDSS RIQs and RLQs 
(see, e.g., McGreer et al. 2009), but they do indicate that 
there is substantial overlap in these selection methods. 

2.1.2. Photometric sample 

The photometric sample of 342 RIQs and RLQs is 
constructed from a parent population of over 1,000,000 
photometric SDSS sources identified as potential quasars 
through the nonparametric Bayesian classification con- 
ducted by Richards et al. (2009) on unresolved SDSS 
DR6 objects. The efficiency of the photometric catalog 
at excluding non-quasar contaminants within the list of 
candidates is high (for example, it is estimated at 97% 
within a subsample of ^500,000 robust UV-excess candi- 
dates; it is lower for high-redshift candidates). We con- 
sider only the most likely quasar candidates by requir- 
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ing the good flag to be > 1 (this measure is determined 
based on several metrics; see Richards et al. 2009 for de- 
tails). Our analysis requires reliable redshifts and lumi- 
nosities, so we discard those sources with more uncertain 
photometric redshifts (probability < 0.5 of lying within 
the given range). Our minimum radio loudness and ra- 
dio luminosity requirements improve efficiency still fur- 
ther, as only a small fraction of the non-quasars in the 
photometric sample would be expected to display suffi- 
cient radio emission to pass these cuts; we expect that 
non-quasar contamination in the matched photometric 
SDSS/FIRST sample is negligible. By utilizing photo- 
metrically selected quasars, it is possible to expand the 
luminosity coverage of the primary sample. Over half 
of the Richards et al. (2009) catalog consists of objects 
fainter than m; = 20.1, which is already a full magnitude 
fainter than the SDSS limit for spectroscopic targeting of 
z < 3 "QSO" candidates (SDSS spectroscopic targeting 
of "HIZ" candidates is limited to objects with m\ < 20.2). 
The classification scheme is calibrated with spectroscopi- 
cally confirmed SDSS quasars, and consequently the opti- 
cal properties of these photometrically identified quasars 
are expected to be consistent with those selected from the 
DR5 Quasar Catalog. This can be verified from Figure 2, 
which shows that the relative colors of the photometric 
RIQs and RLQs are distributed similarly to those of the 
spectroscopic sample, but that the photometric sample 
extends to fainter magnitudes. 

The matching to radio sources and determination of 
X-ray coverage for the photometric quasars is identical 
to the procedure described for the spectroscopic quasars 
in Section 2.1.1 and Appendix A. This process produces 
a candidate list of 427 photometric RIQs and RLQs. 
However, 15 of these objects have SDSS spectroscopic 
redshifts obtained on an MJD of less than 53535; these 
objects were available for inclusion in the Schneider et 
al. (2007) DR5 quasar catalog but were deliberately not 
included therein, and are therefore not appropriate for 
our sampleFI (Recall that no known DR5 spectro- 
scopic quasars are permitted in our photometric sample, 
since spectroscopic data are preferred.) After excluding 
objects rejected from the DR5 Schneider et al. (2007) 
quasar catalog, and additionally culling six GHz-pcakcd 
spectrum objects, the updated candidate list of photo- 
metric RIQs and RLQs contains 406 objects. 

We perform an additional check for those RIQs and 
RLQs with photometric redshifts of z p hot > 1.9 in or- 
der to improve further sample fidelity To our knowl- 
edge, the only significant set of systematically erroneous 
photometric redshifts within the Richards et al. (2009) 
catalog, as established via cross-checking with SDSS 
spectroscopic data, consists of a small fraction of low- 
redshift (z spcc < 1) quasars assigned photometric red- 
shifts of Zphot ~ 2. (The additional and unavoidable ef- 

10 We verified that these objects were properly excluded from 
the Schneider et al. (2007) catalog and from our study. Of 
these 15 objects, 12 (082324.75+222303.2, 100656.46+345445.1, 
101858.54+591127.8, 105829.60+013358.8, 110021.06+401928.0, 
112211.80+431649.7 121026.59+392908.6, 124141.38+344031.0, 
131106.47+003510.0, 132833.56+114520.5 160740.59+254115.7, 
162625.85+351341.4) are included as BL Lacs in the DR5 cata- 
log of Plotkin et al. (2008); the other three (030055.97-002206.5, 
123251.42+123110.9, 133925.47-002705.5) have non-QSO SDSS 
spectra (classified by the SDSS pipeline as "Unknown" or 
"Galaxy" ) . 



feet of increasing redshift uncertainty at very faint mag- 
nitudes is difficult to quantify in the absence of spectro- 
scopic coverage, and we do not consider it here.) While 
such inaccuracies are atypical, it is possible to identify 
many of the low-redshift quasars with z p hot ~ 2 through 
matching to UV observations carried out by the Galaxy 
Evolution Explorer (GALEX; Martin et al. 2005). We 
make use of both redshift-dependent color-color separa- 
tion (D. W. Hogg 2009, personal communication) and 
FUV/NUV (far/near UV; -1350-1750/1750-2750 A) 
band SDSS detection rates (Trammell et al. 2007) to find 
and discard 34 RIQs and RLQs for which the z p hot ~ 2 
photometric redshift is likely inaccurate (and retain an- 
other 9 for which GALEX data correctly indicates an in- 
accurate photometric redshift but for which we are able 
to use available spectroscopic redshifts). We conserva- 
tively also discard a further 25 objects with z p hot > 1.9 
that lack GALEX coverage. Appendix B contains details 
of the use of GALEX data to assess the accuracy of ob- 
jects with z p hot > 1-9. The updated candidate list of pho- 
tometric RIQs and RLQs contains 347 (= 406 - 34 - 25) 
objects, within which the remaining fraction with this 
type of redshift misidentification is only — 1.5%. 

We found 82 (out of 347) photometric RIQs or RLQs 
with SDSS spectroscopic redshifts obtained on an MJD 
of greater than 53535; these objects were not available 
for inclusion in the Schneider et al. (2007) DR5 quasar 
catalog, and thus provide a "blind test" of our selec- 
tion methodology above. After examination of the SDSS 
spectra, only two objecttEH did not show obvious broad 
lines. Two non-quasar object from 82 photometric candi- 
date RIQs and RLQs with SDSS spectra corresponds to 
2.4%, suggesting that non-quasar contamination of our 
optical/radio matched sample is quite low, at least for 
brighter objects. Another three objecttQ show BAL fea- 
tures. Three BAL quasars from 82 photometric RIQs 
and RLQs with SDSS spectra corresponds to 3.7%, sug- 
gesting that BAL contamination of our optical/radio 
matched sample is quite low, at least for brighter ob- 
jects. This fraction is slightly lower than the typical 
fraction of RIQs and RLQs with BALs (e.g., Shankar 
et al. 2008), perhaps because the photometric color- 
selection is less efficient for BAL quasars with redder col- 
ors (with such a tendency reinforced by the requirement 
that A(g - i) < 1). There are then 342 (= 347 - 2 - 3) 
RIQs and RLQs in the photometric sample. 

The photometric redshifts for the remaining 77 pho- 
tometric RIQs and RLQs with SDSS spectra were re- 
placed with their spectroscopic redshifts. Photomet- 
ric redshifts for an additional 25 photometric RIQs and 
RLQs were replaced with spectroscopic redshifts listed 
in the NASA/IPAC Extragalactic Database (NEEEI). 
The ratio of the absolute value of the difference be- 
tween photometric and spectroscopic redshifts to the 
spectroscopic redshift was checked to assess redshift ac- 

11 085448.87+200630.7 and 133818.26+222156.4 have "Un- 
known" SDSS pipeline spectral type and Zp^^ / z S pec = 0.415/4.517 
and 2.035/0.392, respectively, but the spectroscopic redshift is 
highly uncertain. 

12 080447.96+101523.7, 141651.49+185014.1, and 
155259.18+203107.9 appear to be (DRY) BAL quasars; since, 
as discussed in Appendix A, BAL quasars are typically X-ray 
absorbed, these three objects are not appropriate for our sample. 

13 http:/ /nedwww. ipac.caltech.edu/ 
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curacy. The 9 objects already identified as having sub- 
stantially inaccurate photometric redshifts by the prior 
process of matching to GALEX data are not included 
in this comparison. There are only 1/7/24 objects for 
which this ratio exceeds 0.8/0.2/0.1, and these sub- 
stantial/modest/tiny redshift errors are relatively ran- 
dom (the median/mean/standard-deviation of z p hot — 
z spcc = -0.02/ - 0.30/0.59 for the 24 objects with 
Izphot - 2 spcc |/z spcc > 0.1, and 0.01/ - 0.08/0.33 for 
all 93 objects compared). These 1/7/24 objects corre- 
spond to percentages of l.l%/7.5%/25.8% out of the 93 
redshifts checked; applying these percentages to the 240 
(= 342 — 93 — 9) photometric RIQs and RLQs lacking 
spectroscopic coverage suggests substantial/modest/tiny 
redshift errors in 0.8%/5.3%/18.1% of the full photomet- 
ric sample. Not only are the percentages of errors small 
and relatively random, but the impact on the luminosi- 
ties for objects with modestly incorrect photometric red- 
shifts is only ~0.2 (expressed in logarithmic units), less 
than the intrinsic scatter; this should have no appreciable 
impact on our analysis. The luminosities for the photo- 
metric RIQs and RLQs with SDSS spectra or NED z apec 
values are recalculated using the spectroscopic redshifts. 
In no case did this cause the radio luminosity of a previ- 
ously accepted RIQ or RLQ to drop below the minimum 
selection cutoff values of £ T = 31. Median properties of 
the 342 objects in the final photometric sample are pre- 
sented in Table 3. 

2.2. Supplemental samples 

We supplement the primary sample with additional 
RIQs and RLQs chosen to increase coverage of the I — z 
plane, including 93 objects previously observed by Ein- 
stein, 13 high-redshift objects (primarily targeted by 
Chandra), and 31 objects selected from deep-field sur- 
veys. The properties of the deep-field objects are pre- 
sented in Table 2, and Table 3 includes the properties 
of the supplemental samples. The selection methods for 
these additional RIQs and RLQs are by necessity not 
identical to those employed to obtain our primary sam- 
ple, but the optical colors of the supplemental sample 
RLQs are similar to those of the primary sample, as can 
be seen in Figure 2. The additional I — z coverage pro- 
vided by the supplemental samples (Figure 3) helps en- 
sure consideration of the full population of RIQs and 
RLQs and also considerably reduces the i — z degener- 
acy when performing statistical analyses below (but we 
conduct most fitting on both the full and primary-only 
samples). 

2.2.1. Einstein sample 

To increase population of the high-luminosity region 
of the I — z plane, we include the RLQs with Einstein 
observations analyzed by Worrall et al. (1987), as this 
sample includes many luminous RLQs that even today 
do not have high-quality X-ray coverage from other ob- 
servatories. Their sample of 114 RLQs includes objects 
from both the North and South celestial hemispheres and 
has an X-ray detection rate of 89%. Their sample was 
primarily radio-selected at both high and low frequen- 
cies and consequently includes a mix of compact and ex- 
tended radio sources, and their RLQs tend to have higher 
radio, optical, and X-ray luminosities (and also radio- 



loudness values) than the objects in our primary sam- 
ple. We take the radio, optical, and X-ray luminosities 
from Worrall et al. (1987) but translate their values to 
our adopted cosmology. We discard three BAL quasars, 
one compact steep spectrum source, and one GHz-peaked 
spectrum source from the Einstein sample. The relative 
colors of the 56 Einstein non-BAL quasars with SDSS 
magnitudes (of which 44 also have SDSS spectra) are 
plotted on Figure 2 and are similar to the relative colors 
of the primary sample; a Kolmogorov-Smirnov (KS) test 
gives a probability p = 0.36 that the two samples are 
not inconsistent with arising from the same underlying 
distribution. Note that 3C 273 is too bright to fit on this 
plot, and also too bright to be targeted for spectroscopy 
by the SDSS. We independently find 16 objects from the 
full Worrall et al. (1987) sample in our primary sample 
(10 spectroscopic, 6 photometric) and omit these objects 
from the Einstein sample to avoid duplication. Many of 
the other Einstein objects with SDSS coverage do not 
appear in our primary sample for various reasons (e.g., 
they were targeted by SDSS based on radio rather than 
optical properties, or they lack FIRST coverage, or most 
commonly they lack recent high-quality X-ray observa- 
tions). The Einstein sample is then made up of 93 RLQs. 
Of these, 11 were undetected by Einstein, but three of 
these are detected in shallow ROSAT observations, and 
we use these £ x measurements rather than the Einstein 
limits. 

2.2.2. Htgh-redshift sample 

To increase population of the high-redshift region of 
the i — z plane, we include the 15 high-redshift RIQs and 
RLQs tabulated by Bassett et al. (2004) and the 6 high- 
redshift RLQs observed by Lopez et al. (2006). These 
high-redshift objects were selected based on radio flux as 
well as redshift and were typically (15/21) targeted by 
Chandra in "snapshot" observations of ~5 ks (6/21 were 
observed instead by ROSAT; all these are from Bassett et 
al. 2004); all are X-ray detected. The Lopez et al. (2006) 
objects have Southern declinations and are therefore un- 
available to the SDSS/FIRST surveys. Most (11/15) of 
the Bassett et al. (2004) objects have SDSS coverage, and 
most (7/11) of these have SDSS spectra and are identi- 
fied as SDSS quasars. The relative colors of the Bassett 
et al. (2004) RLQs for which we have SDSS magnitudes 
are plotted in Figure 2. One object with A(g — i) > 1 is 
discarded and not shown{3 We independently find 7 of 
the 15 objects from Bassett et al. (2004) in our primary 
sample and for consistency we use our measurements in 
our analysis of these objects. The high-redshift sample 
then consists of 13 (= 15 — 1 — 7 + 6) objects, of which 
12 are RLQs and one is an RIQ. 

Some of the high-redshift RLQs have particularly 
large radio-loudness values, with five having R* > 3.5. 
These objects are referred to as "blazars" by Bassett et 
al. (2004) and Lopez et al. (2006) and are likely viewed 
at lower inclinations than most of our primary sample 
objects (though all of the high-redshift objects are broad- 
line quasars and not BL Lacs). The relatively large frac- 
tion of objects with extreme radio-loudness values within 

14 Additionally, the redshift for 091316.55+591921.6 (z = 5.122) 
is too high to obtain a reliable relative color and so A(g — i) has 
been set to zero for this object. 
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the high-redshift sample should not necessarily be taken 
to be representative of high-redshift RLQs. 

2.2.3. Deep-fields sample 

To increase population of the low-luminosity region of 
the i — z plane, we include RIQs and RLQs identified 
from deep-field surveys; properties of these objects are 
given in Table 2. We select 16 objects by optical color 
(of which 14 have X-ray detections) and include a further 
15 X-ray detected objects known to possess broad-line 
optical spectra. 

We apply a color selection technique which approxi- 
mates that of the primary sample when searching for 
lower luminosity RIQs and RLQs in deep surveys. Our 
general procedure is to utilize the Vanden Berk ct 
al. (2000) color-color selection method to identify poten- 
tial quasars from large catalogs of objects with U BVRI 
photometry. This set of color cuts primarily selects for 
z < 2 UV-excess objects, but also includes additional 
color cuts designed to identify potential quasars at higher 
redshift (z = 2-4). Optical catalogs for the COSMOS 
region are described in Ilbert et al. (2008, 2009); for the 
E-CDF-S they include COMBO-17 (Wolf et al. 2004) and 
MUSYC (Gawiser et al. 2006); for the CDF-N they in- 
clude the Hawaii survey (Capak et al. 2007). We con- 
verted UBVRlIB to SDSS ugriz magnitudes following 
the transformations given by Jester et al. (2005; see their 
Table 1) as calculated for z < 2.1 quasars and use the 
ugriz magnitudes to calculate A(g — i) (discarding any 
heavily reddened objects with A(g — i) > 1) and £ uv . Ac- 
curate photometry is important for calculating colors, 
luminosities, and photometric rcdshifts, and so we addi- 
tionally require mj < 24 within the Chandra Deep Fields; 
since COSMOS has shallower X-ray coverage, we require 
mi < 22.5 for this survey to maintain a reasonable X-ray 
detection fraction. These magnitude limits are factors 
of ^90 and ^20 deeper than the m; < 19.1 limit for 
SDSS targeting of z < 3 "QSO" objects. These opti- 
cal selection criteria do not discriminate with respect to 
X-ray properties. The majority of the selected deep-field 
RIQs and RLQs have spectroscopic redshifts (see Table 2 
for references), and the remainder have accurate photo- 
metric redshifts (references in Table 2) that have been 
derived including UV or IR data where available. 

The resulting optically-selected quasar candidates are 
then matched to radio catalogs, and non-radio-loud ob- 
jects are removed from further consideration. This step 
significantly improves the efficiency of the candidate list 
at excluding non-quasar contaminants. The COSMOS, 
E-CDF-S, and CDF-N surveys have highly sensitive radio 
coverage, with detection limits better than ~50 /iJy at 
1.4 GHz. The VLA 1.4 GHz radio catalogs used for the 
COSMOS, E-CDF-S and CDF-S, and CDF-N fields are 
presented in Schinnerer et al. (2007), Miller et al. (2008) 
[which includes many objects also given in Kellerman 
et al. (2008)], and Biggs & Ivison (2006), respectively. 
Luminous starburst galaxies make up an increasing frac- 
tion of the radio-source population at low radio fluxes 
and luminosities (e.g., Windhorst et al. 1985; Barger et 

15 For COSMOS only, we use available u* BjVjr + i+ magnitudes 
for color selection, then convert the more precise ■u*g+r+i+z+ to 
SDSS ugriz by subtracting median differences. For the somewhat 
brighter broad-line selected COSMOS objects (described below) we 
use SDSS magnitudes directly. 



al. 2007), and so we also impose radio-selection criteria 
designed to screen out starbursts. We require £ r > 31.0 
as for the primary sample and further impose a more 
stringent requirement of R* > 1.3 (e.g., see Figure 8 of 
Barger et al. 2007) upon deep-field candidates, thereby 
decreasing potential starburst contamination of the sam- 
ple with the tradeoff of omitting some genuine radio- 
intermediate deep-field quasars. 

We next match to X-ray catalogs, with any X-ray lim- 
its estimated from sensitivity maps. We make use of 
X-ray point-source catalogs based on Chandra observa- 
tions of the E-CDF-S, CDF-S, CDF-N, and COSMOS 
as presented in Lehmer et al. (2005), Luo et al. (2008), 
Alexander et al. (2003), and Elvis et al. (2009), respec- 
tively. Maximum Chandra effective exposure times are 
~250 ks for the E-CDF-S, ~2 Ms for the CDF-S and the 
CDF-N, and ~160 ks for COSMOS. We use a matching 
radius around the optical position of 2.5", which is large 
enough to account for joint uncertainties in position but 
sufficiently small that no spurious matches are likely (see 
above references). 

As we are interested in characterizing the fundamen- 
tal X-ray emission properties of RIQs and RLQs, it is 
helpful to identify and remove objects with heavy in- 
trinsic X-ray absorption. This is important for the low- 
luminosity deep-field sample, since the fraction of ob- 
scured AGNs is large at low X-ray luminosities and 
decreases to high X-ray luminosities (from w 80% for 
a 2-10 keV luminosity of 10 42 erg s" 1 to w 20% at 
10 45 erg s _1 ; e.g., Hasinger 2008, and see also discussion 
in Brandt & Alexander 2010). Removing objects with 
strong optical reddening (as we do for both the primary 
and the deep-field samples) can also remove many ob- 
jects with X-ray absorption, but for the deep-field sample 
we take the additional step of considering X-ray spectral 
shape (but not X-ray luminosity) as a selection crite- 
rion, as measured using the X-ray hardness ratio [de- 
fined as HR= (H-S)/(H + S), where H and S are the 
2-8 keV and 0.5-2 keV counts, respectively]. We screen 
out sources that are likely absorbed by requiring that the 
hardness ratio satisfy HR < 0; this would correspond to 
a power-law slope of T ~ 1 for no intrinsic absorption, 
~ 2cr from the T ~ 1.55 typical of RLQs (e.g., Page 
et al. 2005). After application of the hardness-ratio cut, 
this color-selection method yields 16 deep-field RIQs and 
RLQs, of which 14 have X-ray detections; 9 are from the 
E-CDF-S, 4 from the CDF-N, and 3 from COSMOS. In 
Appendix C, we briefly comment on interesting aspects 
of some of these RIQs and RLQs. Most (12/16) of the 
deep-field quasars selected in this manner are RLQs with 
R* > 2. 

In addition to the color selection, we also employ opti- 
cal spectra for selection, accepting RIQs or RLQs which 
are known to have broad emission-line spectra. Optical 
spectroscopy specifically targeting X-ray sources is avail- 
able in the CDF-S (based on 1 Ms sources; Giacconi et 
al. 2002), the CDF-N (based on 2 Ms sources; Alexan- 
der et al. 2003), and COSMOS (based on XMM-Newton 
detections; Cappelluti et al. 2009), and is presented in 
Szokoly et al. (2004), Barger et al. (2003) as well as 
Trouille et al. (2008), and Trump et al. (2009), respec- 
tively. We include 15 RIQs and RLQs selected based 
on broad-line emission. Of these, one is from the CDF- 
N and two are from the CDF-S; these three objects have 



radio luminosities of l x < 31 and so were not available for 
consideration as color-selected objects, but they do sat- 
isfy the Vanden Berk et al. (2000) color-color selection 
method for quasar candidates, as well as the hardness 
ratio cut, and so would otherwise have been expected 
to be selected through this process. The remaining 12 
RIQs and RLQs are from the XMM-Newton COSMOS 
survey, which is shallower than the Chandra COSMOS 
survey but covers a wider area; many of these broad-line 
objects do not fall within the Chandra coverage and/or 
have radio luminosities of l x < 31, and so were not con- 
sidered for color selection. Again, all of these broad-line 
objects would have been identified as quasar candidates 
based on their optical colors, and all further satisfy the 
hardness ratio cut. 

In principle this manner of selecting broad-line ob- 
jects known to be X-ray sources might introduce a bias 
toward X-ray bright sources, were similarly optically- 
bright broad-line objects with X-ray limits not also con- 
sidered. However, we are unaware of any broad-line 
objects present in the extensive optical surveys of the 
E-CDF-S or CDF-N which lack X-ray detections, and 
in any event the RIQs and RLQs selected based on op- 
tical spectroscopy could generally have had lower X-ray 
counts by factors of ^10 and still been detected, suggest- 
ing that their relative X-ray brightness is not atypical for 
their optical/UV luminosities. This is supported by the 
observation that the a ox values of the broad- line selected 
objects are similar to those of the color-selected objects. 

The complete deep-field sample then consists of 31 
RIQs and RLQs, of which 29 are X-ray detected. Of 
these, 16 were selected by optical color (of which 14 have 
X-ray detections) and 15 were selected based on possess- 
ing broad-line optical spectra. 

3. SAMPLE CHARACTERISTICS 

Characteristics of the various samples are presented 
in Tables 3 and 4 and illustrated in Figures 1-6. The 
sky coverage of the full sample is indicated in Figure 1, 
a color-magnitude diagram is provided as Figure 2, the 
rcdshift distribution is plotted versus optical/UV lumi- 
nosities in Figure 3 and versus X-ray luminosities in Fig- 
ure 4, the radio characteristics of the primary sample are 
given in Figure 5, and the optical/UV-to-X-ray spectral 
slope as a function of optical/UV luminosity is shown in 
Figure 6. 

3.1. Optical/UV luminosities 

Optical/UV monochromatic luminosities for all pri- 
mary sample objects are calculated at rest-frame 
2500 A from SDSS photometric ugriz PSF magnitudes 
(corrected for Galactic extinction) via comparison to a 
rcdshifted unobscured composite quasar spectrum taken 
from Vanden Berk et al. (2001). This method accounts 
(in a statistical sense) for the flux from emission lines 
as well as typical breaks in the continuum slope. The 
error in this method due to differences in emission-line 
strength or spectral shape in a particular object is ex- 
pected to be less than the inherent uncertainty (^30%; 
e.g., Vanden Berk et al. 2004; Kaspi et al. 2007) due to 
typical quasar optical variability. Only the AGN power- 
law component is included in the luminosity; a typical 
iron-emission "bump" near 2500 A (e.g., Wills et al. 1985; 



Vanden Berk et al. 2001) is subtracted, as is the con- 
tribution from a typical RLQ host galaxy (£ uv = 28.2 
based on an old stellar population as in, e.g., Coleman 
et al. 1980). We verify the accuracy of this approach 
through comparison to ~ 100 RQQs for which Strateva et 
al. (2005) calculated luminosities by fitting SDSS spec- 
tra (after dereddening and correcting for fiber inefficien- 
cies, and also subtracting host-galaxy emission) and find 
close agreement (mean difference of £ uv — £ uv _ s = —0.03 
with a standard deviation of 0.12). By using photometric 
rather than spectroscopic data to compute optical/UV 
luminosities, we can treat the SDSS spectroscopic and 
photometric samples in a consistent fashion. 

Optical/UV luminosities for the Einstein RLQs are 
taken from Worrall et al. (1987), corrected to our adopted 
cosmology. Optical/UV luminosities for the high-redshift 
RLQs are taken from Bassett et al. (2004) and Lopez et 
al. (2006). Optical/UV luminosities for the deep field ob- 
jects are calculated from SDSS ugriz equivalents, deter- 
mined as described in §2.2.3. The full sample (Figure 3) 
achieves dense coverage of the i—z plane, a wide span (al- 
most five decades) in luminosity coverage (~2.5 decades 
at a given redshift even without considering deep-field 
RIQs and RLQs, and ^3.5 including deep-field objects), 
and coverage to z ~ 5. 

3.2. Radio luminosities and radio loudness 

Radio monochromatic luminosities are calculated at 
rest-frame 5 GHz through extrapolation of the observed 
1.4 GHz flux densities. It is desirable to treat the entire 
sample in a uniform fashion, but many sources lack multi- 
frequency radio measurements, or are multi-component 
sources with multi-frequency radio flux densities ob- 
tained at an angular resolution insufficient to distinguish 
a r for the core from any extended emission. Therefore, 
we do not use individual a T values calculated for a partic- 
ular source (see §3.3) to determine the radio luminosity 
of that source, but instead assume a radio spectral in- 
des|3 of a x = —0.9 for lobe emission and a r = —0.3 
for core emission. These values are approximately the 
mean a r for lobe-dominated and core-dominated sources 
observed within the primary sample, respectively (using 
a r calculated from low-frequency radio measurements for 
the lobe-dominated sources). In any event, alternative 
choices of a T produce only small changes in l T since it 
is only necessary to extrapolate over a small frequency 
range. The total £ T is the sum of the core and lobe com- 
ponents, and we also provide ^r iC orc in Table 1. The ra- 
dio monochromatic luminosities within the full sample 
span over four decades, with a median £ x — 32.95. The 
median radio monochromatic luminosity within the pri- 
mary sample is slightly lower, with £ T = 32.83. This 
difference chiefly reflects the high radio luminosities of 
the supplemental Einstein RLQ sample, which may be 
due to some targets being radio selected. The radio lu- 
minosities within the deep-field sample are lower, with 
a median l x = 31.50 (recall that we permit £ x < 31 for 
broad-line selected deep-field objects). 

Radio loudness in our sample is correlated with radio 
luminosity (Figure 5), signifying the influence of beamed 

16 The radio spectral index is defined as a T = 
log (Sh/Sl)/ log (Uh/ui), with v^,v\ and S^, Si the high, low 
frequencies (e.g., 5 and 1.4 GHz) and corresponding flux densities. 
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Figure 3. Sample coverage of the £ U v — z plane, with optical/UV luminosity as measured at rest-frame 2500 A. The filled symbols are 
X-ray detections, and the open symbols are X-ray upper limits. Larger symbols correspond to more radio-loud quasars. Primary sample 
SDSS spcctroscopic/photomctric quasars are plotted in blue/cyan, and supplemental sample Einstem/high-rcdshift/dccp fields RLQs arc 
plotted in green/red/purple. The dotted lines define luminosity bins within which various properties are summarized in Table 4. 



jet emission that comes to dominate the radio emission 
measured by FIRST for objects at low inclinations (or 
with intrinsically powerful jets) but apparently does not 
similarly dominate the optical/UV luminosity (as is also 
suggested by the SDSS broad emission-line equivalent 
widths, which do not suggest significant dilution by a 
featureless jet-linked continuum for these sources). The 
median radio loudness for the full sample is R* — 2.59; it 
is notably higher for the Einstein sample (median R* = 
3.44) and lower for the deep-field objects (median R* = 
2.06, with the lowest permitted value being R* = 1.3). 

3.3. Radio spectral shapes and morphologies 

Although we do not use individually measured radio 
spectral slopes to calculate luminosity, we do distinguish 
between flat-spectrum and steep-spectrum radio sources 



when radio spectral information is available. Following 
Worrall et al. (1987), objects with a r < —0.5 are clas- 
sified as steep spectrum, while objects with a T > —0.5 
are classified as flat spectrum. Flat-spectrum RLQs are 
X-ray brighter than their steep-spectrum counterparts 
(e.g., Worrall et al. 1987), so we also conduct analy- 
ses of X-ray luminosity correlations on separate subsam- 
ples of flat-spectrum and steep-spectrum RLQs. The 
Einstein RLQs arc identified as either flat-spectrum or 
steep-spectrum in Worrall et al. (1987), and we use their 
classification; we do not identify sources from the other 
supplemental samples as either flat-spectrum or steep- 
spectrum. In light of the difficulties inherent in com- 
paring fluxes obtained at widely differing angular reso- 
lutions, spectral indices for primary sample objects with 
multifrequency radio data are computed in the observed 
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Figure 4. Sample coverage of the £ x — 2 plane, with X-ray lu- 
minosity as measured at rest-frame 2 keV. The filled symbols are 
X-ray detections, and the open symbols are X-ray upper limits. 
Larger symbols correspond to more radio-loud quasars. 

frame from the source flux densities (summed over all 
components for cases of resolved extended radio emis- 
sion) . 

To find a T for objects with multifrcqucncy radio cov- 
erage, counterparts for all primary sample objects were 
sought in the Green Bank 6 cm (5 GHz) catalog (Gregory 
et al. 1996), which has a flux-density limit of ~18 mJy 
and covers 0° < S < 75°, in the Texas 82 cm (365 MHz) 
catalog (Douglas et al. 1996), which has a flux-density 
limit of ^250 mJy and covers —35° < 5 < 70°, and in 
the Westerbork 92 cm (325 MHz) Northern Sky Survey 
(Rengelink et al. 1997), which has a flux-density limit 
of ~18 mJy and covers 6 > +30°. All of these surveys 
have significantly lower angular resolution than does the 
FIRST survey, so components resolved by FIRST may be 
blended in these catalogs. Any sources within 90" of the 
optical position are presumed to be associated with the 
quasar, unless a FIRST background source is known to 
be present within the field; for such instances the relative 
positions and fluxes have been evaluated on a case-by- 
case basis and the multi-frequency radio data discarded 
if deemed background contaminated^ 

Green Bank data are prioritized when calculating a r 
since we are most interested in assessing the relation- 
ship of X-ray emission to the radio core (rather than 
extended) emission. The uncertainty on a r is then often 
dominated by the error in the Green Bank flux mea- 
surements, which is generally ~10%. A quasar with a 
typical FIRST radio flux and a borderline a r = —0.5 
radio spectral slope would have an uncertainty on a r of 

17 Because we use the higher-resolution FIRST survey to screen 
for potential background sources, the effective matching radius is 
less than 90". The false-match probability is low even without 
FIRST screening: artificially shifting the declination by one degree 
and rematching returns 5.6% as many matches within 90". How- 
ever, most of these shifted matches are sufficiently distant from 
the core that they would be subject to screening as potential back- 
ground objects as gauged by FIRST data. 



Figure 5. Radio loudness (R*) versus radio luminosity for the 654 
RIQs and RLQs in the SDSS primary sample. Symbol color iden- 
tifies objects as flat-spectrum core-dominated (FSCD; cyan), flat- 
spectrum lobe-dominated (FSLD; blue), unmeasured radio spec- 
trum core-dominated (NSCD; yellow), unmeasured radio spectrum 
lobe-dominated (NSLD; green), steep-spectrum core-dominated 
(SSCD; purple), and steep-spectrum lobe-dominated (SSLD; red). 
Symbols are slightly larger for lobe-dominated objects and open for 
X-ray limits. Radio loudness and luminosity are correlated within 
this sample, albeit with large scatter. Our sample shows the well- 
known tendencies for the most radio-loud objects to be primarily 
core-dominated and flat-spectrum, and for lobe-dominated objects 
to be mostly steep-spectrum; see also discussion in §3.2 and §3.3. 

±0.08. There are 43 objects with -0.58 < a s < -0.42, 
which is 11% of the total number of objects with deter- 
mined a r values, so for most objects the classification 
as flat or steep spectrum is secure. Objects with both 
Green Bank and low-frequency measurements were con- 
sidered more closely. Sources with concave spectra are 
presumed to arise from the emerging dominance of a flat- 
spectrum core and are retained. However, as described in 
Appendix A, the 9 objects with convex spectra are iden- 
tified as possible GPS sources, which have properties not 
shared by RLQs in general and are thus not included in 
the primary sample. 

All primary sample objects are also classified as either 
core-dominated or lobe-dominated, with core-dominated 
objects having a core radio monochromatic luminosity 
(at 5 GHz) greater than half the total radio luminos- 
ityE3 As a consequence of performing radio selection at 
an observed frequency of 1.4 GHz (rather than at, for ex- 
ample, a lower frequency such as 159 MHz or 178 MHz, 
at which the 3C and 4C surveys, respectively, were car- 
ried out), core-dominated, likely low- inclination sources 
are over-represented relative to their presumed parent ra- 
dio population. (Recall, however, that extremely beamed 
objects are mostly excluded from our sample by the op- 
tical selection criteria; for example, objects with feature- 

18 Similar but slightly differing definitions of core-dominated are 
sometimes used; for example, Wills et al. (1995) define core/lobe- 
dominated sources to be those for which the ratio of core/lobe 
emission at rest-frame 6 cm is greater/less than unity. 
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less optical spectra are not included in the SDSS DR5 
quasar catalog.) The primary sample lobe-dominated 
RIQs and RLQs are more likely to have steep radio spec- 
tra, whereas the core-dominated objects are more likely 
to have flat radio spectra (Figure 5). The primary sam- 
ple RLQs with particularly large radio-loudness values 
(R* > 3.5) tend to be core-dominated and have flat radio 
spectra. Such well-known trends are likely in large part 
a consequence of core small-scale jet emission gradually 
overwhelming (steep-spectrum) lobe emission as inclina- 
tion decreases (e.g., Orr & Browne 1982; see also §6). 

3.4. X-ray luminosities 

X-ray counts were measured for all primary sample ob- 
jects using the IDL aper task. X-ray images for objects 
observed with ROSAT or XMM-Newton were down- 
loaded from HEAS ARC0 along with exposure maps and, 
for ROSAT, background images. For objects observed 
with Chandra, the CIAO task dmcopy was used to pro- 
duce full-band images from the pipeline-processed level 
2 event files. Source counts were extracted from ^90% 
encircled-energy aperturesJH3 The energy coverage of 
the utilized images is 0.4-2.4 keV for ROSAT, 0.5-8 keV 
for Chandra, and 0.2-12 keV for XMM-Newton. Back- 
ground counts were determined from the provided back- 
ground image for ROSAT observations and as the me- 
dian of 8 nearby non-overlapping apertures for Chandra 
and XMM-Newton observations. 

Source detection is evaluated by comparison of the ob- 
served aperture counts to the 95% confidence upper limit 
for background alone. Where the number of background 
counts is less than 10 we use the Bayesian formalism of 
Kraft et al. (1991) to determine the limit; else, we use 
Equation 9 from Gehrels (1986). If the aperture counts 
exceed the 95% confidence upper limit we consider the 
source detected and calculate the net counts by subtract- 
ing the background from the aperture counts and then 
dividing by the encircled-energy fraction; else, the source 
is undetected and the 95% confidence upper limit is used. 
Counts were converted to count rates using the furnished 
or generated exposure maps. Count rates were converted 
to observed- frame 2 keV flux densities with PIMMS[H3 in 
all cases assuming Galactic absorption and a power-law 
spectrurrF^I with Y = 1.5 (alternate reasonable choices 
for r have only a few percent impact upon the calcu- 
lated X-ray fluxes), which were then used to determine 
bandpass-corrected rest-frame 2 keV monochromatic lu- 
minosities. 

Where available, data from the Chandra Source Cat- 
alog (Evans et al. 2008) or the XMM-Newton Serendip- 
itous Source Catalog (Watson et al. 2009) were used in 
preference to our own measurements. The X-ray lumi- 
nosities were calculated from catalog broad-band fluxes 



http 



19 High Energy Astrophysics Science Archive Research Center: 

: / /heasarc. gsfc.nasa.gov/ 
A<3 The ~90% encircled-energy apertures were determined 
for Chandra as in Luo et al. (2008), for XM M-Newton as in 
http: / / xmm. esa. int / external / xmm.user .support / documentation/ 
uhb/nodel8.html, and tor HUbA'l as m http://www.mpe.mpg.de/ 
xray / wave /rosat /doc/ tech_reports / rosat_pst.ps. 

21 http:/ /cxc. harvard.edu/toolkit/pimms.jsp 

22 RLQ X-ray spectra can generally be fitted with a single power 
law (e.g., Belsole et al. 2006), although some high-redshift RLQs 
do show evidence of intrinsic absorption (e.g., Cappi et al. 1997; 
Yuan et al. 2006) or spectral curvature (e.g., Tavecchio et al. 2007). 
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Figure 6. Optical/UV-to- X-ray spectral slope a ox as a function of 
< UT ; symbols are coded as in Figure 3. The solid/dotted line is the 
best-fit linear relation for RQQs from Just et al. (2007)/Stcffcn ct 
al. (2006), while the dashed lines show the 25th and 75th percentiles 
for RQQs. Our sample shows the well-known tendency for RLQs 
to be X-ray bright relative to RQQs of comparable optical/UV 
luminosities. 



as given at 0.5-7 keV for Chandra and 0.5-4.5 keV for 
XMM-Newton, with matching performed using radii of 
3" and 10", respectively. It is not possible for this project 
to rely exclusively on catalogs, however, as large "blind 
search" source catalogs must utilize more conservative 
detection thresholds. In addition, it would be difficult 
to determine accurate upper limits based on a cata- 
log non-detection. Our calculated luminosities are in 
good agreement with those derived from catalog data 
(median/mean/standard deviation of the difference in 
4 of 0.037/0.048/0.201 and 0.043/0.079/0.186 for Chan- 
dra and XMM-Newton, respectively), indicating that our 
measurements (detections or upper limits) for those ob- 
jects not included in catalogs should also be accurate. 

The X-ray data for the primary sample of 654 RIQs 
and RLQs are based for 171/202/281 objects on Chan- 
dra/ XMM-Newton/ ROSAT observations. The X-ray de- 
tection rates for Chandra/ XMM-Newton/ ROSAT are 
95%/92%/70%. X-ray data for 115 objects detected with 
Chandra were taken from the Chandra Source Catalog, 
and X-ray data for 160 objects detected with XMM- 
Newton were taken from the XMM-Newton Serendipi- 
tous Source Catalog. In a few cases these source cata- 
logs provide total X-ray coverage deeper than our single- 
observation photometry, and our sample is improved by 
making use of these additional data. No objects are dis- 
carded based on non-inclusion in source catalogs. Many 
of the sources which are not included in the utilized 
source catalogs but are detected by our photometry are 
associated with observations not included in the source 
catalogs, either due to the observation date falling out- 
side of the range covered by the source catalogs or some 
observation parameter (e.g., subarray type) not satisfy- 
ing the requirements for inclusion. 
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Figure 7. X-ray "excess" in RIQs and RLQs relative to comparable RQQs plotted as a function of radio loudness (top) and radio 
luminosity (bottom). The solid line corresponds to a linear multiplicative ratio of one and the dashed line to a ratio of three. The black 
crosses show median values within bins; the arms show the 25th-75th percentile range in each bin (values given in Table 5). Mean values 
are similar to medians and the error on the mean is much smaller than the interquartile range shown. Symbols are coded as in Figure 3. 



Almost all of the ROSAT coverage is serendipitous; 
only 6/281 (2.1%) sources are found within 3' of the 
aimpoint of the best ROSAT observation. The SDSS 
and FIRST surveys served as a source of targets for 
many XMM-Newton and Chandra observations; 34/202 
(16.8%) sources are found within 2' of the aimpoint of 
the best XMM-Newton observation, and 58/171 (33.9%) 
sources are found within 1' of the aimpoint of the best 
Chandra observation. Some of these targeted objects also 
have serendipitous coverage and would be in our sample 
even absent the targeted observations: four of the objects 
targeted by XMM-Newton have off-axis Chandra cover- 
age and two of the objects targeted by Chandra have off- 
axis ROSAT coverage. In total, then, 92/654 (14.1%) of 
the RIQs and RLQs in our primary sample were targeted 
for X-ray observations (of which 91/92 have X-ray de- 
tections) . The analysis of luminosity correlations below 
was carried out including targeted objects to increase the 
size of the primary sample, but results are also provided 



for an "Off-axis" sample of sources with only serendipi- 
tous X-ray coverage and for the "Targeted" sample ex- 
clusively. Differences between the properties of the "Off- 
axis" and "Targeted" subsamples are discussed in §5.4.2. 

There may be cases in which the extraction region used 
to calculate the X-ray luminosity includes both nuclear 
and kpc-scale X-ray jet emission. Many X-ray jets have 
been discovered by Chandra (e.g., Sambruna et al. 2004; 
Marshall et al. 2005), and some would not be resolvable 
with ROSAT or even XMM-Newton. Additionally, in 
some instances the inner knot(s) of an X-ray jet might 
lie inside the Chandra extraction region. Generally even 
the inner knots of RLQs with X-ray jets are observed 
to be only a few percent as bright as the X-ray coreF^I 

23 The X-ray core itself may contain an unresolved component 
of X-ray emission linked to the pc-scale radio jet, of course, but we 
desire such a component be included for our analysis and in any 
case it would be impossible to exclude from simple photometry. 
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Figure 8. Histograms of 4c — ^ x ,rqq (defined as in §4) plot- 
ted in (increasing downward) bins of radio-loudness (left) and 
radio luminosity (right); the bins are identical to those in Ta- 
ble 5. The filled histograms contain X-ray-detected objects 
and the open histograms include upper limits. The vertical 
dashed lines correspond to linear multiplicative factors of 1 
and 3 for the ratio of the X-ray luminosity of RIQs and RLQs 
to comparable RQQs. The solid/dotted curves are Gaussians 
fit to detected/all objects to investigate log-normality; there 
appears to be a tail of objects with high X-ray luminosities 
in the maximum radio-loudness and luminosity bins. 

(e.g., Marshall et al. 2005), which would not significantly 
change the calculated £ x values. The XJET catalo eri 
provides a useful listing of ^100 objects with known ex- 
tended X-ray jet or lobe emission. Out of the primary 
sample of 654 RIQs and RLQs, only 15 are listed in the 
XJET catalog. We examined X-ray images of all 15 ob- 
jects and in no case did it appear likely that the extended 
X-ray emission could significantly increase the measured 
X-ray luminosity. 

4. COMPARISON OF THE X-RAY EMISSION FROM RQQS, 
RIQS, AND RLQS 

Several previous studies have noted the tendency for 
RLQs to be more X-ray bright than RQQs of compara- 
ble optical/UV luminosity (e.g., Zamorani et al. 1981; 
see also discussion and references in §1). We confirm 
that general result, and also quantify the degree to which 
RIQs and RLQs differ in X-ray brightness from RQQs as 
a function of radio loudness and luminosity. Our large 
and carefully constructed sample permits relatively fine- 
grained binning for such measurements. 

The necessity of controlling for optical/UV luminos- 
ity when comparing RIQs and RLQs to RQQs is driven 
by the well-known steepening in RQQs of optical/UV-to- 
X-ray spectral slope with increasing optical/UV luminos- 
ity (e.g., a ox = —0.140 x £ uv + 2.705 as given in Equation 
3 of Just et al. 2007; other studies as listed in §1.1 typ- 
ically find similar results). As RLQs are X-ray brighter 

24 http://hca-www.harvard.edu/XJET/ 



than RQQs at a given £ uv , they have less negative values 
of a ox . The median optical/UV-to-X-ray spectral slope 
for the full sample of RIQs and RLQs is a ox = —1.40, 
and for RIQs/RLQs separately it is a ox = —1.50/ — 1.37. 
The sample of RQQs we use for comparison (from Steffen 
et al. 2006) has a median a ox — —1.52. 

The a ox and £ uv values for the complete sample of RIQs 
and RLQs are plotted in Figure 6, along with the RQQ 
aox(^uv) relations from Just et al. (2007) and Steffen et 
al. (2006), and the 25th and 75th percentiles for a ox for 
RQQs within each £ uv decade (from Table 5 of Steffen et 
al. 2006). The tendency of RQQs at lower luminosities to 
fall below the best-fit a ox (^ uv ) relation appears to be gen- 
uine (e.g., Steffen et al. 2006; Maoz et al. 2007); a larger 
sample is required to determine whether a qualitatively 
similar effect may apply to RIQs and RLQs. Like RQQs, 
RIQs and RLQs also show an anti-correlation between 
a ox and £ uv (a Spearman test gives p = —0.31, with a 
probability less than p = 5xl0 -5 that no correlation 
is present), albeit with a systematic offset toward less 
negative values of a ox . Figure 6 also indicates that the 
degree to which RIQs and RLQs are brighter in X-rays 
than comparable RQQs is dependent upon radio loudness 
(the larger points, representing more radio-loud objects, 
are generally further above the RQQ a ox (£ uv ) relation); 
we now quantify this dependence. 

The "excess" X-ray luminosity from RLQs may be de- 
fined as £ x — ^x,rqq, where we take the expected X-ray 
luminosity for RQQs to be 4,RQQ = 0.709 x £ uv + 4.822 
(Equation 7 from Just et al. 2007) and use £ uv as mea- 
sured for RIQs and RLQsH Fi gure 7 shows excess X-ray 
luminosity for the complete sample of RIQs and RLQs 
as a function of radio loudness and of radio luminosity, 
along with median values and the 25th-75th percentile 
range within bins of R* and £ r . The median values and 
25th and 75th percentiles plotted in Figure 7 are listed in 
Table 5. The mean values are consistent with the median 
values, and the error on the mean is much less than the 
25th-75th percentile range. When expressed in linear 
units, the multiplicative factor by which the X-ray lumi- 
nosity for RIQs and RLQs exceeds that of RQQs ranges 
(25th-75th percentiles) from 0.7-2.8 for RIQs through 
the canonical ~3 for RLQs (e.g., Zamorani et al. 1981) 
to 3.4-10.7 for extremely radio-loud (R* > 3.5) RLQs, 
with a qualitatively similar increase in excess X-ray lumi- 
nosity with increasing radio luminosity. Figure 8 shows 
the distribution of £ x — £ x ,rqq (the shaded histogram is 
detected objects and the open histogram includes X-ray 
upper limits) within the same bins of R* (left) and £ r 
(right) as used to construct Table 5. The distribution of 
excess X-ray luminosity is reasonably well-characterized 
as log- normal (see overplotted Gaussians), with KS test 
probabilities p > 0.13 in all cases. There appears to be 
a tail to brighter X-ray luminosity within the highest 
radio-loudness and luminosity bins (e.g., the percentage 
of objects above 1.645cr from the mean is not 5% but 

25 This probes the relationship between radio and jet-linked 
X-ray emission, presuming the optical/UV luminosity in these 
RIQs and RLQs is disk-dominated and that the disk/corona in 
RIQs and RLQs displays a similar dependence of X-ray emission 
upon £ uv as in RQQs. The first presumption is consistent with 
the apparently undiluted equivalent widths of the broad emission 
lines in the primary sample of RIQs and RLQs; see §6 for further 
discussion of these points. 
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Figure 9. X-ray luminosity as a function of optical/UV luminosity for RIQs and RLQs (a) and for various sub-samples of RLQs (b, c, 
d). The model considered is £ x = ao +b uv x£ uv (where f x and £ uv have been normalized prior to fitting as described in §5). Filled symbols 
are X-ray detections and larger symbols are more radio-loud throughout. Solid lines show best-fit linear models. 



rather 10%/11% within the most radio- loud/luminous 
bin, calculated using the fitted values for a and /i). Sim- 
ilar results hold for the primary sample onl y (T able 6), 
with the mean values of 4 — 4,rqq consistent^] between 
the full and primary samples across R* and £ r bins. 

The excess X-ray luminosity may also be fit directly as 
a function of radio loudness or luminosity, for example 
as 4 - 4,rqq = a + b x R* or 4 - 4,rqq = a + bx£ T , 
where a and b are fitting constants. We carry out such 
a fit using the IDL code of Kelly (2007), which utilizes 
Bayesian techniques that incorporate both uncertainties 
and upper limits. The best-fit models for the full sample 
are £ x - 4 RQQ = (-0.354 ±0.050) + (0.289 ± 0.018) x R* 
and4-4,RQQ = (0.450±0.015) + (0.177±0.014)x(4- 
33.3). Excess X-ray luminosity is more strongly depen- 
dent upon radio-loudness than radio luminosity. Flat- 

26 In only 2/10 cases (the highest radio-loudness and luminos- 
ity bins) is the difference in means greater than the error; the 
mean/standard deviation of the difference in means is 0.02/0.04. 



spectrum RLQs may have excess X-ray luminosity more 
strongly correlated with radio properties than do steep- 
spectrum RLQs, with flat/steep spectrum RLQs having 
coefficients of (0.352 ± 0.039)/(0.274 ± 0.040) xR* and 
(0.204 ±0.027)/(0.088± 0.029) x (4-33.3). The large 
amount of scatter in these relations prevents productive 
consideration of more complex models, but Figures 7 and 
8 do suggest that these linear fits (to log quantities) may 
not adequately capture the apparent slow rise in X-ray 
excess at low radio-loudness or luminosity and the more 
rapid increase at higher R* or £ T values. 

5. PARAMETERIZING THE X-RAY LUMINOSITY OF RIQS 
AND RLQS 

We parameterize X-ray luminosity as a sole function 
of optical/UV luminosity and as a joint function of op- 
tical/UV and radio luminosity for various groupings of 
RIQs and RLQs, and also consider whether an addi- 
tional dependence upon redshift is required. All fit- 
ting is carried out with the IDL code of Kelly (2007); 
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Figure 10. Joint 90% (solid) and 68% (dotted) confidence ellipses 
for fitting X-ray luminosity as a (sole) function of optical/UV lu- 
minosity for various subsamples of RLQs. The model considered is 
= io + b uv x£ uv (where £ x and £ uv have been normalized prior 
to fitting as described in §5). RLQs with larger values of radio- 
loudncss or luminosity or with fiat radio spectra have greater best- 
fit intercepts than those with smaller values of radio-loudness or 
luminosity or with steep radio spectra. 

using the Astronomy Survival Analysis (ASURV) pack- 
age (e.g., Isobe & Feigelson 1990) gives similar results. 
The potential measurement errors in optical magnitudes 
and radio fluxes are generally small (see references to 
SDSS and FIRST in §1), and most objects have sufficient 
X-ray counts that the uncertainties may be assumed to 
be dominated by intrinsic random flux variability (we use 
20%/30%/40% or 0.0792/0.114/0.146 in log units for un- 
certainties on 4/4v/4, motivated by, e.g., §3.5 of Gib- 
son et al. 2008). The luminosities are normalized prior 
to fitting as £ r — 33.3, £ uv — 30.5, £ x — 27.0 (the subtracted 
values are nearly the medians for the sample of RLQs). 
Results are given in Tables 7 and 8 and illustrated in Fig- 
ures 9-17. For any given model fit the quoted parameter 
values are the median of draws from the posterior distri- 
bution and the errors are credible intervals corresponding 
to lcr. See Appendix D for additional discussion of fitting 
methodology. 

5.1. ^x(^uv) 

The first model we consider is X-ray luminosity as a 
sole function of optical/UV luminosity, such as is of- 
ten applied to RQQs (e.g., Steffen et al. 2006; Just et 
al. 2007): £ x = ao + b uv x £ uv . It is most convenient 
to treat £ x as the dependent variable for the purposes of 
fitting, given the X-ray limits (and perhaps also most ap- 
propriate for our predominantly optically-selected sam- 
ple of RIQs and RLQs, whose X-ray properties were not 
considered in selection). Since our analysis is compara- 
tive in nature, it suffices to maintain consistency and so 
we do not also calculate coefficients treating £ uv as the 
dependent variable, nor do we calculate a bisector fit. 
This approach also simplifies analysis when additional 
variables are considered. We first re-fit the RQQs from 
Steffen et al. (2006) using the same procedure we adopt 



for analyzing the RIQs and RLQs, to demonstrate the 
consistency of this method with previous work. Our re- 
sults agree with those of Steffen et al. (2006) for their 
Equation la (see also Appendix D). We also re- fit the 
luminous Einstein RLQs from Worrall et al. (1987) sepa- 
rately, to assess the influence this subsample exerts upon 
the full sample. Finally, we fit the full sample and then 
also fit various groupings of RIQs and RLQs. Results 
for this model are given on the left side of Table 7 and 
plotted in Figure 9. 

The general tendencies revealed by Figure 9 are (9a) 
RIQs have X-ray luminosities only modestly greater than 
those of RQQs of comparable optical/UV luminosities, 
whereas RLQs become increasingly X-ray bright relative 
to comparable RQQs as £ uv increases; (9b) when RLQs 
are subdivided by radio loudness, RLQs with R* > 3 
are more X-ray luminous than those with R* < 3; (9c) 
when RLQs are subdivided by radio luminosity, RLQs 
with l T > 33.3 are more X-ray luminous than those with 
4 < 33.3; (9d) RLQs with flat radio spectra are more 
X-ray luminous than those with steep radio spectra, and 
in particular almost all of the most X-ray luminous RLQs 
(with £ x > 28) have flat radio spectra. There does not 
appear to be any grouping of RIQs or RLQs that con- 
tains objects with X-ray luminosities less than those of 
comparable RQQs at any optical/UV luminosities; this 
is broadly consistent with RIQs and RLQs being similar 
to RQQs but with an "extra" source of X-ray emission 
whose strength depends upon radio properties. Over all 
groupings and models, there is a general tendency for 
the RLQs at the highest optical/UV luminosities to lie 
above their best-fit models (to a degree exceeding any 
possible slight systematic flattening of the slope due to 
the fitting method; see Appendix D), and this structure 
in the residuals suggests that a linear fit (to logarithmic 
quantities) of X-ray luminosity as a sole function of opti- 
cal/UV luminosity is not an adequate model even when 
applied within subgroups of RLQs, at least for particu- 
larly radio-loud, luminous, or flat-spectrum RLQs. 

Joint 68% and 90% confidence ellipses for the various 
fits to this model are plotted in Figure 10. In all panels 
the RQQ result is plotted as a black ellipse for compar- 
ison. It can be seen in (10a) that the confidence region 
for RIQs is near to that of RQQs; the modest radio loud- 
ness and radio luminosity of RIQs generally do not ap- 
pear to enhance substantially their X-ray emission. In 
contrast, the confidence region for RLQs is well sepa- 
rated from that of RQQs, with both a greater model 
intercept and slope. It can also be seen that our sam- 
ple substantially increases the precision with which the 
model parameters can be assessed over that provided by 
previous studies, such as that of Worrall et al. (1987) for 
their Einstein sample of RLQs (recall we use 93 of their 
114 objects). In (10b), (10c), and (lOd) the confidence 
regions for RLQs with R* > 3, £ r > 33.3, and a r > —0.5 
are offset from their weaker RLQ counterparts, which are 
themselves still fully distinct from RQQs. The increased 
X-ray brightness for RLQs with R* > 3 or £ r > 33.3 is 
primarily due to a larger intercept in the modeled re- 
lation. It is possible that flatter radio spectra objects 
may have a stronger dependence of £ x on £ nv , although 
the 90% confidence regions overlap in projection onto the 
slope variable. The trends displayed in (10b) and (lOd) 
are qualitatively similar to those shown by Worrall et 
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Figure 11. X-ray luminosity as a joint function of optical/UV and 
radio luminosity for various sub-samples of RLQs. The model con- 
sidered is £ x = ao + b uv X £ uv + c r X £ T (where I x , £ uv , and £ r have 
been normalized prior to fitting as described in §5). Filled sym- 
bols are X-ray detections and larger symbols are more radio-loud 
throughout. The x-axis is based on the best-fit model for RLQs, 
for which c r ~ 0.58xo uv . The black solid line in panels (b), (c), 
and (d) is the best-fit model for RLQs. 

al. (1987) in their Figure 1. 

5.2. ^ x (^uvi ^r) 

5.2.1. Inclusion of radio luminosity as a fit parameter 

The second parameterization we consider is X-ray lu- 
minosity as a joint lunction of optical/UV luminosity and 
radio luminosity: £ x = a + b uv x £ uv + c r x £ r . The result- 
ing coefficients obtained from fitting various groupings of 
RIQs and RLQs are listed on the right side of Table 7. 
The optical/UV luminosity coefficient is now b uv — 0.506 
for RLQs and 6 UV = 0.432/0.439/0.471 for highly radio- 
loud/radio-luminous/flat-spectrum RLQs, compared to 
the 6 UV = 0.649 for RQQs; this suggests that the appar- 
ently stronger dependence of RLQ X-ray luminosity upon 
£ uv indicated in the previous model of £ x (£ uv ) was actu- 
ally reflecting the influence of radio luminosity, which is 
now explicitly considered. The best-fit model for RIQs 
indicates that radio properties do not strongly infuence 
the X-ray luminosity of RIQs (c r is formally consistent 
with zero). Figure 11 plots £ x versus £ uv + 0.58 x £ r . This 
choice of variables is motivated by the coefficients of the 
4(4v,4) model for RLQs, for which c r is ~ 0.58x6 uv 
(Table 7). Collapsing the £ uv — £ r plane to a single joint 
variable simplifies presentation of the modeling results 
and enables ready comparison of the properties of sub- 
groups of RLQs to those of RLQs as a whole. It can 
be seen from (11a) that the Einstein sample of RLQs 
dominates the highest luminosity region of the full sam- 
ple, which also contains the high-redshift sample objects. 
Conversely, the lowest luminosity region of the full sam- 
ple is strongly influenced by the deep-field sample ob- 
jects, although there are many primary sample photo- 
metric quasars in this region as well. In (lib), (11c), and 
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Figure 12. Joint 90% confidence ellipses for fitting X-ray lumi- 
nosity as a function of radio and optical/UV luminosity. The model 
considered is £ x = ao + b uv X £ uv + c r X £ r (where £ K , £ uv , and £ T have 
been normalized prior to fitting as described in §5). The legend in 
(b) also applies to (a) and the legend in (d) also applies to (c). See 
§5.2 for discussion. 

(lid), data for the same subgroups of RLQs as in Fig- 
ure 9 are shown along with the best-fit model 
for RLQs for comparison. The subgroups of RLQs do 
not deviate strongly from the trend for RLQs in general 
in these coordinates, although the particularly radio and 
optical/UV luminous RLQs from the Einstein sample 
(along with a few high-redshift objects) are still exces- 
sively X-ray bright. 

Joint 90% confidence ellipses (calculated after collaps- 
ing the third dimension for ease of viewing) for the vari- 
ous fits to this model are plotted in Figure 12. The full, 
primary, and Einstein samples are plotted in (12a) and 
(12b), as are the subgroups of RLQs with R* < 3 and 
R* > 3. The parameters for the primary sample of SDSS 
RIQs and RLQs are consistent with those of the full sam- 
ple (not unexpected, since the primary sample makes up 
a majority of the full sample). There is no evidence of 
a statistically significant difference in the X-ray lumi- 
nosity dependence of these subgroups upon optical/UV 
luminositvFI There is suggestive support for Einstein 
RLQs and (relatedly) for RLQs with R* > 3 possess- 
ing a stronger dependence of X-ray luminosity upon ra- 
dio luminosity, but the projected confidence ellipses all 
overlap in (12b). The joint 90% confidence ellipses for 
the 4(4vi4) model applied to subgroups of RLQs di- 
vided by radio luminosity and by radio spectral index 
are plotted in (12c) and (12d), along with the result for 
RLQs in general provided for comparison. Flat-spectrum 
RLQs are X-ray brighter than steep-spectrum RLQs due 
to a larger intercept, but have a consistent dependence 
of X-ray luminosity upon both optical/UV and radio lu- 
minosity. In possible contrast, RLQs with £ r < 33.3 and 
those with £ r > 33.3 share similar best-fit intercepts, but 

27 It appears possible that the R* < 3 objects may have a larger 
b uv coefficient (and a smaller c r coefficient); this could be probed 
with a larger sample. 
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the more radio-luminous RLQs may have a greater/lesser 
dependence ol 4 upon 4/4vj although the conhdence 
ellipses overlap in projection onto c T /b uv . 

5.2.2. A "radio-adjusted" 4(tv) relation for RIQs and 

RLQs 

We consider briefly whether RIQs and RLQs can be 
treated as similar to RQQs but with an additional jet- 
linked contribution to the X-ray luminosity (which could 
imply a consistent disk/coronal structure). If so, then 
after accounting for the influence of radio emission, the 
X-ray luminosity in RIQs and RLQs should be correlated 
with optical/UV luminosity through a relation similar 
to that for RQQs. RIQs display no significant depen- 
dence of X-ray luminosity upon radio luminosity and may 
share the same dependence upon optical/UV luminosity 
as holds for RQQs, although this is not strongly con- 
strained in our data. 

One method of investigating a "radio-adjusted" 4(4v) 
relation could be to set the radio luminosity in the best-fit 
RLQ 4(4v> 4) model to a value representative of RQQs. 
This requires an accurate parameterization of 4(4v) for 
RQQs, a difficult function to evaluate given the inherent 
radio weakness of RQQs. A simple model could have ra- 
dio luminosity proportional to optical/UV luminosity as 
£ r = a + (3 x £ uv ; in this case 4(4v, 4) would transform 
to 4(4v) as (adding in the implicit luminosity normal- 
izations) 4 - 27 = a + c r x (a + 30.5/3 - 33.3) + (b nv + 
Cr/3) x (£ uv - 30.5). White et al. (2007) present a cor- 
relation (their Equation 2) that extends to low radio- 
loudness values (R* < 1) and is equivalent!^! to a — 4.57 
and /3 = 0.85. Adopting these values of a and (3, the 
RLQ fit of 4 = -0.100 + 0.5064v + 0.2924 becomes 
4 = —0.919 + 0.754£ uv ; this may be compared with the 
best-fit RQQ relation of 4 = -0.545 + 0.6494 v (all fit- 
ted coefficients from Table 7). The slope for RLQs in 
the radio-adjusted relation is closer to but slightly larger 
than that for RQQs (the difference is 0.105±0.049). Since 
we have not yet taken the likely beaming of some frac- 
tion of the X-ray emission in RLQs into account, this 
result does not mandate that the disk/corona in RLQs 
is more X-ray efficient at high optical/UV luminosities 
than in RQQs. (We demonstrate in §6 that the RLQ 
4(4v, £r) fit can be reproduced assuming a disk/coronal 
scaling as in RQQs plus a jet component.) The differ- 
ence in intercepts as compared to the RQQ relation may 
be reflective of greater beaming of radio emission in the 
RLQs (although some RQQs may have a boosted com- 
ponent of radio emission; e.g., Miller et al. 1993; Fal- 
cke et al. 1996), in which case the RQQ relation ought 
to have the radio luminosity similarly enhanced prior to 
substitution for a first-order comparison. For illustrative 
purposes, accounting for an additional beaming enhance- 
ment of a factor of 19 (e.g., corresponding to an inclina- 
tion of ~ 12° with 7 = 10.5; see §6) would change the 
transformed RLQ intercept to match the RQQ result. 

5.3. 4(4vi £r: 

5.3.1. Inclusion of redshift as a fit parameter 

We now investigate whether there is a dependency 
upon redshift in addition to the dependencies upon op- 

28 Their radio-loudness has been slightly adjusted as a function 
of optical luminosity; see their Equation 4. 
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Figure 13. Residuals for best-fit results for X-ray luminosity pa- 
rameterized by optical/UV and radio luminosity, plotted against 
redshift and fractional lookback time. The legend in (b) also ap- 
plies to (a) and the legend in (d) also applies to (c). A single object 
may be a member of two categories in a given panel, in which case 
the residual values from the separate fits to each category are each 
plotted at the redshift of that object. There is no apparent redshift 
dependence. See §5.3 for discussion. 

tical/UV and radio luminosity. Although the £ x (£ uv ,£ r ) 
residuals (see Figure 13) do not show any obvious red- 
shift dependence™ it is in principle possible that some 
of the apparent luminosity dependence might actually 
be driven by redshift evolution. We test this possibil- 
ity by including the redshift dependence as a parameter 
when modeling, by fitting 4(4v,4 z )? where again the 
luminosities are normalized prior to fitting. The redshift 
dependence is put in terms of log(l + z) or lookback frac- 
tion [t z — 1— age(z)/age(z = 0)]. The best- fit model for 
the full sample is 4 = (-0.098±0.051) + (0.482±0.036) x 
£ uv + (0.273 ± 0.019) x l T + (0.000 ± 0.128) x log(l + z); 
when RLQs only are considered, the best-fit model is 
4 = (-0.109 ± 0.052) + (0.503 ± 0.042) x £ uv + (0.293 ± 
0.027) x £ T + (0.025 ± 0.131) x log(l + z). When the 
redshift dependence is expressed instead in terms of 
the fractional lookback time r z , the best-fit model for 
the full sample is 4 = (-0.064 ± 0.072) + (0.488 ± 
0.034) x4v + (0.273 ±0.019) x £ r + (-0.052 ±0.109) x r z ; 
when RLQs only are considered, the best-fit model is 
4 = (-0.065 ± 0.074) + (0.511 ± 0.042) x £ uv + (0.292 ± 
0.027) x 4+ (-0.053 ± 0.111) x t z . In all these cases the 
difference between the coefficient for the redshift term 
and zero is not statistically significant, and the joint 68% 
and 90% confidence ellipses with the optical/UV and ra- 
dio luminosity coefficients include zero (Figure 14). The 
redshift coefficients for the log(l±z)/r z models when fit- 
ting the primary sample, spectroscopic sample, or Ein- 
stein sample alone are -0.500 ± 0.167/-0.364 ± 0.140, 
-0.351 ± 0.245/-0.337 ± 0.202, or 0.775 ± 0.443/0.282 ± 

29 The possible tendency for RLQs at high redshift to have 
positive residuals is not necessarily a redshift effect; the best-fit 
ix(^uv,ii) model appears to underpredict X-ray emission for par- 
ticularly radio-loud or radio luminous RLQs such as these (see also 
§2.2.2). 
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Figure 14. Consideration of redshift dependence. The first model 
investigated is £ x = ao + b uv x£ uv + c r x£ r + d z xlog(l + z) (where 
^x, ^uv, and £ T have been normalized prior to fitting as described 
in §5); the second model investigated replaces log(f + z) with the 
fractional lookback time r z . Joint 90% (solid) and 68% (dotted) 
confidence ellipses are plotted for the full sample (red) and for the 
sample of RLQs (blue). The difference between the coefficients for 
redshift dependence and zero is not statistically significant. 



0.379, respectively. The indicated redshift dependence 
is marginal (3.0/2.6, 1.4/1.7, or 1.7/0.7<r), and the di- 
rection of the potential trend is inconsistent between the 
primary and Einstein samples; recall also that Worrall 
et al. (1987), using independent methods, found no sig- 
nificant redshift dependence among the Einstein objects. 
The superior coverage of the t — z plane provided by the 
full and RLQ samples leads us to favor those results. 

We find no evidence for a dependence upon redshift for 
the X-ray emission properties of the full sample or RLQs 
alone. The degree to which otherwise comparable RIQs 
or RLQs at different redshifts could differ in X-ray lumi- 
nosity may be constrained via the value and la errors on 
the best-fit coefficient to the redshift term. For the full 
sample, the coefficient of (0.000 ±0.128) x log(l + z) sug- 
gests a maximum evolution (la, i.e., using 0.000+0.128 
for the coefficient) in X-ray luminosity between redshift 
z = and z = 5 of £ x ,z=o — ^ x ,z=5 = 0.100, or a ratio in 
linear units of 1.26. This suggests that the X-ray lumi- 
nosity of otherwise comparable RIQs or RLQs has not 
changed by more than 30% over z = — 5. For RLQs 
only, a similar analysis suggests a maximum change of 
<35%. The r z best-fit coefficients similarly suggest a 
maximum redshift-driven change in X-ray luminosity for 
the full sample or for RLQs of <30% over z = — 5. 
Although the fraction of quasars that are RLQs is de- 
pendent upon redshift (e.g., Jiang et al. 2007), the X-ray 
properties of individual RIQs and RLQs do not appear 
to differ strongly when comparing objects at low versus 
high redshift. Apparently the cosmic evolution in the 
efficiency of generating RLQs does not substantially im- 
pact RLQ structure post-formation. 



5.3.2. Dependence of relative X-ray brightness on redshift 

An alternative manner of visualizing the potential red- 
shift dependence of the X-ray luminosity of RIQs and 
RLQs and of searching for redshift evolution as a func- 
tion of radio-loudness is to construct a multiwavelength 
color-color plot and group objects within separate red- 
shift bins. Figure 15a shows ct ox plotted versus R* for the 
full sample divided into low, medium, and high redshift 
(z < 1, 1 < z < 2, and z > 2, respectively). Objects be- 
come increasingly X-ray bright (with less negative values 
of a ox ) as they become increasing radio-loud, as previ- 
ously shown in Figure 7, but there is no obvious difference 
in this trend with redshift. Fitting a ox (R*) does suggest 
that the slope steepens slightly with redshift, but this is 
offset by a decreased intercept (Figure 15c; see also Fig- 
ure 5 of Lopez et al. 2006). It is preferable to take the 
dependence of X-ray luminosity upon optical/UV lumi- 
nosity into account when constructing the optical/UV- 
to-X-ray color. We calculate Aa ox = a ox — a ox (^ U v) 
by taking a ox (^uv) from the Just et al. (2007) relation 
for RQQs (their Equation 3). Figure 15b shows Aa ox 
plotted versus R* for the full sample grouped into low, 
medium, and high redshift bins. There is no obvious 
stratification by redshift, but here fitting Aa ox (i?*) sug- 
gests that the slight increase in slope with increasing 
redshift may result in a slight increase in relative X-ray 
brightness, at least for RLQs with R* > 2 (Figure 15d). 
It is possible that a disproportionate number of objects 
in the high-redshift bin with low line-of-sight inclinations 
could produce this result. The best-fit model parameters 
are given in Table 8. Two-dimensional KS tests indicate 
that the distribution of a ox — R* is inconsistent between 
different redshift bins (low-mid, mid-high, or low- high), 
but that of Aa ox — R* is not (p>0.12 for all cases). 

It has been suggested that X-ray jets for which in- 
verse Compton scattering of cosmic microwave back- 
ground photons (IC/CMB; e.g., Tavecchio et al. 2000) 
provides the primary emission component should become 
the dominant source of RLQ X-ray emission at high red- 
shift (e.g., Rees & Setti 1968; Schwarz 2002). We esti- 
mate the enhancement in total X-ray luminosity from 
this process for sources at redshifts of 0.69/1.44/2.48 
(the medians within the low/medium/high groupings) 
and for a jet-linked fraction j of the nuclear X-ray emis- 
sion (prior to accounting for the redshift dependence 
of IC/CMB X-ray emission) of 1%, 10%, and 99%. 
This is expressed in terms of an increase in Aa ox as 
0.384 x log [j(l + z) 4 + (1 — j)], where the coefficient con- 
verts the increase in £ x to a flattening of the optical/UV- 
to-X-ray spectral slope. The results are plotted (in in- 
creasing order of jet dominance) at radio-loudness val- 
ueJ^I of 1.5, 3, and 4.5 in Figure 14d. It can be seen 
that the (1 + z) 4 dependence of X-ray jet-linked IC/CMB 
emission would lead for these parameters to a stronger 
splitting with redshift than is observed. We find that 
the X-ray luminosities within our sample of RIQs and 
RLQs are unlikely to include significant contributions 
from X-ray IC/CMB jet emission. In principle it is pos- 

30 These R* values increase with increasing jet dominance, as 
is qualitatively expected, but the precise numerical association be- 
tween R* and X-ray jet dominance is model-dependent (as explored 
in §6.3) and so the chosen R* values should be regarded as illus- 
trative. 
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Figure 15. Consideration of relative X-ray brightness as a function of radio-loudness. The full sample of RIQs and RLQs is shown in (a) 
for a ox (R*) and in (b) for Aa ox (if), with blue, green, and red colors indicating low, medium, or high redshift (a < 1, 1 < z < 2, and 
z > 2, respectively). Best-fit trend lines are shown in (c) and (d) with the la errors on the slopes indicated by dashed lines. The expected 
enhancement in X-ray brightness from the (1 + z) 4 dependence of IC/CMB jet-linked X-ray emission is also plotted in (d) for comparison; 
the set of points at if =1.5, 3, and 4.5 correspond to models in which the jet-linked fraction of the X-ray continuum (prior to considering 
redshift) is 1%, 10%, and 99%, respectively. 



sible that extremely bright X-ray jet features could be 
routinely found outside of the extraction regions we use 
to measure the X-ray luminosity, but to our knowledge 
such sources are extremely rare. This is consistent with 
the general lack of observed high-redshift RLQs in which 
the X-ray jet outshines the core (e.g., Bassett et al. 2004; 
Lopez et al. 2006). 

5.4. Additional considerations 

We briefly examine four additional topics relevant to 
the analysis of luminosity correlations: the impact of the 
selection method, the influence of targeted sources, the 
normality of the variables, and the potential for spurious 
or inaccurate results due to luminosity dispersion effects. 



5.4.1. Optical versus radio selection 

As discussed in §2.1.1, the median properties of a 
quasi-radio-selected sample of RIQs and RLQs (con- 
structed from SDSS spectroscopic quasars by requiring 
the "FIRST" target flag to be set) are similar to those of 
"QSO/HIZ" targeted SDSS spectroscopic quasars (with 
substantial overlap between these samples). The depen- 
dence of X-ray luminosity upon optical/UV and radio 
luminosities is also similar for the "FIRST" and for the 
"QSO/HIZ" samples, as may be seen from the best-fit re- 
lations given in Table 7. Joint 68% and 90% confidence 
contours from fitting the "FIRST" and the "QSO/HIZ" 
samples to £ K — ao + b uv x £ uv + c r x £ x are provided in 
Figure 16, where it can be seen that the parameter values 
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Figure 16. Comparison of a quasi-radio-selected (red ellipses) 
sample of RIQs and RLQs to an optically-selected (blue ellipses) 
sample (both drawn from SDSS/FIRST data; there is >90% over- 
lap), illustrated with joint 90% (solid) and 68% (dotted) confidence 
ellipses for the model i x = ao + b uv X £ uv + c r X £ r . The axes are 
scaled to match Figure 12. The RIQs and RLQs targeted by SDSS 
as FIRST sources have properties consistent with those targeted 
due to optical colors. 

are consistent. By excluding the small fraction of objects 
with A(g — i) > 1 we omit strongly dust-reddened RLQs, 
which may be associated with young quasars (e.g., Urru- 
tia et al. 2008, 2009). Within the relative color range we 
accept, there does not appear to be a difference in the 
X-ray properties of SDSS RIQs and RLQs selected by 
color compared to those selected due to radio emission. 
Presumably the X-ray emission mechanisms are likewise 
similar, and so the results from §5 and the modeling in 
§6 are generally applicable. 

5.4.2. Influence of targeted objects 

Only 14.1% of RIQs and RLQs within the primary 
sample were targeted for X-ray observations (§3.4), with 
the remainder possessing serendipitous off-axis cover- 
age. Inclusion of the targeted sources does not bias the 
results of the luminosity correlation analysis, although 
there are some minor differences in the properties of tar- 
geted and off-axis objects. Compared to those objects 
observed serendipitously, the targeted RIQs and RLQs 
are at somewhat lower redshifts (median z of 1.05 versus 
1.45) but are also brighter (median m, of 17.94 versus 
19.54) and so have modestly higher optical/UV lumi- 
nosities (median £ uv of 30.74 versus 30.33). The targeted 
RIQs and RLQs are also somewhat more radio-loud (me- 
dian R* of 3.09 versus 2.44) and -40% X-ray brighter 
relative to optical/UV luminosity (median £ x — I x ,kqq of 
0.53 versus 0.38, with £ x — ^ x ,rqq as defined in §4). Pre- 
sumably, targets were preferentially selected from sources 
already known to be X-ray bright. 

The results of parameterizing X-ray luminosity for the 
off-axis objects are similar to those obtained for the full 
and primary samples (Table 7) . The differences between 
the best-fit parameters for the slopes for the off-axis sam- 
ple compared to the full or primary samples are essen- 
tially zero(3 Fitting the targeted objects separately indi- 
cates they possess a weaker dependence of X-ray luminos- 
ity upon optical/UV luminosity and perhaps a stronger 




31 Differences between the off-axis and full/primary samples are 
(-0.045 ± 0.044)/(0.023 ± 0.048) for the 4(<„) model 



Ab uv 

and A6 UV = (0.046 ± 0.051)/(0.048 ± 0.052), Ac 
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Figure 17. Evaluation of normality for residuals from fitting X- 
ray luminosity as a joint function of optical/UV and radio lumi- 
nosity (RIQs and RLQs, top panel; RLQs, middle panel; fits from 
Table 7) and as a sole function of optical/UV luminosity (RQQs, 
bottom panel; data from Steffen et al. 2006). The open histograms 
include limits and are fit with Gaussians as indicated by the dotted 
curves, while the filled histograms are detections only and are fit 
with Gaussians as indicated by the solid curves. 

dependence upon radio luminosity than do the full or pri- 
mary samples, although their parameters are consistent 
with those of the spectroscopic sample. 

5.4.3. Log-normality of residuals 

A presumption that considered variables are normally 
distributed is inherent in the chosen method of analy- 
sis (and indeed in most parametric modeling), and can 
be checked by examining the residuals from the best- 
fit models for X-ray luminosity in our sample of RIQs 
and RLQs. Histograms of the residual X-ray luminos- 
ity for RIQs and RLQs, RLQs alone, and RQQs are 
provided in Figure 17. While it is not clear a priori 
that the RIQs and RLQs considered here ought to show 
normally-distributed residual X-ray luminosity (in log- 
arithmic units), since there are distinct X-ray emission 
mechanisms for RIQs and RLQs (disk/corona and jet- 
linked) as opposed to a single dominant X-ray emission 
mechanism (disk/corona) in RQQs, it may be seen in 
Figure 17 that the distributions closely match the over- 
plotted best-fit Gaussians. Unfortunately the presence 
of limits in our sample makes it difficult to utilize stan- 
dard tests for normality. We examine the distribution 
of detected objects and also the distribution of all ob- 
jects, treating limits as detections for this calculation 
only. The distribution of residual X-ray luminosity is 
acceptably characterized as normal, with KS test proba- 
bilities p > 0.18 in all cases. While this is best regarded 
as suggestive rather than conclusive, it does indicate that 
our methodology is unlikely to produce biased results due 
to substantial underlying non-normality. 

There appears to be a tail of X-ray weak objects within 
the RQQ sample (17c), some of which might be low- 
rcdshift BAL RQQs. The sample of RIQs and RLQs 
(17a) also appears to show an X-ray weak tail, perhaps 
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less prominently for RLQs only (17b). These may again 
be BAL objects, which decrease in percentage as radio- 
loudness increases (e.g., Shankar et al. 2008). However, 
in contrast to the situation for RQQs, there may also be 
an X-ray bright tail for RIQs and RLQs (for example, 
the percentage of objects above 1.645tr from the mean 
is not 5% but rather 8%/8% for (RIQs+RLQs)/RLQs, 
calculated using the fitted values for a and /i; it is 6% 
for RQQs). These X-ray bright objects are generally 
members of the set of particularly radio-loud or radio- 
luminous RLQs, and hence may be jet-dominated at 
X-ray frequencies. 

5.4.4. Robustness of luminosity correlations 

The degree of inherent scatter in the luminosities can 
affect the results of correlation studies. For example, 
Yuan et al. (1998) demonstrated that when the opti- 
cal/UV luminosity scatter a uv is significantly larger than 
the X-ray luminosity scatter cr x , a fit to a ox (£ uv ) over 
a limited range of £ nv can indicate an anti-correlation 
where none necessarily exists. Relatedly, <t uv 3> cr x can 
give a fitted slope for the £ x (£ uv ) relation less than unity 
even if these luminosities are actually proportional (when 
the considered luminosity ranges are small). In gen- 
eral, examined luminosity correlations are significantly 
less likely to yield spurious or inaccurate results with a 
sample spanning a large range in luminosities relative to 
the observed dispersion (e.g., see discussion in §5 of Just 
et al. 2007), as does our sample. 

It is possible to assess the luminosity scatter through 
consideration of the standard deviation of the residuals 
about the £ x (4v,4) relation. For the primary sample, 
these values are 0.34/0.71/1.30 for £ x /£ uv /£r, with corre- 
sponding 5th-95th percentile luminosity ranges of 25.81- 
27.62/29.40-31.26/31.47-34.46 (treating X-ray limits as 
detections for this calculation only). This fit is performed 
treating £ x as the dependent variable, so the residu- 
als about this fit are larger than the inherent scatter. 
We use a Monte Carlo approach to determine the corre- 
sponding cr x /a uv /cr T values. For simplicity, we simulate 
^uv,mod as uniformly distributed (n = 654) between 29.50 
and 31.25, with 4, mod = 4v,mod + 2.5 (with random 
log-normal scatter of <7r» ) and ^ Ximo d calculated from 
£x{£uv, £r)- Then random log-normal scatter cr x /a uv /a v 
is added to the model luminosities and the standard devi- 
ations of the residuals and the luminosity ranges are cal- 
culated. A good match to the observed data is provided 
with a x /a uv /a r = 0.28/0.32/0.420 and cr R . = 0.55 gives 
a radio-loudness standard deviation of ~0.75 in agree- 
ment with the primary sample. 

We assess the impact of such inherent scatter by fit- 
ting the simulated data (including X-ray censoring as 
described in §6.1.3). This scatter acts to depress the de- 
pendence of X-ray luminosity upon optical/UV and radio 
luminosities, so that fitted values of b nv and c r are slightly 
lower than those used in the input £ x (£uv,£r) relation. 
The observed best-fit values of & uv — 0.48 and c r ~ 0.26 
for the primary sample would result in fitted coefficients 
of 0.39 and 0.23, respectively, and may themselves be 
obtained for input coefficients of ~0.6 and ~0.3, respec- 

32 These are larger than the uncertainty typically associated 
with variability, plausibly suggesting a spread in physical properties 
among RIQs and RLQs. 



tively. Further, data simulated with an input model of 
£x{£\iv:£r) that is then fit as £ x (^uv 

) produce fitted val- 
ues for & uv greater than the input coefficient (e.g., fitted 
6 UV = 0.67±0.03 for input 6 uv /c r =0.6/0.3), as suggested 
in §5.2.1. 

A contribution to the X-ray luminosity linked to the 
radio luminosity appears to be mandated for RIQs and 
RLQs. For the derived inherent scatter, it is not possible 
to reproduce the observed best-fit £ x (£ U v,£r) relation for 
the primary sample with an input model using £ x (£ uv ) 
as for RQQs. With the luminosity ranges of our primary 
sample, the & uv — 0.65 observed for RQQs may result 
from an input coefficient of ~0.85 (1.0 is excluded; see, 
e.g., §3.5 of Strateva et al. 2005). Data simulated with 
an input model of £ x tx 0.85 x £ uv that is then fit as 
£*(£uv,£r) produce fitted coefficients of & uv = 0.56 ±0.03 
and c r = 0.09 ± 0.02. If £ r is initially taken directly 
from £ uv (i.e., ctr* = 0), the fitted coefficients are 6 UV = 
0.47 ±0.03 and c r = 0.22 ±0.03, but the scatter in radio- 
loudness is then significantly lower than observed. (Using 
4 = a + X 4 V with a = 4.57 and /3 = 0.85 as in §5.2.2 
gives similar results.) 

We note for completeness that more extreme val- 
ues of er x / 'a uv I '<7 r could distort luminosity correlations 
(e.g., for cr x /CT uv /cr r = 0.28/0.64/0.42 or 0.28/0.32/0.84 
rather than the derived 0.28/0.32/0.42, input parame- 
ters of 6 U v/c r =0.6/0.3 give fitted coefficients of 0.22/0.33 
or 0.56/0.15, respectively). However, as demonstrated 
above, these considerations do not apply to our sample. 
Additionally, since much of our analysis is concerned with 
comparing the results of (similarly conducted) fits to dif- 
ferent groupings of RIQs and RLQs, any modest sys- 
tematic skewing of correlations does not impact our con- 
clusions regarding the relative differences in how X-ray 
luminosity is dependent on optical/UV and radio lumi- 
nosities within these subgroups. 

6. A PHYSICAL MODEL FOR X-RAY EMISSION IN RIQS 
AND RLQS 

It is of interest to evaluate the physical basis for the 
correlations discussed in §5. It is widely theorized that 
the increasing X-ray brightness of RIQs and RLQs with 
increasing radio loudness or luminosity is driven by a 
source of nuclear X-ray emission that is directly or in- 
directly powered by the radio jet (e.g., see discussion in 
§1), but the precise nature of this linkage is not clearly 
understood. We make use of previous results from the lit- 
erature to simulate a population of RIQs and RLQs with 
radio and optical properties consistent with observations, 
and then test competing models for the X-ray emission 
through comparing the properties of the simulated data 
sets to observations. The structure and parameters of 
this modeling are given in Table 9. 

We adopt a physical model that contains emission con- 
tributions from the core and lobes at radio frequencies 
(making the common assumption that the core is dom- 
inated by the small-scale radio jet, but see also Bell & 
Comeau 2010), from the disk and small-scale jet at op- 
tical/UV frequencies, and from the disk/corona and an 
additional "jet-linked" component at X-ray frequencies. 
Most of the parameters for the radio and optical emis- 
sion components in this model are fixed by prior work; 
likely values for a few free parameters were determined 
through comparison to our sample data (see Table 9 for 
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details). The X-ray disk/corona emission is presumed to 
scale with the optical/UV disk emission as established for 
RQQs, as seems reasonable based on the results of §4 and 
§5. The Doppler beaming factor is 8 = 7~ 1 (1 — /Jcos^) -1 
(e.g., Worrall & Birkinshaw 2006). 

We consider three possibilities for the X-ray jet-linked 
emission: 

• In model A, the X-ray jet-linked emission is pro- 
portional to the intrinsic radio-jet emission (prior 
to applying beaming) and is also itself unbeamed. 

• In model B, the X-ray jet-linked emission shares 
the beaming factor S that applies to the radio-jet 
emission. 

• In model C, the X-ray jet- linked emission has a 
lower bulk Lorentz factor j x (and thus a lower 
beaming factor <5 X at low inclinations) than the 
radio-jet emission. 

The models A, B, and C correspond in a general sense 
to cases in which the jet-linked X-ray emission originates 
in an additional accretion flow structure, or co-spatially 
with the radio-jet emission, or within the small-scale jet 
but predominantly in a less-relativistic region, respec- 
tively. Note that we are not attempting to model the 
X-ray emission in individual sources (such an approach 
is unproductive with these data). We instead conduct a 
statistical study to examine different plausible and repre- 
sentative physical scenarios, and present the consequent 
implications including, for example, the implied fraction 
of RIQs and RLQs for which the X-ray emission is jet- 
dominatcd. 

In this section, we first describe the various compo- 
nents of the basic model, then match the simulated lu- 
minosities to the observed primary sample data, and 
then consider additional observational constraints upon 
the models and determine the dominant source of X-ray 
emission in each case. 

6.1. The model components 

6.1.1. Radio emission: core and lobes 

We follow the general unification scheme described by 
Jackson & Wall (1999), in which FR I and low-excitation 
emission-line FR II radio galaxies are the parent pop- 
ulation of BL Lacs, while high-excitation emission-line 
FR II radio galaxies are the parent population of RLQs 
(cf. Donoso et al. 2009), but update the beaming model 
and luminosity function to reflect more recent consen- 
sus. The radio source population is presumed to be de- 
scribed by the luminosity function presented by Willott 
et al. (2001), which is based on low-frequency (151 MHz 
and 178 MHz) data and is thus relatively unbiased to- 
ward beaming, and is computed at luminosities suffi- 
ciently high that contamination from star-forming galax- 
ies is negligible (Willott et al. 2001 estimate it to be only 
~2% even at their 0.1 Jy limit). The simulated objects 
include radio galaxies and quasars; we select for the lat- 
ter simply by requiring the inclination to be < 60° (we 
also require 9 > 5.8° to remove highly beamed objects, 
which our sample attempts to exclud . A sample of 

33 An indication that ~ 5° is reasonable is provided by applying 
the orientation measure from Wills & Brotherton (1995) to the 



RIQs and RLQs is synthesized with redshift and luminos- 
ity distributions drawn from the luminosity function, and 
with randomly assigned orientations (uniform in sin 6) . 
The intrinsic core prominence (the ratio of core-to-lobe 
radio flux at low frequencies; i.e., unaffected by beaming) 
is taken from the Bayesian modeling of FR II sources car- 
ried out by Mullin & Hardcastle (2009) and is simulated 
including intrinsic scatter in core power based on their 
best-fit model; we also take the typical bulk Lorentz fac- 
tor 7 = 10.5 for core emission from their work. The 
parameters we adopt are those from Table 5 of Mullin & 
Hardcastle (2009) for the model excluding low-excitation 
emission-line objects. 

The observed radio characteristics for a given simu- 
lated source are calculated for the inclination of that 
source. The core emission from the small-scale jet is 
boosted by S 2 ~ a ' (e.g., Worrall & Birkinshaw 2006). The 
1.4 GHz flux densities for the lobes and core of each sim- 
ulated source are determined assuming a t — —0.9 for 
lobe emission and a r = —0.3 for core emission (consis- 
tent with the methodology used in §3.2 to calculate £ T for 
the sample sources from the observed FIRST flux densi- 
ties). A limit of 1 mJy was imposed to match the FIRST 
catalog detection limit, and components with simulated 
flux densities below this limit were dropped from further 
consideration. Faint and diffuse lobes may not always 
register as FIRST catalog sources; for only ^5% of the 
primary sample objects with lobes is the integrated lobe 
flux below 2.7 mJy, so we additionally do not include 
any contribution from lobes with simulated fluxes below 
2.7 mJy in the calculated radio luminosities. 

6.1.2. Disk- dominated optical/UV emission 

Optical emission from RLQs can be generated by 
quasi-thermal emission from the accretion disk and by 
nonthermal (e.g., synchrotron) emission from the small- 
scale jet (e.g., Wills et al. 1995). In our sample, the 
ordinary equivalent widths of the broad emission lines 
place an upper limit upon the degree to which a fea- 
tureless jet-linked component can contribute to the opti- 
cal/UV emission. However, it is not possible in practice 
to determine simply the fraction of optical/UV emission 
that is disk-linked in individual objects from either the 
optical/UV spectra^ or from the radio multifrequency 
data0 The observational indications that disk emission 
remains dominant in broad-line RIQs and RLQs there- 
fore guide our models but do not provide sufficient moti- 
vation to fix simply the fraction of optical/UV emission 
that is disk-linked. 

most radio-loud group of primary sample sources. When we also 
model the optical/UV emission (§6.1.2) we find that a cutoff within 
4° — 7° is required to match the observed distributions of both £ T 
and£ uv . The adopted 8 > 5.8° helps provide a simulated £ x (£ U v, (-i) 
relation (§6.1.3) matching that observed for the primary sample. 

34 Such an approach would have to account for the known trends 
in emission-line strength with luminosity (Baldwin effects; e.g., 
Baldwin 1977; Osmer & Shields 1999) that occur even absent op- 
tical jet emission, but there is still significant scatter in broad line 
strength between otherwise similar objects and so the accuracy of 
this method is fundamentally limited for individual objects. 

35 In principle, the observed radio core flux could be extrapo- 
lated to the optical/UV band using the measured radio spectral 
index. However, even ignoring the uncertainty in the spectral in- 
dex, any intervening spectral break (which is common, even in a 
simple synchrotron context; e.g., Worrall & Birkinshaw 2006) es- 
sentially destroys the accuracy of the extrapolation. 
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We estimate the jet-linked optical/UV component from 
the jet radio emission presuming a ro = —0.8, reflective of 
the tendency for the synchrotron spectrum of jet emission 
to steepen at higher frequencies (e.g., Worrall & Birkin- 
shaw 2006). The disk- linked component is dependent 
upon physical processes not directly incorporated into 
our modeling (most notably the accretion rate); fortu- 
nately, there are radio characteristics which are thought 
to be likewise sensitive to such processes, and so it is 
possible to estimate the optical disk luminosity from the 
already- mo deled radio properties. In particular, Willott 
et al. (1999) find a correlation between low-frequency 
radio power and narrow-line emission in radio galaxies 
which they argue indicates an underlying dependence 
upon accretion rate (to feed the jet and furnish sufficient 
ionizing photons incident upon the narrow line region; 
Willott et al. 1999 also comment that direct illumina- 
tion from the jet or jet-cloud collisions are generally of 
secondary importance but are relevant in some individ- 
ual objects). Motivated by such relationships, we set the 
disk to have a monochromatic luminosity at rest-frame 
2500 A that is a fixed fraction of the intrinsic (unbeamed) 
jet monochromatic luminosity at rest-frame 5 GHz. A log 
offset of —1.8, or a fraction of 1.6%, provides a reasonable 
match to observed £ uv ; values of —1.7 to —2.0 (2-1%) 
are viable. The correlation between low-frequency radio 
power and narrow-line emission observed by Willott et 
al. (1999) has a la scatter of 0.5 dex; we introduce into 
the model a similar spread in disk optical/UV emission 
at a given intrinsic jet luminosity by adding random nor- 
mal scatter with a = 0.5 dex to the disk optical/UV 
emission. We emphasize that this methodology is not 
reflective of a direct physical link between jet and disk 
emission but rather captures the influence of additional 
processes, particularly the accretion rate, on both jet and 
disk emission. 

For most of the simulated sources the disk emission 
dominates over that from the jet at optical/UV wave- 
lengths (by factors of a few to several hundred), and it 
is only for particularly low-inclination objects that the 
jet emission contributes significantly to the optical/UV 
emission. This is consistent with the observed opti- 
cal/UV spectra for our sample RIQs and RLQs. 

6.1.3. Dual X-ray emission components 

Chandra and XMM-Newton observations of FR II ra- 
dio galaxies show that their X-ray spectra often contain 
two components, which can be interpreted as emission 
from both a disk/corona and a jet-linked component 
(e.g., Evans et al. 2006; Hardcastle et al. 2009). The 
X-ray spectra of powerful RLQs, on the other hand, can 
typically be fit with a single power-law model; this sug- 
gests they are dominated by jet emission, generally in- 
ferred to be inverse Compton radiation (e.g., Belsole et 
al. 2006). Guided by prior work and the results discussed 
in §4 and §5, we include in our model X-ray emission from 
a disk/corona and also from the small-scale jet. Our sam- 
ple includes objects of modest radio loudness and lumi- 
nosity, including RIQs; since these observed properties 
are dependent on intrinsic power but also viewing angle, 
we expect our sample to include objects at inclinations 
intermediate between radio galaxies and luminous RLQs, 
allowing us to investigate in our simulations the increas- 
ing jet contribution to the X-ray continuum as inclination 
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Figure 18. Simulated populations of RIQs and RLQs com- 
pared with the observed primary sample. The models differ 
principally in the degree of beaming of the jet-linked X-ray 
emission: A is unbeamed, B is beamed with 7 X = 10.5 as 
for the radio jet, and C is beamed to a lesser degree with 
7 X = 2.0. The left-hand side shows histograms of radio (red), 
optical/UV (green), and X-ray (blue) luminosity, for obser- 
vations (dashed) and simulation (solid). The right-hand side 
shows the a ox — relation for the primary sample (black 
crosses) and for the simulation, with objects color-coded by 
inclination (blue is 5.8° < 9 < 10°, green is 10° < 9 < 20°, 
and red is 20° < 9 < 60°). The solid line is the RQQ relation 
from Just et al. (2007). 

decreases (with scatter reflective of intrinsic variance in 
unbeamed radio core luminosities). 

It is not possible to leave both the disk/corona and jet- 
linked X-ray emission as simultaneously free parameters 
in the absence of additional observational input (such 
as high-quality X-ray spectra) that most of our sources 
lack. We assume the X-ray emission associated with 
the disk/corona is related to the optical disk emission 
in the same manner as for RQQs (note that this does 
not necessarily require identical accretion structure) and 
can therefore be calculated using Equation 7 from Just 
et al. (2007) and the simulated optical/UV disk emission. 
Random log normal scatter with a = 0.3 dex is added to 
the disk/coronal X-ray luminosities to mimic the scatter 
in the RQQ relation. 

The jet-linkec£3 X -ray emission is unlikely to be syn- 
chrotron emission from the same population of electrons 
as generates the radio-jet emission (in contrast to the 
apparent situation for lower power FR I jets; e.g., Chi- 
aberge et al. 2000), as this process does not appear to 
generate sufficient X-ray emission to match observations 
(e.g., Landt et al. 2008). Inverse Compton processes 
must be considered; sources of seed photons include emis- 
sion from the jet (self-Compton, or SSC), radiation from 

36 The use of the term "jet-linked" with respect to model A 
should be understood to refer to the correlation between X-ray 
and intrinsic (unbeamed) radio-jet emission; it is possible but not 
required that this X-ray component is produced in a jet. 
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the central engine (external Compton, or EC), and the 
cosmic microwave background (IC/CMB). The IC/CMB 
model has a strong dependence upon redshift, which con- 
flicts with the lack of redshift dependence observed in our 
sample (see §5.3), and in any case is unlikely to produce 
as many seed photons as can the nucleus on the relevant 
parsec-scale or smaller distances (e.g., Schwartz 2002). 
Comparison of the SSC and EC processes suggests the 
latter dominates this close to the central engine (e.g., 
Sokolov & Marscher 2005). 

The precise physical parameters governing X-ray jet 
emission (such as the size of the emission region or the 
magnetic-field strength) likely vary significantly from ob- 
ject to object. Since we are interested in general trends 
rather than the specific details of a given source, we 
assume a "standard" intrinsic ratio between radio and 
X-ray jet-linked emission, which we set to obtain con- 
sistency with the observed primary sample X-ray lumi- 
nosities. The intrinsic jet-linked X-ray emission is then 
modified in models B and C for the observer based on the 
bulk velocity and inclination of the jet, with presumed 
dominant EC emission boosted by S 3 ~ 2a * (an additional 
factor of 1 — a beyond the radio synchrotron emission; 
e.g., Dermer 1995), with a x — —0.3 for model B taken to 
match the radio value and a x = —0.5 for model C taken 
to match the typical energy index of RLQs (e.g., Page et 
al. 2005). X-ray censoring is added to the simulated sam- 
ple in quasi-random fashion, with fainter objects more 
likely to be labeled as upper limits. For simplicity, this is 
accomplished through matching the X-ray detection frac- 
tion within the top/bottom half of the primary sample 
ranked by £ x . We randomly select 88%/77% of objects 
with simulated £ x > 26.8/ < 26.8 as detections, and then 
the simulated detection rate matches that of the primary 
sample overall and within each X-ray luminosity bin. 

6.2. Comparison to observed luminosities 

In principle, this modeling process will produce a simu- 
lated population of PJQs and RLQs that can be directly 
compared to the observed primary sample data. How- 
ever, the nature of our sample introduces a few com- 
plicating effects that need to be considered. One im- 
portant parameter is the limiting magnitude for inclu- 
sion within the primary sample. Since we prefer to re- 
tain as many observed objects as possible, we do not 
impose a magnitude limit on the observed objects. How- 
ever, only ~5% of primary sample RIQs and RLQs have 
mi > 19.95/20.13 and m K > 20.60/21.08 for z < 1/z > 1 
spectroscopic and photometric quasars, respectively, and 
so these values are used as effective magnitude cutoffs 
for the simulations. The numbers of simulated "spectro- 
scopic" and "photometric" objects (distinguished in the 
model only by the applied magnitude cutoffs) are set to 
match those of the primary sample. 

The normalization (i.e., number of objects) of the sim- 
ulated population should be able to be set solely based 
on the sky coverage fraction of the observed sample, but 
this is difficult to evaluate for the various X-ray missions 
and their complex instrumental fields of view. Matched 
total numbers of simulated and sample objects are not 
critical for conducting a comparison (it is the distribu- 
tion of luminosities that is of interest); nonetheless, it is 
convenient for examining the results of the simulations 
if the simulated population is of the same size as the 



observed sample. Since we are interested in how X-ray 
luminosity is related to radio and optical/UV luminosi- 
ties and unconcerned with the sky density of (selected) 
RIQs and RLQs, we simulate a large population and then 
draw from it to obtain an equal number of simulated 
and observed objects. We also simulate variability due 
to the radio/optical observations being carried out non- 
simultaneously with the X-ray observations by adding 
random normal scatter to the simulated £ r /£ uv values 
(after calculating all interdependent quantities) with am- 
plitudes of 0.0792/0.114 (the same uncertainties as used 
with fitting; see §5), but inclusion of this minor effect 
does not significantly impact the modeling. 

Another relevant effect is the redshift-dependent effi- 
ciency of the SDSS color-selection targeting algorithm 
and of the photometric quasar classification algorithm. 
The median redshift of the simulated population is 
1.79/1.88 for the spectroscopic/photometric magnitude 
cutoffs, higher than that of the observed primary sample 
(for which the median redshift is 1.41). Although the red- 
shift dependencies of quasar color-selection techniques 
presumably explain most of the discrepancy between the 
observed and simulated median redshifts, there are a 
few more factors that could act to decrease the median 
redshift of the observed primary sample. For example, 
there is a bias toward lower redshifts for the minority of 
sources that were X-ray targeted (rather than serendip- 
itous) and a possible bias against spectroscopic identifi- 
cation of fainter objects for which measuring broad emis- 
sion lines is more difficult. Fortunately, our goal of eval- 
uating X-ray luminosity as a function of radio and opti- 
cal/UV luminosity does not depend on redshift (e.g., see 
§5.3), and so it is not necessary to account for all the 
various selection effects that may influence the difference 
in redshift distribution. We match in redshift and num- 
ber simultaneously by drawing from the large simulated 
population those objects with redshifts close to the ob- 
served sample. Note that the underlying basic model 
could be falsified at this stage were the simulated radio 
luminosities (specified primarily by the adopted luminos- 
ity function, the adopted intrinsic jet/lobe properties, 
and the random inclination) to disagree with observa- 
tion. However, a KS test gives a probability of p = 0.21, 
indicating that the simulated and observed samples can- 
not be considered to differ significantly in their radio lu- 
minosity distributions. The median radio luminosity for 
the primary sample is £ T = 32.83 with a standard de- 
viation of 0.89, while that of the simulated sample is 
£ r = 32.80 with a standard deviation of 0.79. It would 
also be possible at this stage to falsify the nature of the 
presumed correlation between intrinsic radio jet power 
and optical/UV disk luminosity if the distribution of 
simulated optical/UV luminosities were to disagree with 
observation (the median can always be made to agree 
by adjusting the ratio, but the distribution reflects the 
lack of beaming in the optical/UV luminosities, in con- 
trast to the radio). However, a KS test indicates that 
the simulated and observed samples are not inconsistent 
(p = 0.41) in their optical/UV luminosity distributions. 
The median optical/UV luminosity for the primary sam- 
ple is £ uv = 30.36 with a standard deviation of 0.57, 
while that of the simulated sample is £ uv — 30.33 with 
a standard deviation of 0.63. The radio-loudness dis- 
tributions are also consistent (KS probability p — 0.15) 
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Figure 19. Simulated populations of RIQs and RLQs compared with the observed primary sample. The models differ principally 
in the degree of beaming of the jet-linked X-ray emission: A is unbeamed, B is beamed with 7 X = 10.5 as for the radio jet, and 
C is beamed to a lesser degree with y x = 2.0. The left-hand side shows "excess" X-ray luminosity relative to comparable RQQs 
(see Figure 7a) for the observed sample (small black crosses) and for the simulation, with objects color-coded by inclination 
(blue is 5.8° < 6 < 10°, green is 10° < 9 < 20°, and red is 20° < 8 < 60°). The large black crosses give the median and 
interquartile range for the primary sample. The right-hand side shows joint 90% confidence ellipses for X-ray luminosity fit as 
a function of optical/UV and radio luminosity, for the primary sample (red), the full sample (black), and each model (blue). 
The dotted blue ellipses are for models CI. 5 and C3.0 with 7 X = 1.5 and 3.0. See §6.3 for discussion. 



with observed/simulated median values of 2.50/2.46 and 
standard deviations of 0.74/0.70. These radio and opti- 
cal/UV properties are the same for each of models A, B, 
and C. 

For any of models A, B, or C to constitute viable de- 
scriptions of the X-ray emission from RIQs and RLQs, 
it is necessary (but not sufficient) that the distribution 
of simulated X-ray luminosities not disagree with obser- 
vation. The X-ray disk/corona emission is completely 
specified by the optical/UV disk emission, and so it is 
the additional jet-linked X-ray emission that must be ac- 
curately modeled. The median value of the overall sim- 



ulated X-ray luminosity can always be made to agree 
with observations by adjusting the ratio between the jet- 
linked X-ray component and the intrinsic radio-jet lumi- 
nosity, but the distribution is strongly influenced by the 
degree to which the jet-linked emission is beamed. The 
log offsets between intrinsic radio-jet luminosity (prior to 
applying beaming) and intrinsic X-ray luminosity of the 
jet-linked component (prior to applying beaming to mod- 
els B and C) required to match simulated with observed 
median X-ray luminosities are —5.35, —7.25, and —7.25 
for models A, B, and C, respectively. Random log- normal 
scatter is added to these offsets; we find empirically that 
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a somewhat better match to the 4 scatter observed for 
the primary sample is obtained if the magnitude of this 
scatter is greater at lower simulated redshifts (0.7/0.3 
for z < 1/z > 1), but this could potentially be related to 
the increased fraction of X-ray bright targeted objects at 
lower redshifts within the primary sample (§5.4.2), so we 
do not ascribe physical significance to this result. The 
X-ray bulk Lorentz factor in model C is 7 X = 2.0; for 
reference, this corresponds to 6 K = 3.6 for an inclination 
of 9 = 5° (rather than 6 = 11.4 as with the 7 = 10.5 
that applies to the radio jet). For model C, X-ray bulk 
Lorentz factors of 7 X ~ 1 — 3 can also match the primary 
sample X-ray luminosity distributions, but as shown in 
§6.3, 7 X = 2.0 produces the best agreement to the ob- 
served 4(4v,4) relation. KS tests (treating upper lim- 
its as detections, but the simulated censoring matches 
the data) give probabilities of p = 0.37 and p = 0.30 for 
models A and C, respectively, indicating that the simu- 
lated and observed samples are not inconsistent in their 
X-ray luminosity distribution for these models. However, 
the KS test probability for model B is p = 0.0052; the 
simulated distribution of X-ray luminosities differs signif- 
icantly from that observed. In particular, a tail to large 
X-ray luminosity values in the simulated population for 
model B does not agree with the data, which appear to 
show at most a modest skew toward higher X-ray lumi- 
nosities (see also §5.4.3). This tail in the model is more 
pronounced at X-ray than radio frequencies because the 
dependence upon the 5 beaming factor is greater (by 1 — a 
in the exponent) for the X-ray emission. (Model B would 
be an even worse match to the primary sample if the 4 
medians were equalized.) These results may be seen in 
the left-hand panels of Figure 18. 

All three models obviously produce (by design) en- 
hanced X-ray emission relative to RQQs of comparable 
optical/UV luminosity, as can be seen in the right-hand 
panels of Figure 18 and as given in Table 10. However, 
the manner in which they do so differs. In these fig- 
ures the simulated objects are color-coded by inclination 
(blue for 5.8° < 9 < 10°, green for 10° < 6 < 20°, 
and red for 20° < 9 < 60°). For model A, there is es- 
sentially no inclination dependence to the X-ray excess, 
since the X-ray luminosity is unbeamed. For model B, 
there is a strong inclination dependence, with the low- 
est inclination RLQs being extremely X-ray bright (to 
a degree that does not match observations). Note that 
these simulated low-inclination objects are still expected 
to be disk-dominated at optical/UV frequencies (as also 
applies to models A and C) and are thus not blazars, 
so the comparison to the observed primary sample (from 
which all identified blazars have been removed) is prop- 
erly matched. For model C, there is a definite but modest 
dependence upon inclination. Simultaneous considera- 
tion of the distribution of radio, optical/UV, and X-ray 
luminosities, following the analysis performed in §4 and 
§5 for observed RIQs and RLQs, provides additional in- 
sight helpful to evaluating the feasibility of the three 
models. 

6.3. Modeling results 

We now assess the accuracy with which these models 
(generated as described in §6.1, §6.2, and Table 9) repro- 
duce the observed dependence of X-ray luminosity upon 
radio and optical/UV luminosities. We first examine how 



the "excess" X-ray luminosity (X-ray luminosity minus 
that of a RQQ of matched optical/UV luminosity) for 
the simulated RIQs and RLQs depends upon radio loud- 
ness for the three models. In Figure 19, left-hand side, 
the simulated objects are again color-coded by inclina- 
tion, and the primary sample data are plotted along with 
the median and interquartile range measurements from 
Figure 7a. It may be observed that model A does not 
appear to capture accurately the rise in X-ray brightness 
with increasing radio loudness. This is because in model 
A the X-ray luminosity is unbeamed, whereas the radio 
luminosity is beamed, and so the most radio-luminous 
and radio-loud objects tend to have low inclinations (are 
color-coded blue) but are not X-ray brighter than higher 
inclination objects. Conversely, model B overpredicts the 
degree to which X-ray brightness depends on increasing 
radio loudness, or inclination. It is only for model C that 
the X-ray brightness compared to RQQs may be seen to 
increase with radio loudness in a manner analogous to 
that observed. 

Modeling X-ray luminosity as a joint function of op- 
tical/UV and radio luminosity, as done for the observed 
RIQs and RLQs in §5.2, makes clear the differing man- 
ner in which the three models are dependent on opti- 
cal/UV and radio luminosity (see Table 10). The right- 
hand side of Figure 19 shows the joint 90% confidence 
ellipses for the primary (red) and for the full (black) 
observed samples, along with the calculated result for 
each model (blue). For model A, the X-ray luminos- 
ity is strongly dependent on optical/UV luminosity but 
only weakly dependent on radio luminosity. The best- 
fit relation for the population simulated with model A 
is 4 = -0.142 + 0.730 x £ uv + 0.100 x 4 This result 
reflects the underlying dependence of both the optical 
disk emission (and thereby the disk/corona X-ray com- 
ponent) and the "jet-linked" X-ray component on the in- 
trinsic (unbeamed) radio-jet luminosity. The parameter- 
ization of X-ray luminosity upon optical/UV and radio 
luminosity in model A is inconsistent with the fits to the 
primary and full samples. For model B, the X-ray lumi- 
nosity is strongly dependent upon radio luminosity but 
only weakly dependent on optical/UV luminosity. This 
result is due to the increased dominance of the X-ray 
continuum by jet-linked emission for a large fraction of 
the simulated sources, particularly those at low inclina- 
tions. The best-fit relation for the population simulated 
with model B is 4 = 0.099 + 0.195 x £ uv + 0.746 x 4 
This reflects the mutual beaming of the radio and X-ray 
jets, but these parameters are also inconsistent with ob- 
servations. For model C, the X-ray luminosity is de- 
pendent upon both optical/UV and radio luminosities 
in a manner similar to that observed, as illustrated 
by the overlapping confidence ellipses. The best-fit re- 
lation for the population simulated with model C is 
4 = -0.098 + 0.526 x 4 v + 0.279 x 4 This result reflects 
the milder beaming with 7 X = 2.0 that places model C 
in some sense intermediate between model A and model 
B. For comparison, model C with -f x = 1.5/3.0 (models 
C1.5/C3; best-fit relations given in Table 10 and shown 
as dotted blue ellipses in Figure 19) bracket the 7 X = 2.0 
result but do not provide as close a match to observa- 
tions. 

The simulated populations may be utilized to evalu- 
ate the dominant emission component at each frequency, 
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Figure 20. Simulated populations of RIQs and RLQs. The 
models differ principally in the degree of beaming of the jet- 
linked X-ray emission: A is unbeamed, B is beamed with 
7 X = 10.5 as for the radio jet, and C is beamed to a lesser 
degree with 7 X = 2.0. In all panels objects are color-coded by 
inclination (blue is 5.8° < 9 < 10°, green is 10° < 9 < 20°, 
and red is 20° < 9 < 60°). The left-hand side shows X-ray 
luminosity components versus optical/UV luminosity com- 
ponents, while the right-hand side shows X-ray luminosity 
components versus radio luminosity components. 

and in particular to examine the degree to which X-ray 
luminosity is dominated by the jet-linked component, as 
a function of inclination, for the different models (see 
Figure 20). In nearly all cases, the optical/UV lumi- 
nosity is disk-dominated. For most sources, and for all 
the low-inclination sources, the radio luminosity is jet- 
dominated. In each of models A, B, and C, the X-ray 
luminosity is dominated by disk/corona emission for a 
sizeable fraction of objects; there is agreement that de- 
spite the greater median X-ray luminosities of RIQs and 
RLQs relative to RQQs, there are many individual ob- 
jects for which a disk/corona analogous to that in RQQs 
is sufficient to produce the majority of observed X-ray 
emission. For models B and C, the degree to which the 
X-ray emission is jet-dominated is correlated with the 
degree to which the radio emission is jet-dominated (to 
a lesser degree in model C, for which the beaming fac- 
tor is lower). For the preferred model C, the majority of 
objects with inclinations less than 20° have X-ray lumi- 
nosities dominated by jet-linked emission. 

Despite the general nature of the scenarios considered 
under models A, B, and C, these simulations provide 
insight into the degree to which the "excess" X-ray lu- 
minosity in RIQs and RLQs is jet-linked and beamed. 
It seems likely that some beaming must be present, as 
model A does not satisfactorily match observations; this 
conclusion is consistent with the results of previous work 
(e.g., Browne & Murphy 1987; Shastri et al. 1993). It 
also appears that the jet-linked X-ray emission is beamed 
to a lesser degree than affects the radio jet, as model 
B does not satisfactorily match observations; this re- 



sult has also been suggested by other work (e.g., Dou 
& Yuan 2008). The success of model C suggests that the 
majority of the jet-linked X-ray emission is not produced 
cospatially with the radio jet emission, but likely primar- 
ily originates in a somewhat slower jet region. The pre- 
ferred model C implies that high signal-to-noise X-ray 
spectroscopy of RIQs or RLQs established (e.g., from 
radio properties) as possessing intermediate inclinations 
should result in some such objects showing X-ray spectra 
in which an unabsorbed disk/corona-linked component 
and a jet-linked component are simultaneously apparent. 
This is consistent with observational results (although 
additional X-ray studies would help establish this point); 
for example, Sambruna et al. (2006) describe three RLQs 
with concave broken power law X-ray spectra, which they 
attribute to dual jet/disk-corona components. 

Future work with significantly larger samples and 
denser coverage of the I — z plan^HH could potentially 
investigate whether the small-scale jet has velocity struc- 
ture (as in a fast-spine/slow-sheath model; see, e.g., 
Jester et al. 2006b), most simply through considering 
two X-ray jet components with differing bulk Lorentz 
factors. Modeling deficiencies that are unlikely to im- 
pact these results but could be addressed to improve the 
precision of comparison to data include limited consid- 
eration of selection effects and ad-hoc normalization of 
simulated sample size. 

6.4. Alternative model 'parameters 

Although we have not conducted an exhaustive search 
of the multidimensional parameter space, the presented 
models (A, B, C, CI. 5, and C3.0) together cover a large 
range of plausible possibilities, and from this set model 
C best matches observations. It is of interest to eval- 
uate whether a modified set of models would still yield 
qualitatively similar results. Here, we briefly describe 
the effects of altering selected model parameters in a 
reasonable fashion. Minor adjustments to the log off- 
sets between intrinsic radio-jet luminosity and intrinsic 
jet-linked X-ray luminosity are generally made without 
comment below; other model parameters are as given in 
§6.1, §6.2, and Table 9, unless otherwise stated. 

The credible interval for 7 from the adopted Mullin 
& Hardcastle (2009) model ranges from 3.20 to 14.05. 
The higher 7 produces an acceptable match to the pri- 
mary sample radio and optical/UV luminosities (KS 
probabilities of 0.59/0.26/0.13 for t r /i uv /R*) and can 
match the X-ray luminosity of models A and C but 
not B (KS probabilities of 0.21/< 0.01/0.49 for A/B/C). 
The lower 7 does not easily give a good match to both 
the radio and optical/UV luminosities. For a ratio of 
optical/UV disk emission to unbeamed radio emission 
of —2.0 (logarithmic units) and a minimum inclination 
of 4°, the KS probabilities are 0.059/0.075/0.0076 for 
£ r /£ uv / R* . The simulated median X-ray luminosities 
may be matched to observations with resulting KS prob- 
abilities of 0.35/< 0.01/0.10 for models A/B/C. In this 
scenario the majority of the simulated objects have in- 
clinations with 9 > 20°, but the X-ray brightness of 
the lower-inclination objects cause the increase in "ex- 

37 We are currently carrying out Swift snapshot observations of 
high-luminosity RLQs (Jianfeng Wu et al., in preparation) that 
will help fill in the I — z plane. 
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cess" X-ray luminosity with increasing radio-loudness for 
model B and (to a lesser degree) model C to be more 
rapid than observed. For model C under this scenario, 
the dependence of X-ray luminosity upon optical/UV 
and radio luminosity is weaker and stronger (b uv = 0.31 
and c r = 0.55), respectively, than holds for the primary 
sample. 

Willott et al. (2001) consider two additional models 
for their radio luminosity function in addition to the 
version we use (see their §3.2). If instead a population 
of RIQs and RLQs is simulated based on their model 
which contains a steeper high-redshift decay (their Equa- 
tion 11), then an acceptable match to the primary sample 
may be obtained (KS probabilities of 0.24/0.12/0.42 for 
&t ft-iiv / R*) with a ratio of optical/UV disk emission to 
unbeamed radio emission of —1.9. The KS probabili- 
ties for matching simulated versus observed X-ray lumi- 
nosity are 0.065/< 0.01/0.35 for models A/B/C. Model 
B is again not a good match. If instead the Willott et 
al. (2001) model which becomes constant at high redshift 
(their Equation 12) is used to simulate a population of 
RIQs and RLQs, the match to the primary sample is 
again acceptable (KS probabilities of 0.082/0.19/0.26 for 
ii/iuv/R*)- The KS probabilities for matching simulated 
versus observed X-ray luminosity are 0.10/< 0.01/0.20 
for models A/B/C. Model B remains a poor match. Ad- 
ditionally, the median simulated redshift (prior to match- 
ing redshifts to the primary sample) is 1.95/2.22 for 
the spectroscopic/photometric magnitude cutoffs, higher 
than that of the primary sample (1.41) and also higher 
than that simulated with the version of the Willott et 
al. (2001) model we use (1.79/1.88). 

We also consider a scenario in which the X-ray jet emis- 
sion is beamed as 5 2_Qx (similarly to the radio jet emis- 
sion) rather than as S 3 ~ 2a *. This might apply if, for 
example, the X-ray jet emission were to be dominated 
by synchrotron emission from a secondary population of 
high-energy electrons. This change only affects models B 
and C. Maintaining a x = —0.3/ — 0.5 and adjusting the 
X-ray jet offsets to —6.45/ — 6.50, the KS probabilities 
are 0.065/0.76 for models B/C. With less extreme X-ray 
beaming model B can match the observed distribution 
of X-ray luminosities. However, model B still predicts a 
greater fraction of X-ray-bright objects than is observed, 
and the increase in "excess" X-ray luminosity with in- 
creasing radio-loudness is also more rapid than observed. 
Model C remains a superior explanation of the data. For 
model C under this scenario, the dependence of X-ray lu- 
minosity upon optical/UV and radio luminosity is mod- 
erately stronger and weaker (b uv = 0.58 and c r = 0.23), 
respectively, than holds for the primary sample. 

7. SUMMARY 

The primary results of our analysis of the nuclear X-ray 
properties of RIQs and RLQs are the following: 

1. Sample: We have compiled a sample of 188 RIQs 
and 603 RLQs (primarily by matching optically selected 
SDSS quasars to the FIRST survey, taking into account 
extended radio emission) with high-quality archival X- 
ray coverage by Chandra, XMM-Newton, or ROSAT. 
The full sample is almost unbiased with respect to X-ray 
properties, has a high (85%) detection rate, and can be 
utilized to investigate the nature and origin of X-ray 
emission in RIQs and RLQs. The sample size is signifi- 



cantly larger than earlier studies of RLQs and provides 
superior coverage of the luminosity-redshift plane. 

2. Trends: We calculate the ratio of X-ray luminosity 
in RIQs and RLQs relative to that of comparable RQQs 
and determine how this ratio increases with both radio 
loudness and luminosity. This ratio of "excess" X-ray 
luminosity ranges from ~0.7-2.8 for RIQs through the 
canonical ~3 for RLQs to > 10 for strongly radio-loud 
(R* > 4) or luminous (£ T > 35) objects. We also carry 
out fits to the X-ray luminosity as a function of both 
optical and radio luminosity, which are useful for deter- 
mining the "typical" X-ray luminosity for any given RIQ 
or RLQ. The X-ray emission in RIQs is not significantly 
dependent upon radio luminosity. We quantify the man- 
ner in which RLQs become more X-ray luminous for large 
R* or l t or for flat radio spectra. Finally, we do not find 
any significant redshift dependence in the properties of 
RIQs and RLQs (implying, e.g., that IC/CMB jet-linked 
emission does not contribute substantially to the nuclear 
X-ray continuum). 

3. Models: We conduct Monte Carlo simulations based 
on a low-frequency radio luminosity function to which 
we add a randomly inclined relativistic jet. The opti- 
cal disk emission is successfully modeled as correlated 
with the intrinsic (unbeamed) radio-jet luminosity, and 
the X-ray emission contains both disk/corona and jet- 
linked components. We consider three models for the 
non-coronal core X-ray luminosity component, and con- 
clude that the jet-linked X-ray emission is likely beamed 
but with a lesser bulk Lorentz factor than applies to 
the radio-jet emission. The alternative possibilities of 
unbeamed X-ray emission and of X-ray emission with 
7 ~ 10.5 as for the radio jet do not appear to match ade- 
quately the observed data. For the preferred model, the 
radio emission is mostly jet dominated, the optical/UV 
emission is almost exclusively disk dominated, and the X- 
ray emission is split between disk/corona and jet-linked 
components with the jet becoming increasingly dominant 
at low inclinations. 

We thank the anonymous referee for many construc- 
tive suggestions that improved this paper. We thank 
B. Luo, D. A. Rafferty, and Y. Xue for assistance with 
photometric redshifts and for making their catalogs avail- 
able prior to publication. We thank M. Eracleous and 
B. Kelly for useful discussions. BPM, WNB, and DPS 
thank NASA ADP grant NNX10AC99G for support; 
DPS also thanks NSF grant AST06-07634. Funding for 
the SDSS and SDSS-II has been provided by the Alfred 
P. Sloan Foundation, the Participating Institutions, the 
National Science Foundation, the U.S. Department of 
Energy, the National Aeronautics and Space Adminis- 
tration, the Japanese Monbukagakusho, the Max Planck 
Society, and the Higher Education Funding Coun cil fo r 
England. The SDSS Web Site is |http://www.sdss.org/| 

REFERENCES 

Adelman-McCarthy, J. K., et al. 2007, ApJS, 172, 634 

Alexander, D. M., et al. 2003, AJ, 126, 539 

Antonucci, R. 1993, ARA&A, 31, 473 

Avni, Y., & Tananbaum, H. 1986, ApJ, 305, 83 

Barthel, P. D. 1989, ApJ, 336, 606 

Baldwin, J. A. 1977, ApJ, 214, 679 

Barger, A. J., et al. 2003, AJ, 126, 632 



29 



Barger, A. J., Cowie, L. L., & Wang, W.-H. 2007, ApJ, 654, 764 
Barger, A. J., Cowie, L. L., & Wang, W.-H. 2008, ApJ, 689, 687 
Bassett, L. C, Brandt, W. N., Schneider, D. P., Vignali, C, 

Chartas, G., & Garmire, G. P. 2004, AJ, 128, 523 
Becker, R. H., White, R. L., & Helfand, D. J. 1995, ApJ, 450, 559 
Becker, R. H., White, R. L., Gregg, M. D., Brotherton, M. S., 

Laurent-Muehleisen, S. A., & Arav, N. 2000, ApJ, 538, 72 
Becker, R. H., et al. 2001, ApJS, 135, 227 
Bell, M. B., & Comeau, S. P. 2010, Ap&SS, 325, 31 
Belsole, E., Worrall, D. M-, & Hardcastle, M. J. 2006, MNRAS, 

366, 339 

Best, P. N., Kauffmann, G., Heckman, T. M., & Ivezic, Z. 2005, 

MNRAS, 362, 9 
Biggs, A. D., & Ivison, R. J. 2006, MNRAS, 371, 963 
Blundell, K. M., & Kuncic, Z. 2007, ApJ, 668, L103 
Brandt, W. N., & Alexander, D. M. 2010, Proceedings of the 

National Academy of Science, 107, 7184 
Brinkmann, W., Yuan, W., & Siebert, J. 1997, A&A, 319, 413 
Brotherton, M. S., Laurent-Muehleisen, S. A., Becker, R. H., 

Gregg, M. D., Telis, G., White, R. L., & Shang, Z. 2005, AJ, 

130, 2006 

Browne, I. W. A., & Murphy, D. W. 1987, MNRAS, 226, 601 
Capak, P., et al. 2007, ApJS, 172, 99 
Cappelluti, N., et al. 2009, A&A, 497, 635 

Cappi, M., Matsuoka, M., Comastri, A., Brinkmann, W., Elvis, 
M., Palumbo, G. G. C, & Vignali, C. 1997, ApJ, 478, 492 

Chiaberge, M., Capetti, A., & Celotti, A. 2000, A&A, 355, 873 

Cirasuolo, M., Celotti, A., Magliocchetti, M., & Danese, L. 2003, 
MNRAS, 346, 447 

Coleman, G. D., Wu, C.-C, & Weedman, D. W. 1980, ApJS, 43, 
393 

Condon, J. J., Cotton, W. D., Greisen, E. W., Yin, Q. F., Perley, 

R. A., Taylor, G. B., & Broderick, J. J. 1998, AJ, 115, 1693 
Delia Ceca, R., et al. 2004, A&A, 428, 383 
Dermer, C. D. 1995, ApJ, 446, L63 

Donoso, E., Best, P. N., & Kauffmann, G. 2009, MNRAS, 392, 
617 

Dou, L.-M., & Yuan, W.-M. 2008, Chinese Journal of Astronomy 

and Astrophysics, 8, 653 
Douglas, J. N., Bash, F. N., Bozyan, F. A., Torrence, G. W., & 

Wolfe, C. 1996, AJ, 111, 1945 
Elvis, M., et al. 2009, ApJS, 184, 158 

Evans, D. A., Worrall, D. M., Hardcastle, M. J., Kraft, R. P., & 

Birkinshaw, M. 2006, ApJ, 642, 96 
Evans, I. N., Evans, J. D., Fabbiano, G., Glotfelty, K. J., 

McCollough, M. L., McDowell, J. C, Primini, F. A., & Rots, 

A. H. 2008, Proceedings of the SPIE, 7016, 53 
Falcke, H., Sherwood, W., & Patnaik, A. R. 1996, ApJ, 471, 106 
Fanaroff, B. L., & Riley, J. M. 1974, MNRAS, 167, 31P 
Galbiati, E., et al. 2005, A&A, 430, 927 

Garofalo, D., Evans, D. A., & Sambruna, R. M. 2010, MNRAS, 
406, 975 

Gawiser, E., et al. 2006, ApJS, 162, 1 
Gehrels, N. 1986, ApJ, 303, 336 
Giacconi, R., et al. 2002, ApJS, 139, 369 

Gibson, R. R., Brandt, W. N., & Schneider, D. P. 2008, ApJ, 685, 
773 

Gibson, R. R., et al. 2009, ApJ, 692, 758 
Gopal-Krishna, & Wiita, P. J. 2000, A&A, 363, 507 
Green, P. J., et al. 2009, ApJ, 690, 644 

Gregory, P. C, Scott, W. K., Douglas, K., & Condon, J. J. 1996, 

ApJS, 103, 427 
Haardt, F., & Maraschi, L. 1991, ApJ, 380, L51 
Hall, P. B., Gallagher, S. C, Richards, G. T., Alexander, D. M., 

Anderson, S. F., Bauer, F., Brandt, W. N., & Schneider, D. P. 

2006, AJ, 132, 1977 
Hardcastle, M. J., Evans, D. A., & Croston, J. H. 2009, MNRAS, 

396, 1929 
Hasinger, G. 2008, A&A, 490, 905 

Ilbert, O., et al. 2008, Astronomical Society of the Pacific 

Conference Series, 399, 169 
Ilbert, O., et al. 2009, ApJ, 690, 1236 
Isobe, T., & Feigelson, E. D. 1990, BAAS, 22, 917 
Ivezic, Z., et al. 2002, AJ, 124, 2364 

Ivezic, Z., et al. 2004, Astronomical Society of the Pacific 

Conference Series, 311, 347 
Jackson, C. A., & Wall, J. V. 1999, MNRAS, 304, 160 



Jester, S., et al. 2005, AJ, 130, 873 

Jester, S., Kording, E., & Fender, R. 2006a, VI Microquasar 

Workshop: Microquasars and Beyond, 1, 26 
Jester, S., Harris, D. E., Marshall, H. L., & Meisenheimer, K. 

2006b, ApJ, 648, 900 
Jiang, L., Fan, X., Ivezic, Z., Richards, G. T., Schneider, D. P., 

Strauss, M. A., & Kelly, B. C. 2007, ApJ, 656, 680 
Just, D. W., Brandt, W. N., Shemmer, O., Steffen, A. T., 

Schneider, D. P., Chartas, G., & Garmire, G. P. 2007, ApJ, 

665, 1004 

Kaspi, S., Brandt, W. N., Maoz, D., Netzer, H., Schneider, D. P., 

& Shemmer, O. 2007, ApJ, 659, 997 
Kellermann, K. I., Sramek, R., Schmidt, M., Shaffer, D. B., & 

Green, R. 1989, AJ, 98, 1195 
Kellermann, K. I., Fomalont, E. B., Mainieri, V., Padovani, P., 

Rosati, P., Shaver, P., Tozzi, P., & Miller, N. 2008, ApJS, 179, 

71 

Kelly, B. C. 2007, ApJ, 665, 1489 

Kelly, B. C, Bechtold, J., Siemiginowska, A., Aldcroft, T., & 

Sobolewska, M. 2007, ApJ, 657, 116 
Kraft, R. P., Burrows, D. N., & Nousek, J. A. 1991, ApJ, 374, 344 
Labita, M., Treves, A., & Falomo, R. 2008, MNRAS, 383, 1513 
Landt, H., Padovani, P., Giommi, P., Perri, M., & Cheung, C. C. 

2008, ApJ, 676, 87 
Laor, A. 2000, ApJ, 543, Llll 
Laor, A., & Behar, E. 2008, MNRAS, 390, 847 
Lehmer, B. D., et al. 2005, ApJS, 161, 21 
Lonsdale, C. J., et al. 2003, PASP, 115, 897 
Lopez, L. A., Brandt, W. N., Vignali, C, Schneider, D. P., 

Chartas, G., Garmire, G. P. 2006, AJ, 131, 1914 
Lu, Y., Wang, T., Zhou, H., & Wu, J. 2007, AJ, 133, 1615 
Lu, Y., Wang, T., Zhou, H., & Wu, J. 2010. larXivTT0"0 2.0636 
Luo, B., et al. 2008, ApJS, 179, 19 
Luo, B., et al. 2009, ApJ, 695, 1227 
Luo, B., et al. 2010, ApJS, 187, 560 
Lynden-Bell, D. 1969, Nature, 223, 690 
Mackay, C. D. 1971, MNRAS, 154, 209 
Maoz, D. 2007, MNRAS, 377, 1696 
Marshall, H. L., et al. 2005, ApJS, 156, 13 
Martin, D. C, et al. 2005, ApJ, 619, LI 

McGreer, I. D., Helfand, D. J., & White, R. L. 2009, AJ, 138, 
1925 

Meier, D. L. 2001, ApJ, 548, L9 
Menou, K., et al. 2001, ApJ, 561, 645 

Metcalf, R. B., & Magliocchetti, M. 2006, MNRAS, 365, 101 
Mignoli, M., et al. 2005, A&A, 437, 883 

Miller, P., Rawlings, S., & Saunders, R. 1993, MNRAS, 263, 425 
Miller, N. A., Fomalont, E. B., Kellermann, K. I., Mainieri, V., 

Norman, C, Padovani, P., Rosati, P., & Tozzi, P. 2008, ApJS, 

179, 114 

Miller, B. P., Brandt, W. N., Gibson, R. R., Garmire, G. P., & 

Shemmer, O. 2009, ApJ, 702, 911 
Mullin, L. M., & Hardcastle, M. J. 2009, MNRAS, 398, 1989 
Norris, R. P., et al. 2006, AJ, 132, 2409 
O'Dea, C. P. 1998, PASP, 110, 493 

Orr, M. J. L., & Browne, I. W. A. 1982, MNRAS, 200, 1067 
Osmer, P. S., & Shields, J. C. 1999, Astronomical Society of the 

Pacific Conference Series, 162, 235 
Padovani, P., Allen, M. G., Rosati, P., & Walton, N. A. 2004, 

A&A, 424, 545 

Page, K. L., Reeves, J. N., O'Brien, P. T., & Turner, M. J. L. 

2005, MNRAS, 364, 195 
Plotkin, R. M., Anderson, S. F., Hall, P. B., Margon, B., Voges, 

W., Schneider, D. P., Stinson, G., & York, D. G. 2008, AJ, 135, 

2453 

Plotkin, R. M., et al. 2010, AJ, 139, 390 

Rafferty, D. A., Brandt, W. N., Alexander, D.M. , Xue, Y. Q., 

Bauer, F. E., Lehmer, B. D., Luo, B., & Papovich, C. 2010, 

ApJ, submitted 
Rawlings, S., & Saunders, R. 1991, Nature, 349, 138 
Rees, M. J., & Setti, G. 1968, Nature, 219, 127 
Rengelink, R. B., Tang, Y., de Bruyn, A. G., Miley, G. K., 

Bremer, M. N., Roettgering, H. J. A., & Bremer, M. A. R. 

1997, A&AS, 124, 259 
Richards, E. A., Kellermann, K. I., Fomalont, E. B., Windhorst, 

R. A., & Partridge, R. B. 1998, AJ, 116, 1039 



30 



Richards, E. A., Fomalont, E. B., Kellermann, K. I., Windhorst, 

R. A., Partridge, R. B., Cowie, L. L., & Barger, A. J. 1999, 

ApJ, 526, L73 
Richards, G. T., et al. 2001, AJ, 121, 2308 
Richards, G. T., et al. 2002, AJ, 123, 2945 
Richards, G. T., et al. 2009, ApJS, 180, 67 
Rovilos, E., Gcorgakakis, A., Georgantopoulos, I., Afonso, J., 

Koekemoer, A. M., Mobasher, B., & Goudis, C. 2007, A&A, 

466, 119 

K., Maraschi, 
Urry, C. M.. 



L., Tavecchio, F., 
Chartas, G. 2004, 



Tavecchio, F., Maraschi, L., & 



& de 



1999, 



Foltz, C. B. 1992, 



Sambruna, R. M., Gambill, J. 
Cerutti, R., Cheung, C. C. 
ApJ, 608, 698 
Sambruna, R. M., Gliozzi, M. 

Foschini, L. 2006, ApJ, 652, 146 
Scheuer, P. A. G. 1995, MNRAS, 277, 331 
Schinnerer, E., et al. 2007, ApJS, 172, 46 
Schneider, D. P., et al. 2007, AJ, 134, 102 
Schwartz, D. A. 2002, ApJ, 569, L23 

Shankar, F., Dai, X., & Sivakoff, G. R. 2008, ApJ, 687, 859 
Shankar, F., Sivakoff, G. R., Vestergaard, M., & Dai, X. 2010, 

MNRAS, 401, 1869 
Shastri, P., Wilkes, B. J., Elvis, M., & McDowell, J. 1993, ApJ, 

410, 29 

Shields, G. A. 1978, Nature, 272, 706 

Siemiginowska, A., LaMassa, S., Aldcroft, T. L., Bechtold, J., & 

Elvis, M. 2008, ApJ, 684, 811 
Sokolov, A., & Marscher, A. P. 2005, ApJ, 629, 52 
Spergel, D. N., et al. 2007, ApJS, 170, 377 
Steffcn, A. T., Strateva, I., Brandt, W. N., Alexander, D. M., 

Koekemoer, A. M., Lehmer, B. D., Schneider, D. P., & Vignali, 

C. 2006, AJ, 131, 2826 
Stocke, J. T., Morris, S. L., Weymann, R. J. 

ApJ, 396, 487 

Strateva, I. V., Brandt, W. N., Schneider, D. P., Vanden Berk, 

D. G., & Vignali, C. 2005, AJ, 130, 387 
Suchkov, A. A., Hanisch, R. J., Voges, W., & Heckman, T. M. 

2006, AJ, 132, 1475 
Szokoly, G. P., et al. 2004, ApJS, 155, 271 
Tavecchio, F., Maraschi, L., Sambruna, R. M., & Urry, C. M. 

2000, ApJ, 544, L23 
Tavecchio, F., Maraschi, L., Ghisellini, G., Kataoka, J., Foschini, 

L., Sambruna, R. M., & Tagliaferri, G. 2007, ApJ, 665, 980 
Tozzi, P., et al. 2009, ApJ, 698, 740 

Trammell, G. B., Vanden Berk, D. E., Schneider, D. P., Richards, 
G. T., Hall, P. B., Anderson, S. F., & Brinkmann, J. 2007, AJ, 
133, 1780 

Trouille, L., Barger, A. J., Cowie, L. L., Yang, Y., & Mushotzky, 
R. F. 2008, ApJS, 179, 1 

A. DETAILS OF PRIMARY SAMPLE RADIO/OPTICAL/X-RAY MATCHING PROCESS 

We here describe the details of the selection, screening, and matching process used to generate the spectroscopic sample. Many 
of these steps are identically used in the construction of the photometric sample. Together, the spectroscopic and photometric 
samples make up the primary sample of RIQs and RLQs used throughout this paper. 

The SDSS quasar survey is ~95% complete for unresolved sources (for objects categorized as quasars on the basis of optical 
spectra) and ~89% complete overall (Vanden Berk et al. 2005); the DR5 quasar catalog of Schneider et al. (2007) consists of 
quasars with strong broad emission lines (of width >1000 km s _1 ) and by design does not include BL Lacs with featureless 
spectra (e.g., Plotkin et al. 2008, 2010). It does, however, include broad-absorption line (BAL) quasars, and it is now clear (e.g., 
Becker et al. 2000, 2001; Menou et al. 2001) that a small fraction (<8% at high radio luminosities; e.g., Shankar et al. 2008) of 
RLQs have BALs and are, like BAL RQQs, relatively X-ray weak due to intrinsic absorption (e.g., Brotherton et al. 2005; Gibson 
et al. 2008; Miller et al. 2009). Since we are interested in intrinsic emission and not absorption, we exclude from our sample any 
objects included in the SDSS BAL catalog of Gibson et al. (2009) and also any RLQs we are aware of as possessing BALs not 
observable in the SDSS spectrum (e.g., J100726.10+124856.2, the z = 0.24 object PG 1004+130 identified as a BAL RLQ by 
Wills et al. 1999). Through this process 22 BAL quasars are identified and dropped from further consideration. Contamination 
from low-redshift (z < 1.7) RLQs not identifiable as having BALs due to SDSS bandpass limitations (e.g., spectroscopic quasars 
lacking coverage of the 1400-1550 A region, which includes the strong C IV BAL transition) is estimated to be ~2% of the entire 
sample[ff| and a small number of X-ray weak outliers would not materially impact our results below. Objects with apparent 
heavy intrinsic reddening are sometimes also associated with X-ray absorption (e.g., Hall et al. 2006). The quantity A(g — i) 
(plotted in Figure 2) is the g — i color of an object (corrected for Galactic extinction) less the median g — i quasar color at the 
redshift of that object; negative A(g — i) values correspond to objects bluer than the norm and positive A(g — i) values indicate 
redder colors (e.g., Richards et al. 2001). For reference, the median relative color for our sample of spectroscopic SDSS RIQs 
and RLQs is A(g — i) = 0.07 (there is a general tendency for RLQs to have somewhat redder colors than RQQs; e.g., Labita et 
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al. 2008). Almost all non-reddened objects have A(g — i) < 1 while intrinsically reddened (and potentially absorbed) objects 
form a tail in the distribution that extends past this value (e.g., Hall et al. 2006), so objects with A(g — i) > 1 are culled from 
the sample. This cut removes 8 objects that would otherwise be included (an additional 4 of the already excluded BAL quasars 
also have A(g — i) > 1). Finally, 3 objects possess radio spectra (described in §3.3) indicating they are GHz-peaked spectrum 
(GPS; e.g., O'Dea 1998) sources, which are generally considered to be young sources and have X-ray properties that are not 
necessarily representative of RLQs in general (e.g., Siemiginowska et al. 2008); these GPS sources are therefore culled and not 
considered further. 

Extended radio sources (with jets or lobes) may have separate entries in the FIRST source catalog for each detected component, 
and lobe-dominated sources do not always possess a detected core component. It is therefore necessary to match to the entire 
environment surrounding each quasar to recover all associated emission (e.g., Best et al. 2005; Lu et al. 2007). Very few radio 
sources that lack a core component have lobe-to-lobe angular sizes greater than 90" (e.g., de Vries et al. 2006), so we search for 
radio emission within 90" of each optically-selected quasar. We considered as candidates for radio-detected quasars those objects 
with a radio source within 2", referred to hereafter as the core, or two or more radio sources within 2" < 6 < 90" of which 
at least one pair has component angular orientations relative to the core differing by more than 90° (where a 180° separation 
would indicate components on directly opposite sides of the core), classified as potential lobes pending additional review. The 
potential lobes were then matched to the SDSS photometric catalog, and if they had an optical counterpart (unless that apparent 
counterpart was a spectroscopically-classified star) they were flagged as intruding background sources and not considered further 
in determining the radio characteristics of that quasar; if this process left only one potential lobe and no core then the object 
was eliminated as not radio-detected. All objects with remaining potential lobes were then examined visually to screen further 
for misclassifications (for example, unrelated background double-lobed radio sources are readily identifiable by eye as intruding 
but are difficult to categorize correctly automatically). A few cases in which an obviously-associated extended component had 
been mistakenly flagged as background due to a randomly coincident SDSS optical source were also corrected. Visual screening 
also identified four objects with likely associated radio components outside the 90" consideration radius (J093200. 08+553347. 4, 
J112956.53+364919.2, J094745.14+072520.6, and J142735.60+263214.5, with maximal lobe offset from the core of 93", 99", 
117", and 148", respectively); since we wish to take into account all lobe emission from objects already included in the sample, 
these components were added to the total radio fluxes for these objects. Visual screening also confirmed that the FIRST catalog 
contains entries for all obvious components in all sources included in the sample [with only one exception: we added a core 
component to J170441. 37+604430.5 (3C 351) based on examination of the FIRST and other radio maps]. The FIRST survey is 
not particularly sensitive to large-scale regions of low surface brightness, so diffuse lobes may not be detected or may have their 
fluxes underestimated (e.g., Lu et al. 2010). The NRAO VLA Sky Survey (NVSS; Condon et al. 1998) provides lower angular 
resolution coverage that complements the high angular resolution imaging of FIRST; 1.4 GHz fluxes from the NVSS catalog 
were used for a handful of obviously extended FIRST sources. 

The 312 RIQs and RLQs in the spectroscopic sample all possess sensitive X-ray coverage from Chandra, XMM-Newton, 
or ROSAT . For Chandra, we checked all public non-grating ACIS-S or ACIS-I observations with exposures longer than 1 ks 
and off-axis angles of less than 12' (note that each one of the ACIS CCDs covers an area of approximately 8.4' by 8.4'). For 
XMM-Newton, we considered only those observations with exposure times greater than 1 ks and off-axis angles of less than 15'. 
For ROSAT , we initially considered unfiltered PSPC observations with exposure times greater than 2 ks and off-axis angles less 
than 40'. All X-ray observations were examined to screen out cases in which the quasar did not fall on the detector, was too 
close to another bright source, was located within an instrumental artifact, or was otherwise unsuitable for further analysis. 
Due to their greater X-ray sensitivity, Chandra or XMM-Newton observations were prioritized over ROSAT observations. In 
cases where both Chandra and XMM-Newton provide coverage, the observation with the greatest X-ray signal-to-noise ratio 
(estimated based on exposure time, off-axis angle, detector efficiency, and average background) was used. ROSAT observations 
were ranked by X-ray signal-to-noise ratio, and lower-quality observations were discarded. The threshold value for discarding 
ROSAT observations was chosen so as to provide a large sample size while maintaining a relatively high detection fraction. In 
cases where multiple ROSAT observations were available, the observation with the greatest X-ray signal-to-noise ratio was used. 
Despite their lower sensitivity, the ROSAT observations are useful for this project, as they provide a large area of sky coverage 
beyond that available solely from Chandra and XMM-Newton observations. Note that this method is unbiased with respect to 
intrinsic quasar properties (including X-ray brightness) as it is based only on the expected quality of the X-ray observations, so 
it should not lead to any biases or systematic errors in our analyses. 

B. USING GALEX DATA TO IMPROVE PHOTOMETRIC REDSHIFTS 

We here describe the process through which GALEX observations were used to identify and discard 34 RIQs and RLQs with 
likely incorrect photometric redshifts and a further 25 objects for which the accuracy of the photometric redshifts could not be 
checked (see also §2.1.2). 

Out of the 406 photometric RIQs and RLQs under consideration prior to incorporation of GALEX data, 151 (37%) have 
Zphot > 1-9, of which 121 have GALEX coverage. From these 121, we use a color-color cut and a joint-detection cut to identify 
43 objects with questionable photometric redshifts. For these 43 objects, 37/29/23 (86%/67%/53%) are flagged by the color- 
color/joint-detection/both methods (so 14/6, or 33%/14%, are flagged exclusively by the color-color/joint-detection method). 
From these 43 objects, 9 have spectroscopic redshifts (in all cases confirming the photometric redshift was incorrect) and are 
retained with corrected redshifts and luminosities. The remaining 34 objects are discarded. Of the 30 objects with z p hot > 1.9 
lacking GALEX coverage, 5 have spectroscopic redshifts (in all cases indicating the photometric redshift is correct); the remaining 
25 objects are conservatively discarded, sacrificing some completeness to minimize potential contamination. Additional details 
are provided below. 

The color-color index (winuv — rn T ) — 2.5(m g — m r ) cleanly separates 2 spcc < 1 from z spac > 1.9 objects (D. W. Hogg 2009, 
personal communication). We flag the 37 objects with photometric redshifts z p hot > 1-9 but (w-nuv — rn r ) — 2.5(m g — m r ) < 1.5, 
of which we estimate ~31 are actually low-redshift quasars. From the catalog of SDSS DR3 quasars matched to GALEX data 
provided by Trammell et al. (2007), of the quasars with both FUV and NUV detections the fraction with z spcc > 1.9 is 6.9% 
(or 7.7%/8.6% for objects with m\ > 19/20; the median m\ for the photometric sample is 20.0). In contrast, within the primary 
sample, ~23% of the quasars with both FUV and NUV detections have z p hot > 1-9. We flag the 29 objects with NUV and 
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FUV detections and z p hot > 1-9, of which we estimate ~18 are actually low-redshift quasars. Most (23/29) of the objects with 
both FUV and NUV detections are already flagged for rejection by the color-color cut, so this process finds 6 new objects with 
questionable photometric redshifts. After combining methods we generate a rejection list of 43 objects with photometric redshifts 
z > 1.9 but GALEX properties more characteristic of low-redshift quasars; we estimate ~(30-35)/43 are indeed low-redshift, 
and only ~10% of genuine z > 1.9 quasars are lost through these cuts. The GALEX detection rate of low-redshift SDSS quasars 
is sufficiently high (Trammell et al. 2007) that virtually all interlopers with GALEX coverage are expected to be identifiable 
through this process. 

We checked the rejection list of 43 objects (with photometric redshifts z > 1.9 but GALEX properties more characteristic 
of low-redshift quasars) for SDSS DR6 or NED spectroscopic coverage, finding 9 objects with spectroscopic redshifts available 
for comparison. From these, 7 objects have incorrect photometric redshifts z > 1.9 but spectroscopic redshifts z < 1 (of 
which 6/5/4 objects were flagged for rejection by color-color /joint-detection/both cuts.) The other 2 objects also have (less 
dramatically) incorrect photometric redshifts, with spectroscopic redshifts ~0.5 below the photometric estimate (these objects 
are flagged for rejection by both methods). No objects with spectroscopically confirmed photometric redshifts of z > 1.9 are on 
this rejection list. We correct the redshifts (and luminosities) of these 9 objects with spectroscopic coverage and then remove 
them from the rejection list, which therefore consists of 34 objects. Of the 30 objects with photometric redshifts z > 1.9 that 
lack GALEX coverage, 5 also have spectroscopic redshifts and are retained, while the remaining 25 (of which ~7 are likely 
actually low-redshift) are conservatively added to the rejection list. After discarding objects on the rejection list, the updated 
candidate list of photometric RlQs and RLQs contains 347 (= 406 — 34 — 25) objects, within which the remaining fraction with 
this type of redshift mismdcntification is only ~1.5%. 

C. NOTES ON INDIVIDUAL DEEP-FIELD OBJECTS 

Here we present brief commentary on selected RlQs and RLQs from the CDF-N and the E-CDF-S; the interested reader is 
referred to the provided references for additional detail. 

123538.51 + 621643.0 is included in the XMM-Newton Bright Serendipitous Survey (Delia Ceca et al. 2004). Galbiati et 
al. (2005) describe it as a non-blazar AGN (with R* — 1.22 and a T = 2.5 calculated from 1.4 to 8.5 GHz) and find that the 
X-ray spectrum can be adequately fit by a power-law model with T — 1.96i 0B and no intrinsic absorption. 

123649.62+620737.8 is near the borderline of our X-ray hardness-ratio cut and may be mildly obscured, although it is 
nevertheless relatively X-ray bright. It is classified as a type-2 AGN by, e.g., Padovani et al. (2004), using the definition of 
Szokoly et al. (2004). Trouille et al. (2008) characterize the optical spectrum as showing high excitation lines (i.e., not an 
obvious broad-line object), but its A(g — i) value is —0.17, bluer than the typical quasar at that redshift. 

123746.63+621739.0 is listed in the radio catalog of Richards et al. (1998) as having a radio flux density of 11.1 fijy at 8.5 
GHz, and is also in the radio catalog of Richards et al. (1999) with a flux density of 998 /iJy at 1.4 GHz, in agreement with 
the flux density from Biggs & Ivison (2006). The optical counterpart has quasar- like colors but lacks an optical spectrum that 
would permit definitive classification. 

033115.03—275518.5 has a lobe-dominated radio morphology and shows double-peaked optical broad-line emission structure. 
Details may be found in Luo et al. (2009) . 

033124-86— 275207 .1 has two other optical objects within 2.5", but the selected counterpart is secure. It is described in 
Rovilos et al. (2007) as a type-1 AGN based on its soft X-ray spectrum. 

033139.49-274119.6 is listed in the ACTA radio catalog of Norris et al. (2006) as having a 1.4 GHz flux density of 0.2 mjy, in 
agreement with the flux density of 206 pJy from Miller et al. (2008). This object is detected in the Spitzer Wide- Area Infrared 
Extragalactic Survey (Lonsdale et al. 2003) and has optical quasar-like colors but lacks an optical spectrum that would permit 
definitive classification. 

033208.66- 2747344 is a bright RIQ with a broad-line optical spectrum (e.g., Tozzi et al. 2009). Lehmer et al. (2005) describe 
the X-ray spectrum (containing ~20000 counts) as having an effective photon index of T = 1.73 (calculated from the hardness 
ratio), and Steffen et al. (2006) quote a radio-loudness of R* = 1.14. 

033210.91-2744U.9 is a RLQ with a broad-line optical spectrum (e.g., Tozzi et al. 2009). Wang et al. (2007) find that the 
X-ray spectrum (containing ~1000 counts) can be well-fit by a power-law model with Y = 1.88T0.09 and no intrinsic absorption. 

033211.63-273725.9 has COMBO-17 colors consistent with those of a template quasar spectrum. Lehmer et al. (2005) 
describe the X-ray spectrum as having an effective photon index of F = 1.44 (calculated from the hardness ratio), and Steffen et 
al. (2006) quote a radio-loudness of R* = 1.29. This source is detected in the Spitzer Wide- Area Infrared Extragalactic Survey. 

033227.01-274105.0 is a broad-line quasar (e.g., Tozzi et al. 2009). Galbiati et al. (2005) quote R* = 2.05 and a r = -0.53 
and fit the X-ray spectrum with a power-law plus blackbody model (r = 2.04 ± 0.15 and kT — 0.15 ± 0.06) with no intrinsic 
absorption. 

D. VERIFICATION OF FITTING METHODOLOGY 

We fit X-ray luminosity as a function of various parameters using the IDL code of Kelly (2007), which takes into account 
censoring of the dependent variable and uncertainties in all variables. This is an advance over methods that consider either 
censoring (e.g., ASURV) or errors (e.g., the IDL program fitexy) but not both simultaneously. Here, we illustrate the fitting 
technique through application to the RQQ sample of Steffen et al. (2006), for which results using an alternative method (ASURV) 
have already been reported. 

In this sample of RQQs the uncertainties are dominated by intrinsic variability (see, e.g., §3.5 of Gibson et al. 2008), as 
is also the case for the RIQs and RLQs that form the primary focus of this work. The uncertainties we use throughout 
are based on observed data, but we investigate here the degree to which the fitting results are sensitive to alternative values 
for the uncertainties. Figure 21a shows 4(C) computed for RQQs for our standard errors (solid line), for a fit resulting 
from doubling the standard errors (dashed line), and for a fit resulting from halving the standard errors (dotted line). The 
coefficients for the standard error fit are 4 = (—0.545 ± 0.023) + (0.649 ± 0.021) x £ uv ; for the double errors fit they are 
4 = (-0.540 ± 0.023) + (0.673 ± 0.022) x £ uv ; for the half errors fit they are 4 = (-0.546 ± 0.023) + (0.644 ± 0.021) x £ uv . 
When comparing the coefficient values it must be kept in mind that these coefficients are probabilistically derived, so a large 
number of fits would be required to quantify most accurately the impact of using different error values. However, it is clear 



33 



30 
29 

28 

27 

26 

25 

24 
30 

29 
28 
27 
26 
25 
24 




Stondord errors 
Double errors 
Holf errors 



(a) 



l l l l | l ¥ | i i i i i i i i i | l 



l l | l 



l l I | i i i i 



t r dep, Kelly 
i uu dep, Kelly 
Bisector, Kelly 
(. dep, ASURV 
t m dep, ASURV 
Bisector. ASURV 




(b) 



' 1 ' " ' <■ ' 



— 1 — i — i — i — i — i — i — i — i — I — i — i — i — i — i — i — i—i — i — I — i — i — i — i—i — i — i- 



28 



29 



30 



31 



32 



uv 

Figure 21. Demonstration of fitting methodology for f x (f m ) applied to RQQs. The top panel illustrates that varying the 
uncertainties in the variables within a reasonable range has only a minimal impact upon the best-fit correlation. The bottom 
panel shows that using the code of Kelly (2007) provides results consistent with those given by ASURV, and illustrates the effect 
of considering £ uv rather than £ x as the dependent variable, or of using the bisector line. See Appendix D for details. 



that varying the uncertainties within reasonable bounds does not dramatically alter the parameters of the best-fit line (i.e., the 
credible intervals of the parameters overlap). 

The best-fit model obtained through using the Kelly (2007) code agrees with that calculated by Steffen et al. (2006) using 
ASURV and treating £ x as the dependent variable (their Equation fa), which is £ K = (-0.546±0.023) + (0.642±0.021) x£ uv using 
our normalization convention. As described in §5.1, it seems reasonable on both practical and statistical grounds to consider £* 
as the dependent variable throughout. However, other approaches are possible; we investigate the effect of considering £ uv as the 
dependent variable and of using the bisector best-fit line. When fitting £ U v(£x) and inverting the results, the slope will generally 
be steeper than for the fit conducted with £ x as the dependent variable. A complicating factor is that it is no longer possible to 
consider the censoring in £ x . However, treating all upper limits as X-ray detections, we find £ x = —0.475 + 0.812 x £ uv , which 
may be compared with the Steffen et al. (2006) result (their Equation lb) of £ x = —0.400 + 0.815 x £ uv . For reference, the best-fit 
lines calculated using the IDL program fitexy are £ x = —0.500 + 0.695 x £ uv and £ x — —0.491 + 0.736 x £ uv treating £ K and £ uv 
in turn as the dependent variable. Calculating the bisector using Table 1 of fsobe et al. (1990) gives £ x = —0.493 + 0.727 x £ uv 
for our fits, or £ K — —0.478 + 0.721 x £ uv from Steffen et al. (2006; their Equation lc). The bisector slope is slightly higher, but 
it may be seen in Figure 21b that the change in the distance from the best-fit line for any given point is at most ~0.2 over the 
span of the variables. When we discuss "outliers" or "positive residuals" we are referring to greater distances from the best-fit 
line. We note that the best-fit lines may not always appear to run through the highest density of points when upper limits are 
present; for example, the RIQ fit to £^(£ U v) in Figure 9a may appear different than would be the case were all of the upper 
limits detections. 
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As a final consistency check, we also performed every fit given in Table 7 with ASURV0 In all cases, the parameters agreed (in 
the sense that the la intervals overlapped) with those computed using the code of Kelly (2007). The median differences between 
the best-fit parameters obtained using the code of Kelly (2007) and those obtained with ASURV (and the standard deviations 
thereof) for the £^(£ U v) model are —0.0009 (0.0085) for clq an d 0.0262 (0.0082) for b uv ; without considering the uncertainties, 
ASURV apparently computes slightly flatter slopes for these datasets (recall that the typical values for b uv are 0.6-0.8 in these 
fits, so this systematic difference is 3-4%). The median differences (and standard deviations) for the £ x (£ U v,£r) model are 
0.0006 (0.0049) for a , 0.0200 (0.0130) for 6 UV , and -0.0088 (0.0094) for c r . For illustrative purposes, if the Kelly (2007) code 
is run with an uncertainty on £ uv of 0.05 rather than 0.114 (keeping the uncertainties on £ T and £^ at 0.0792 and 0.146), the 
average change in 6 uv /c r is —0.022/0.0079, producing coefficients matching the ASURV results (with a similar impact on 6 UV in 
the 4(4v) relation). 

We conclude that using only £^ as the dependent variable is acceptable for our comparative study and that using the IDL 
code of Kelly (2007) to take both censoring and uncertainties into account is advantageous. 



39 For computational purposes, for the full, primary, off-axis, , r , , , cr , n ... A CT Tmr 

i tit ^ i i . i i ■ i i , ■ r^r* \ ■ i \ randomly selected ~50(J objects to lit with AbUKV. 

and KLQ subsamples only (which each contain >500 objects) we 
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Table 1 

Primary sample RIQs and RLQs 



Name a 


z 


m\ 


A( 9 - i) 








Dct c 


R* 






^r,corc 


SDSS° 


X-ray f 


000442.18+000023.3 


1.008 


19.02 


-0.02 


31.92 


30.28 


26.91 


1 


1.64 


-1.30 


0.00 


31.92 


S 


X 


000622.60-000424.4 


1.038 


19.49 


0.21 


34.97 


30.05 


27.29 


1 


4.91 


-1.06 


-0.72 


34.97 


s 


c 


001130.40+005751.8 


1.492 


20.08 


-0.05 


33.96 


30.13 


27.17 


1 


3.84 


-1.13 


-0.22 


33.96 


s 


X 


004230.40-092202.7 


2.085 


19.97 


0.01 


34.01 


30.51 


26.31 


1 


3.50 


-1.61 


0.00 


34.01 


p 


c 


004413.72+005141.0 


0.941 


18.44 


0.06 


32.64 


30.42 


26.61 


1 


2.22 


-1.46 


0.00 


<31.29 


s 


X 


005009.81-003900.6 


0.728 


19.89 


0.63 


31.68 


29.36 


26.25 


1 


2.32 


-1.20 


0.00 


31.68 


s 


C 


005905.50+000651.6 


0.719 


17.46 


-0.07 


34.42 


30.62 


27.37 


1 


3.80 


-1.24 


-0.40 


34.42 


s 


c 


012401.76+003500.9 


1.850 


20.23 


0.26 


34.12 


30.26 


26.90 


1 


3.86 


-1.29 


-0.70 


34.12 


s 


c 


012528.84-000555.9 


1.077 


16.47 


-0.09 


34.59 


31.36 


27.73 


1 


3.23 


-1.39 


0.22 


34.59 


s 


r 


012734.57-000523.8 


1.598 


20.37 


0.14 


33.32 


30.04 


26.74 


1 


3.28 


-1.27 


0.00 


33.32 


s 


r 



Note. — Tabic 1 is provided in its entirety in the online edition of the journal. A portion is shown here for guidance as to its form and content. Some 
values here and in the following tables are given to precision greater than the uncertainty to avoid future accumulated round-off error. 
a Name is J2000 from optical RA/Dec and is SDSS DR5 for spectroscopic quasars or SDSS DR6 for photometric quasars. 

k The radio luminosity £ r is calculated from both core and lobe components as described in §3.2, while £ riCorc is for the core only (limits indicate no FIRST 

core and arc calculated for 1 mjy). 

c 1 — X-ray detection; — X-ray upper limit. 

a r is calculated between 1.4 and 5 GHz where Green Bank fluxes arc available and is otherwise based on low-frequency (Texas and Westerbork) surveys; 
zero-valued entries indicate a lack of multi-frequency radio data. 

° S — spectroscopic, from the DR5 quasar catalog of Schneider ct al. (2007); P — photometric, from the DR6 quasar catalog of Richards ct al. (2009). 
C/c — Chandra: X/x — XMM-Newton; r — ROSAT. Lowercase indicates our measurements, uppercase from catalogs as described in §3.4. 
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Table 2 

RIQs and RLQs selected from deep surveys 



Name 1 


z 










R* 






ReF 




Color selected 


033115.03-275518.5 


1.368 


21.43 


32.71 


29.60 


26.48 


3.11 




-1.20 


1 


033124.86-275207.1 


1.328 


21.09 


33.18 


29.31 


25.72 


3.87 




-1.38 




033139.49-274119.6 


2.215 


23.12 


31.36 


29.29 


<25.04 


2.06 


< 


-1.63 




033208.66-274734.4 


0.543 


18.78 


31.03 


29.26 


25.94 


1.77 




-1.27 


2 


033210.91-274414.9 


1.605 


22.95 


32.26 


29.10 


26.08 


3.16 




-1.16 


3 


033211.63-273725.9 


1.636 


18.99 


32.36 


30.72 


26.85 


1.64 




-1.49 




033227.01-274105.0 


0.737 


19.09 


32.31 


29.95 


26.34 


2.36 




-1.39 


3 


033302.67-274823.1 


3.021 


23.96 


31.38 


29.31 


26.51 


2.06 




-1.08 




033310.20-274841.9 


1.034 


22.66 


32.68 


28.58 


25.71 


4.10 




-1.10 


3 


123529.38+621256.4 


2.413 


23.21 


31.24 


29.37 


26.11 


1.88 




-1.25 


5 


123538.51+621643.0 


0.712 


19.70 


31.50 


29.53 


26.50 


1.97 




-1.16 


5 


123649.62+620737.8 


1.610 


23.51 


31.25 


28.85 


26.17 


2.40 




-1.03 


6 


123746.63+621739.0 


2.316 


23.36 


32.12 


29.25 


<25.20 


2.87 


< 


-1.56 




095958.54+015254.6 


1.019 


21.86 


31.02 


28.56 


24.80 


2.46 




-1.44 




100046.92+020726.6 


1.210 


21.70 


31.81 


28.97 


25.37 


2.85 




-1.38 




100114.86+020208.9 


0.969 


20.63 


31.97 


29.31 


26.00 


2.66 




-1.27 


8 


Broad-line selected 


033225.17-274218.8 


1.617 


23.37 


30.70 


28.82 


25.78 


1.87 




-1.17 


4 


123704.11+620755.4 


1.253 


22.07 


30.44 


29.11 


25.54 


1.32 




-1.37 


7 


123707.51+622148.0 


1.451 


23.04 


30.45 


28.92 


26.20 


1.53 




-1.04 


7 


095821.65+024628.2 


1.403 


19.05 


32.36 


30.48 


26.75 


1.88 




-1.43 


8 


095838.47+022439.3 


1.161 


21.86 


30.34 


28.71 


26.23 


1.63 




-0.95 


8 


095908.32+024309.6 


1.317 


18.43 


33.32 


30.66 


27.16 


2.66 




-1.35 


8 


095921.31+024412.4 


1.004 


20.46 


31.57 


29.42 


25.53 


2.15 




-1.49 


8 


100114.86+020208.8 


0.969 


21.00 


31.97 


29.27 


26.18 


2.70 




-1.19 


8 


100129.83+023239.0 


0.826 


20.54 


30.36 


28.80 


25.54 


1.55 




-1.25 


8 


100205.03+023731.5 


0.519 


19.05 


30.68 


29.25 


26.23 


1.43 




-1.16 


8 


100213.42+023351.7 


1.143 


21.76 


30.25 


28.94 


25.68 


1.32 




-1.25 


8 


100228.82+024016.9 


3.144 


21.27 


31.56 


29.69 


26.76 


1.87 




-1.12 


8 


100230.06+014810.4 


0.626 


19.65 


30.22 


28.30 


25.51 


1.92 




-1.07 


8 


100240.93+023448.4 


1.677 


21.98 


30.74 


29.35 


25.89 


1.40 




-1.33 


8 


100249.33+023746.5 


2.124 


19.75 


32.84 


30.62 


26.37 


2.22 




-1.63 


8 



a Name is .T2000 from optical RA/Doc. 

b X-ray luminosities are from Chandra observations except for the broad-line selected COSMOS objects, 
for which the X-ray luminosities are from XMM-Newton observations. See §2.2.3 for details. 
c References for spectroscopic rcdshifts: 1 — Luo ct al. (2009); 2 — Mignoli et al. (2005); 3 — Silverman 
ct al., in prep; 4 = Szokoly ct al. (2004); 5 = Barger ct al. (2003); 6 = Bargcr et al. (2008); 7 = Trouille 
ct al. (2008); 8 = Trump ct al. (2009). Unlablcd redshift values are photometric from Luo ct al. (2010; 
CDF-S X-ray detections); Rafferty ct al. (2010; X-ray limits and E-CDF-S detections); Xue et al. (2010; 
CDF-N); and Ilbert ct al. (2009; COSMOS). 
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Table 3 

Sample characteristics 



X-ray sources Median properties 

Sample N N dct 100 z t r t uv £x R* a ox 



Full, primary, and supplemental samples 



Full 


791 


674 


85 


1 


.40 


32 


.95 


30. 


.42 


26 


.79 


2 


.59 


-1 


.40 


Primary 


654 


547 


84 


1 


.41 


32. 


.83 


30 


36 


26 


73 


2. 


.50 


-1 


.40 


Spectroscopic 


312 


274 


88 


1 


.32 


33. 


.00 


30. 


49 


26 


83 


2 


.46 


-1 


.41 


QSO/HIZ 


200 


186 


93 


1 


.20 


33 


.14 


30 


67 


26 


97 


2 


.50 


-1 


.41 


Photometric 


342 


273 


80 


1 


.47 


32 


.72 


30 


23 


26 


64 


2 


.52 


-1 


.39 


Off-axis 


562 


456 


81 


1 


.45 


32 


73 


30 


33 


26 


.70 


2 


.44 


-1 


.40 


Targeted 


92 


91 


99 


1 


.05 


33 


65 


30 


.74 


27 


04 


3 


.09 


-1 


.39 


Supplemental 


137 


127 


93 


1 


.32 


34 


33 


30 


96 


27. 


.38 


3 


11 


-1 


.36 


Einstein 


93 


85 


91 


1 


.02 


34 


63 


31 


15 


27. 


.55 


3 


.44 


-1 


.39 


High-z 


13 


13 


100 


4 


.31 


34 


69 


31 


42 


27 


70 


2 


.90 


-1 


.34 


Deep fields 


31 


29 


94 


1 


.32 


31 


.50 


29. 


27 


26 


.04 


2 


06 


-1 


.26 


FIRST 


180 


164 


91 


1 


.27 


33. 


16 


30. 


73 


26 


99 


2 


.41 


-1 


.43 



PJQs and RLQs 



RIQs 


188 


142 


76 


1.42 


32.05 


30.41 


26.50 


1.67 


-1.50 


RLQs 


603 


532 


88 


1.38 


33.30 


30.44 


26.89 


2.88 


-1.37 


Groupings of RLQs 


R* < 3 


343 


291 


85 


1.42 


32.87 


30.34 


26.69 


2.47 


-1.41 


R* > 3 


260 


241 


93 


1.33 


34.08 


30.59 


27.16 


3.46 


-1.31 


£ T < 33.3 


301 


251 


83 


1.26 


32.68 


30.11 


26.53 


2.50 


-1.37 


£ r > 33.3 


302 


281 


93 


1.63 


34.04 


30.82 


27.22 


3.31 


-1.37 


a r < -0.5 


181 


160 


88 


1.31 


33.59 


30.48 


26.93 


3.18 


-1.38 


a r > -0.5 


221 


207 


94 


1.26 


33.66 


30.61 


27.11 


3.10 


-1.34 



Note. — The columns N, N^ ct , and 100 — give the number of sources, the number of X-ray detected sources, 
and the percentage of sources with X-ray detections, respectively. Radio, optical/UV, and X-ray monochromatic 
luminosities £ r , £ uv , and £ x are expressed as logarithmic quantities with units of erg s -1 Hz -1 at rest-frame 5 GHz, 
2500 A, and 2 kcV, respectively. The radio loudness is R* — £ r — £ uv and the optical/UV-to- X-ray spectral slope 
is a ox — 0.384 x (£ x — £ uv ). RIQs have 1<R*< 2 and RLQs have R* > 2. Medians that include limits were 
determined with ASURV. Details regarding the various samples arc given in §2. 
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Table 4 

RIQ and RLQ properties as a function of I, 



Statistic 


<29.5 


29.5-30 


£ nv bin 
30-30.5 


range 
30.5-31 


31-31.5 


>31.5 








X- 


■ray sources 








N 

N dct 
100 


69 
57 
83 


119 

98 
82 


257 
217 
84 


210 

176 
84 


92 
86 
93 


44 
40 
91 


Rcdshift 


Mean 
25th % 
50th % 
75th % 


0.87 
0.47 
0.71 
1.15 


1.04 
0.74 
1.03 
1.33 


1.36 
1.05 
1.39 
1.69 


1.72 
1.09 
1.70 
2.20 


2.09 
1.27 
1.99 
2.90 


2.71 
2.01 
2.44 
3.52 


^uv 


Mean 
25th % 
50th % 
75th % 


29.12 
28.85 
29.27 
29.38 


29.80 
29.67 
29.82 
29.92 


30.27 
30.16 
30.28 
30.40 


30.73 
30.61 
30.73 
30.85 


31.21 
31.11 
31.19 
31.31 


31.78 
31.60 
31.73 
31.91 




Mean 
25th % 
50th % 
75th % 


25.77 
25.38 
25.75 
26.23 


26.35 a 
26.06 
26.38 
26.57 


26.67 
26.42 
26.72 
26.94 


26.97 
26.72 
27.03 
27.21 


27.52 
27.27 
27.53 
27.78 


27.83 
27.48 
27.71 
28.13 




Mean 
25th % 
50th % 
75th % 


-1.28 
-1.46 
-1.25 
-1.16 


-1.32 
-1.43 
-1.32 
-1.25 


-1.39 
-1.46 
-1.38 
-1.28 


-1.45 
-1.52 
-1.44 
-1.34 


-1.42 
-1.53 
-1.42 
-1.31 


-1.52 
-1.62 
-1.54 
-1.45 




Mean 
25th % 
50th % 
75th % 


2.66 
1.92 
2.66 
3.27 


2.52 
2.03 
2.43 
2.90 


2.51 
1.97 
2.47 
3.01 


2.67 
2.06 
2.64 
3.37 


2.96 
2.33 
3.14 
3.71 


2.67 
2.19 
2.77 
3.24 



Note. — Properties here arc for the full sample of RIQs and RLQs. Quantities and 
units arc as defined in Table 3. The £ uv bins are plotted in Figure 3. 



a Estimate of the mean is biased because it was computed by ASURV treating the censored 
first point in this bin as a detection. 
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Table 5 









Excess X-ray luminosity 


as a function of R* and t T 


: full sample 








Statistic 


<2.0 


2.0-2.5 


R* bin range 
2.5-3.0 


3.0-3.5 


> 3.5 


< 32.2 


£ r bin range 
32.2-32.9 32.9-33.6 


33.6-34.3 


> 34.3 


N 

N dct 
100 %t 


188 
142 
76 


175 
145 
83 


168 
146 
87 


141 
128 
91 


119 
113 

95 


176 
138 
78 


198 
156 
79 


187 
163 
87 


108 
102 
94 


122 
115 
94 


^x - ^x.RQQ 


Mean 
Error 
25th % 
50th % 
75th % 


0.11 
0.04 
-0.16 
0.16 
0.45 


0.29 
0.03 
0.08 
0.30 
0.52 


0.41 
0.03 
0.26 
0.44 
0.62 


0.58 
0.03 
0.35 
0.57 
0.77 


0.77 
0.05 
0.53 
0.71 
1.03 


0.22 
0.04 
-0.13 
0.29 
0.59 


0.22 
0.04 
0.05 
0.29 
0.52 


0.42 
0.03 
0.22 
0.41 
0.62 


0.55 
0.04 
0.34 
0.53 
0.72 


0.76 
0.04 
0.51 
0.67 
0.94 


R* 


Mean 
Error 
25th % 
50th % 
75th % 


1.62 
0.02 
1.42 
1.67 
1.86 


2.26 
0.01 
2.14 
2.28 
2.39 


2.75 
0.01 
2.62 
2.74 
2.88 


3.25 
0.01 
3.11 
3.25 
3.38 


3.85 
0.03 
3.62 
3.78 
4.00 


1.88 
0.04 
1.52 
1.86 
2.14 


2.27 
0.04 
1.88 
2.29 
2.58 


2.77 
0.04 
2.50 
2.79 
3.07 


3.15 
0.05 
2.88 
3.20 
3.48 


3.61 
0.04 
3.35 
3.65 
3.87 


It 


Mean 
Error 
25th % 
50th % 
75th % 


31.98 
0.05 
31.56 
32.05 
32.36 


32.62 
0.05 
32.17 
32.58 
32.95 


33.14 
0.05 
32.79 
33.18 
33.47 


33.75 
0.06 
33.34 
33.70 
34.23 


34.38 
0.07 
33.93 
34.51 
34.94 


31.69 
0.03 
31.44 
31.82 
32.03 


32.54 
0.01 
32.36 
32.54 
32.70 


33.23 
0.01 
33.07 
33.23 
33.39 


33.91 
0.02 
33.72 
33.92 
34.06 


34.77 
0.03 
34.46 
34.72 
34.97 



Note. — The "excess" X-ray luminosity is defined as £ x — ^ x ,rqq ? where ^ X) rqq — 0.709^ uv + 4.822 from Just et al. (2007). Other quantities 
and units are as defined in Table 3. Note that the first R* bin is the sample of RIQs. The inner bins in R* increase by linear factors of ~3; the 
inner bins in l r increase by linear factors of ~5. These values arc plotted in Figure 7. 



Table 6 

Excess X-ray luminosity as a function of R* and £ r : primary sample 



Statistic 


<2.0 


2.0-2.5 


R* bin range 
2.5-3.0 


3.0-3.5 


> 3.5 


< 32.2 


t T bin 
32.2-32.9 


range 

32.9-33.6 


33.6-34.3 


> 34.3 


N 

N dct 


171 
125 


154 
126 


147 
128 


111 
102 


71 

66 


154 
118 


187 
145 


172 
150 


90 
84 


51 
50 


100 %t 


73 


82 


87 


92 


93 


77 


78 


87 


93 


98 


^x - ^x.RQQ 


Mean 
Error 
25th % 
50th % 
75th % 


0.07 
0.04 
-0.20 
0.13 
0.41 


0.30 
0.03 
0.12 
0.30 
0.52 


0.42 
0.03 
0.26 
0.44 
0.62 


0.56 
0.03 
0.34 
0.54 
0.74 


0.67 
0.06 
0.53 
0.65 
0.97 


0.20 
0.04 
-0.13 
0.23 
0.57 


0.21 
0.04 
0.03 
0.29 
0.52 


0.44 
0.03 
0.23 
0.41 
0.63 


0.54 
0.04 
0.34 
0.53 
0.71 


0.70 
0.04 
0.52 
0.65 
0.87 


R* 


Mean 
Error 
25th % 
50th % 
75th % 


1.62 
0.02 
1.41 
1.67 
1.84 


2.26 
0.01 
2.14 
2.28 
2.39 


2.74 
0.01 
2.61 
2.74 
2.87 


3.25 
0.01 
3.12 
3.25 
3.38 


3.82 
0.03 
3.61 
3.77 
3.92 


1.86 
0.04 
1.51 
1.84 
2.11 


2.25 
0.04 
1.87 
2.28 
2.56 


2.77 
0.04 
2.50 
2.80 
3.07 


3.21 
0.05 
2.95 
3.27 
3.49 


3.66 
0.05 
3.45 
3.70 
3.85 




Mean 
Error 
25th % 
50th % 
75th % 


32.06 
0.04 
31.69 
32.07 
32.38 


32.57 
0.04 
32.17 
32.54 
32.92 


33.08 
0.04 
32.76 
33.13 
33.41 


33.58 
0.06 
33.26 
33.61 
34.04 


34.09 
0.09 
33.66 
34.32 
34.67 


31.78 
0.03 
31.55 
31.87 
32.04 


32.54 
0.01 
32.36 
32.52 
32.70 


33.23 
0.02 
33.07 
33.22 
33.39 


33.89 
0.02 
33.70 
33.90 
34.04 


34.61 
0.04 
34.38 
34.57 
34.74 



Note. - 



Quantities defined as in Table 5. 
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Table 7 

Correlations with X-ray luminosity 



= do + b uv x £ uv £ x = ao + b uv x £ uv 4-c r x4 

Sample cto ^uv &o b uv c r 



Full, primary, and supplemental samples 



Full 


-0. 


144+0-017 
14 -0.016 


0. 


7aQ +0.026 
<° y -0.025 


-0 


nq „+0.014 
U9 °-0.015 


0. 


4H9 +0.031 
4 °^-0.031 


0. 


97 o+0.019 

■ a ' ,J -o.oi9 


Primary 


-0. 


lqo +0.018 
±yu -0.018 


0. 


791 +0.031 
- ' z± -0.031 


-0 


116 +0.016 
1±D -0.016 


0. 


ion+0.033 
48U -0.034 


0. 


nco+0.021 
- zo -0.021 


Spectroscopic 


-0. 


1Q o+0.024 
lyo -0.024 





640+0 042 
o4u -0.043 


-0 


1n ,+0.021 
luo -0.020 





oon+0 041 

•J° u _0.042 





9Qfi +0.025 
zyu -0.024 


QSO/HIZ 


-0. 


ico+0.030 
luo -0.030 


0. 


Rri+0.049 
DJ± -0.048 


-0 


na o+0.025 
U °°-0.025 


0. 


400+ 048 
4uu -0.047 





OO4+0-030 
zy -0.031 


Photometric 


-0. 


i R9 + 0.027 

1DZ -0.028 


0. 


s1fi +0.050 
Olu -0.048 


-0 


,11+0.026 
1±1 -0.027 


0. 


fi1 1+0.057 
D1± -0.058 





91 ,+0.035 
- zi -0.034 


Off-axis 


-0. 


9n o+0.020 
zuo -0.020 


0. 


744+0-036 
' 4 +-0.036 


-0 


1 QQ + 0.019 

iz,3 -0.020 


0. 


crno+0.039 
Jzo -0.040 


0. 


9 io+0.023 
z4J -0.025 


Targeted 


-0. 


n o q +0.045 
UJy -0.043 





c-co+0.064 
OOJ -0.065 


-0 


nsQ +0.035 
uoy -0.037 





ooc+0.059 

O,:>0 -o.o59 





o 9n +0.044 
JZU -0.044 


Supplemental 


+0. 


n9 o+0.044 
UZO -0.044 





ooQ+0 045 
5ZO -0.044 


-0 


017 +0.041 
u± ' -0.041 





49 i+0.079 
■ 4y± -0.079 


0. 


9 ko+0.051 
z O^_0.052 


Einstein 


+0. 


og+0.072 
uoy -0.075 





7Q o+0.087 
,yo -0.087 


-0 


9OQ+0.072 





479 +0.089 
■ 4,z -0.090 





4Q7 +0.066 
4U -0.068 


FIRST 


-0. 


iOfj+0.035 
iyD -0.034 


0. 


bt)8 -0.063 


-0 


nQS +0.028 
■ UyS -0.030 


0. 


4nq +0.060 
+ uy -0.059 





9Q 7+0.032 
zy ' -0.032 



RQQs, RlQs, and RLQs 



RQQs 

RfQs 

RLQs 



-0 545+0- 023 

-0 447+ 037 
-0.038 

-0 057+ ' 017 



0.6491°;™ 
0.567tg;»? 
0.8311°;°^ 



-0.477 
-0.100 



+0.212 
0.215 
0.016 



0.5881°;^ 
0.506+°°?! 



-0.025l°;i 77 
0.2921°;™ 



Groupings of RLQs 



R* < 3 


n ,qi +0.020 
-U.191_ 020 


0. 


771 +0.032 
' '1-0.032 


-0 149+0- 032 
u.ria_ 32 





620 +o. 100 
uzu -0.100 


0. 


144+0.088 
144 -0.091 


R* > 3 


+0.112t°;°|J 


0. 


o 44 +0.036 
- 04 -0.034 


n 177 +0.057 
U - i "-0.058 


0. 


4 o 9 +0.081 
■ 4oz -0.085 


0. 


nqo+0.073 
oyo -0.070 


tr < 33.3 


-0.180i°;°i 





7 o Q +0.045 
loy -0.046 


-0 098+ 040 
u.uac5_ 040 


0. 


61Q +0.061 
u -0.062 


0. 


9n r:+0.065 
ZUO -0.065 


e T > 33.3 


+0 r)70+°'° 28 
i-U.U(U_ Q 028 


0. 


fis7 +0.045 
°S<_0.044 


-0 170+ 034 
u - 1,u -0.035 


0. 


4 oo+0.047 
■^ oy -0.048 


0. 


oon+0.040 
ooy -0.040 


a r < -0.5 


-o.09it°;^ 





fi c 9 +0.045 
DOZ -0.045 


-0 175+0-025 
u -l'°-0.025 


0. 


000 + 0.061 

JDJ -0.061 





25O+0 041 
zoy -0.041 


a r > -0.5 


+0.066t°;°l 


0. 


snn+o° 39 

8UU -0.040 


-0 041 +°' 026 
U.U41_ 026 


0. 


471 +0.052 
+ ' i -0.050 


0. 


o,i 9 +o.o39 

• 34z -0.040 



Note. — Fitting is done with the IDL code of Kelly (2007), which utilizes Baycsian techniques that incorporate 
both errors and upper limits. Errors are assumed to be dominated by intrinsic random flux variability; sec, e.g., §3.5 of 
Gibson ct al. (2008). The luminosities are normalized prior to fitting as £ r — 33.3, £ uv — 30.5, £ x — 27.0. For any given 
model fit the quoted parameter values are the median of draws from the posterior distribution and the errors are la. 
These results arc plotted in Figures 9-13. 
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Table 8 

a ox and Aa ox as functions of R* 



Sample N z a R* a aj ox /Aa ox a Intercept a Slope b 



«ox = a + b X R* 



2 < 1 

1< z < 2 
z > 2 


213 
406 
172 


0.69 
1.44 
2.49 


2.68 
2.54 
2.56 


-1.33 
-1.39 
-1.49 


-1 556+ ' 043 
i.OOO_ .o 43 

i ccQ+0.029 
l.DOO_ 02g 

1 SR9+ - 048 
-i.»b2_ 048 


79+0.0l5 

101+ 0011 
u - lui -0 .011 

152+ 017 
U - 1JZ -0.017 








Aa ox = a - 








2 < 1 

1 < z < 2 
z>2 


213 
406 
172 


0.69 
1.44 
2.49 


2.68 
2.54 
2.56 


0.18 
0.17 
0.16 


_n mi +0.040 

U.1U1_q Q41 
n n „+0.027 

-u.±is_ 027 

-0 234+ 047 
u.zoi_ 04g 


097+ ' 014 
u - uy ' -0.014 

106+ 010 
u.iuo_ 010 

n ic;c+0.017 
u.ioo_ 017 



Note. — Fitting as described for Table 7. Uncertainties for the fitting were propagated from 
those assumed for luminosities. 

a Median values from within the subsample. The first three rows give a ox values and the last three 
rows give Aq ox values, where Act ox corresponds to a ox — a x,RQQ with q ox ,rqq — —0.140 X £ uv + 
2.705 (Equation 3 from Just ct al. 2007; sec also §4). 



Table 9 

Description of model components 



Parameter 


Value a 


Commcnt b 


z 




Quasi-random 


From W01 RLF, then use values close to primary sample to account for rcdshift selection effects. 


e 




Random 


Drawn from uniform distribution in sin 6. 


^min 




5.8°t 


Minimum permitted inclination; values between 4° — 7° are viable. 


Radio 


P, <>>, 7 




-5.61, 1.38, 10.49 


From Table 5 of MH09. Both p and <r p are expressed as natural logs. 


^r, lobes 




Random 


From W01 RLF, then shifted to 5 GHz using a T = -0.9 


^rjetO 




Calculated 


jcto from £ r i b es and p with scatter from u p . 






Calculated 


4jet = 4,jeto'x S 2 - a ', 5 = 7"! (1 - f3 cosd)- 1 , a r = -0.3. 


/r,min 




1, 2.7 mjy 


Flux density limit at 1.4 GHz for core (FIRST 5<r), effective limit for lobes ( < 5% of lobes below). 


Optical 


PdjO; 0" d j() 




-1.8t,0.5 


Logarithmic ratio of disk to intrinsic (unbeamed) radio jet; correlation motivated by W99 results. 


^uv,disk 




Calculated 


^uv.disk = 4,jeto + PdjO with scatter from o- d j 






Calculated 


Extrapolated from £ r jct using a TO = —0.8. Observations restrict optical jet flux. 


^i,max 




19.95/20.13, 


Cutoff mi for z < 1/z > 1; from depth (<5% objects above) of spectroscopic, photometric samples. 






20.60/21.08 




X-ray 


■^x, corona 




Calculated 


Determined from £ uv .disk following the RQQ relation of J07. 


Pj,A, Pj,B, 


Pj.C 


-5.35t,-7.25t,-7.25t 


Logarithmic ratio of X-ray jet-linked emission to intrinsic radio jet for models A, B, and C. 


CTdx, crj 


0.3, 0.7t/0.3t 


Scatter in £ x , corona (from RQQ < x (< nv ) relation) and in pj for z < 1/z > 1 for models A, B, and C. 


7x,A, 7x,B 


. 7x,C 


1.0, 10.49, 2.0 1 " 


Bulk Lorentz factor governing beaming of X-ray jet-linked emission for models A, B, and C. 


^x.jct 




Calculated 


Determined from £ r ,jcto and pj, with boosting of <5 x ~ 2c * x , a x = —0.3/ — 0.5 for B/C. 



Note. — Sec §6 for modeling details. 

a Only those numbers indicated with a dagger symbol were varied to best match the simulated to the observed data; all other numbers were taken from 
indicated references or fixed by observations. 

b References here arc W01: Willott ct al. (2001); MH09: Mullin & Hardcastlc 2009; W99: Willott ct al. (1999); J07: Just ct al. (2007). 
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Table 10 

Model X-ray luminosities and <,(<„», £ r ) 





X-ray 


luminosities 


£ x = ao ■ 


4" t^uv X f uv + 


c r X l T 




Sample 


^x 


a 


KS p a 






Cr 


Primary 


26.80 


0.54 




-0 Hfi+ 016 

U.11D_q 01 g 


U.480_ 034 


0.263 


+0.021 
-0.021 


Model A 


26.74 


0.58 


0.37 


-0.1421H1? 


n 7 o n +0.030 


0.100 


+0.022 
-0.022 


Model B 


26.72 


0.75 


< 0.01 


+0.099l°° 2 2 ° 


105+0035 
u.rao_ 035 


0.746 


+0.027 
-0.027 


Model C 


26.75 


0.55 


0.30 


-o.o98i8:8i? 


526+ ' 028 

U.OZD_q 027 


0.279 


+0.021 
-0.021 


Model CI. 5 


26.78 


0.55 


0.76 


-o.io5i8;8^ 


503+0.028 
u.o»o_ 02g 


0.210 


+0.022 
-0.022 


Model C3.0 


26.78 


0.59 


0.10 


-o.oool8;8i 6 7 


400+ 028 


0.435 


+0.022 
-0.022 



Note. — Models are described in §6. Fitting as described in Table 7. Upper limits are included in the 
simulated sample at a rate consistent with those observed in the primary sample. 

a Kolmogorov-Smirnov test comparing distribution of model £ x values to that of the primary sample 
(computed treating observed and simulated upper limits as detections, as are the median £ x values here). 
The probability p-valucs given indicate that for all but Model B the simulated X-ray luminosities are not 
inconsistent with the observed X-ray luminosities. 



